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A METHOD  FOR  CALCULATING  THE  NATURAL  FREQUENCIES 
OF  CONTINUOUS  BEAMS r FRAMES.  AND  CERTAIN 
TYPES  OF  PLATES 


I.  INTRODUCTION 


it  Object  and  Scope  of  Investigation. 

The  purpose  of  this  report  is  to  present  & method  for  calculating 
the  undamped  natural  frequencies  of  fl crural  vibration  of  elastic  structures. 
The  method  is  applicable  to  continuous  beams  on  rigid  or  flexible  supports, 
to  rigid  jointed  plane  frameworks,  and  to  certain  types  of  continuous  plates. 
The  beams  may  have  any  number  of  spans  of  arbitrary  length  and  any  condition 
of  restraint  at  the  far  ends.  In  general,  the  frames  are  assumed  to  be  fix- 
ed against  lateral  displacement,  but  consideration  is  also  given  to  symmetri- 
cal, single-bay,  multi-story  frames  free  to  undergo  sidesway.  The  plates  are 
assumed  to  be  simply  supported  along  two  opposite  edges  and,  in  one  direction, 
continuous  over  a series  of  rigid  supports  transverse  to  the  simply  supported 
edges.  The  mass  of  the  members  composing  the  structure  is  assumed  to  be  uni- 
formly distributed  along  each  member.  A system  which  has  distributed  mass 
and  elasticity  has  an  infinite  number  of  natural  frequencies.  With  the  method 
presented  herein  one  is  capable  of  determining  all  natural  frequencies  as  well 
as  the  corresponding  natural  modes  of  vibration  of  a system.  The  assumptions 
made  in  the  analysis  are  those  of  the  ordinary  theory  of  flexure  of  beams  and 
of  medium-thick  plates. 

Knowledge  of  the  natural  frequencies  cf  structures  is  Important  for 
the  analysis  and  design  of  structures  subjected  to  time-dependent  forces. 
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This  knowledge  ie  particularly  significant  in  the  case  of  stationary  periodic 
forces  such  as  those  resulting  from  rotating  machinery.  If  the  operating  fre» 
quency  of  the  machinery  is  sufficiently  close  to  one  of  the  natural  frequencies 
of  the  structure  supporting  it,  violent  vibrations  will  ensue  which,  in  the 
absence  of  dissipative  forces,  may  attain  extremely  large  amplitudes.  In 
order  to  avoid,  by  proper  design,  the  destructive  condition  of  resonance,  it 
is  necessary  to  have  a workable  method  for  predicting  the  natural  frequencies 
of  structures.  It-  has  been  the  object  of  this  investigation  t©  attempt  to 
meet  this  need. 

The  problem  of  calculating  natural  frequencies  of  dynamic  systems 
has  been  ths  subject  of  discussion  for  a long  period  of  time.  The  natural 
frequencies  of  single  span  members  having  different  boundary  conditions  have 
been  investigated  rather  exhaustively;  yet,  comparatively  little  has  been 
done  for  multiple  member  systems.  Except  for  structures  consisting  of  only 
a few  members,  the  classical  method  of  determining  natural  frequencies  be- 
comes so  cumbersome  that  it  tends  to  be  entirely  useless  for  practical  pur- 
poses. Several  considerably  more  efficient  methods  have  been  developed,  but 
these  seem  to  be  applicable  to  limited  typ-as  of  structures. 

Among  the  available  methods,  Mudrsk's  method  (1),  (2),  (3)  which 
utilizes  the  effective  stiffness  criterion  for  determining  natural  frequencies, 
is  by  far  the  most  efficient.  This  method  has  been  applied  only  to  continu- 
ous beams  cn  rigid  supports  and  apparently  is  not  capable  of  extension  to 
continuous  frames  involving  closed  panels.  The  natural  frequencies  of  con- 
tiimwus  beams  on  intermed.1  ate  flexible  supports  may  be  best  determined  by  the 
method  developed  by  Lee  and  Seibel  (4)>  (5)*  This  method  utilizes  principles 

S 

Numbers  in  parentheses,  unless  otherwise  identified,  refer  to  the  Bibli- 
ography at  the  end  of  this  report. 
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that  are  not  vel.1  knovn  to  engineers;  furthermore,  it  presupposes  that  the 
natural  Hodes  of  vibration  of  the  beam  without  the  intermediate  axipporta  are 
known.  This  assumption  restricts  seriously  the  range  of  applicability  of  this 
procedure.  The  determination  of  the  natural  frequencies  of  continuous  frame a, 
involving  closed  panels,  has  apparently  been  attempted  only  by  use  of  the 
classical  method  (6),  which,  as  already  stated,  is  too  laborious  for  practical 
applications. 

The  method  described  herein  is  a generalization  of  Bolder' s method 
(7)  for  calculating  the  natural  frequencies  of  torsional  vibration  of  shafts. 
It  utilizes  well  knovn  engineering  principles  and,  like  Holler’s  method,  it 
is  reduced  to  a routine  scheme  of  computation  which,  when  repeated  a suf- 
ficient number  of  times,  will  give  the  natural  frequencies  of  the  system  to 
any  desired  degree  of  accuracy.  Holier’ s method  has  been  applied  to  the  de- 
termination of  the  natural  frequencies  of  flexural  vibration  of  beams,  fire 
by  Jtyklestad  (8)  and  later  by  Pro  hi  (9),  Rankin  (10),  and  Beilin  (11).  In 
these  studies,  distributed  masses  were  assumed  to  be  lumped  at  a number  of 
stations  along  the  length  of  the  beam  while  the  portion  of  the  beam  between 
these  stations  vas  assumed  to  be  massless.  In  the  method  to  be  presented, 
the  mass  is  assumed  to  be  uniformly  distributed  along  each  member  of  the 
structure.  Abrup.t.  changes  in.  the  magnitude  of  the  distributed  mass  or  of  ‘ • 
the  flexural  rigidity  within  a member  say  be  treated  by  assuming  that  the 
member  is  supported  by  a flexible  support  of  sero  stiffness  at  the  point  of 
the  discontinuity. 

The  principles  underlying  the  method  are  presented  in  Chapter  II 
of  this  report.  This  chapter  also  includes  definitions  of  the  several  physi- 
cal quantities  which  are  necessary  in  the  analysis.  These  quantities  are  tha 
dynamic  stiffnesses  and  the  dynamic  carry-over  factors  which  are  analogous  to 
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those  introduced  by  Hardy  Cross  in  connection  with  the  method  of  moment  distri- 
bution (12).  Extensive  tables  of  numerical  values  of  these  quantities  are  pre- 
sented in  Appendix  a,  while  the  derivation  of  the  governing  equations  is  given 
in  Appendix  Bo  With  these  tabulated  values,  the  calculations  required  in  the 
application  cf  tha  method  to  particular  problems  ere  simplified  immensely. 

The  application  of  the  method  to  continuous  beams  on  rigid  supports 
is  discussed  in  Chapter  III.  In  addition,  several  alternate  methods  of  analysis 
are  considered  and  the  range  of  applicability  and  the  relative  merits  of  each 
are  discussed. 

Chapter  IV  presents  the  extension  of  the  method  to  frames  without 
Sidesway.  For  continuous  beams,  a single  procedure  is  capable  of  determining 
ell  possible  natural  frequencies.  For  frames,  however,  this  procedure  may 
fail  to  detect  those  natural  frequencies  for  which  only  a portion  of  the 
structure  vibrates  while  the  rest  remains  stationary.  A technique  for  over- 
coming this  difficulty  has  been  developed  and  is  presented  also  in  Chapter  IV. 

In  the  study  of  frames,  the  effect  of  permanent  axial  forces  is  neglected. 

The  concluding  section  of  Chapter  IV  is  devoted  to  a discussion  of  the  manner 
in  which  this  effect  may  be  taken  into  account.  Also,  it  is  pointed  out  that 
problems  of  framework  instability  ere  special  esses  of  the  more  general  problem 

of  the  vibration  of  axially  compressed,  pqmbpre*  

• • 

In  Chapter  V,  the  method  is  extended  to  beams  continuous  over  sup- 
ports tnat  are  flexible  instead  of  rigid.  The  resistance  to  deformation  of 
the  supports  is  represented  by  an  equivalent  set  of  mutually  indepenuent  de- 
fies tic-nai  and  rotational  springs.  It  is  shown  that  the  method  may  be  modi- 
fied readily  to  include  the  influence  of  concentrated  rigid  messes,  of  con- 
centrated sprung  masses,  and  of  au  elastic  Subgrade  of  the  Winkler  type. 


Chapter  VI  shows  tho  application  of  the  method  to  sunanetrical,  one- 
bsy,  multi-story  building  frames  which  are  free  to  undergo  side sway® 

Chapter  VII  is  concerned  with  the  extension  of  the  method  to  con- 
tinuous plates  having  two  opposite  edges  simply  supported.  The  pertinent 
expressions  for  dymanic  stiffness  and  dynamic  carry-over  factor  for  plates 
are  presented  in  Appendix  8.  It  is  also  shown  that  numerical  values  of  these 
quantities  may  be  obtained  from  available  tables  of  stiffness  and  carry-over 
factor  for  compressed  plates. 

Table  II  in  Appendix  A gives  influence  coefficients  for  calculating 
the  natural  frequencies  of  vibration  of  systems  composed  of  bars.  It  is 
pointed  out  that  Mttller  -Breslau^  principle  of  influence  lines  is  applicable 
in  the  case  of  steady-state  forced  vibrations,  so  that  these  coefficients  may 
be  interpreted  also  as  coefficients  for  dynamic  fixed-end  moments  produced  by 
a concentrated  harmonic  force* 

Appendix  C includes  a brief  account  of  the  manner  in  vhioh  the  in- 
formation presented  in  this  report  may  be  used  in  the  analysis  of  the  steady- 
state  forced  vibration  of  frames. 

For  convenience  of  reference,  a detailed  outline  of  the  steps  in- 
volved in  the  application  of  the  method  to  each  class  of  problems  is  in- 
cluded in  each  chapter.  In  addition,  several  numerical  examples  are  given 
to  illustrate  the  application  of  the  method  and  to  indicate  convenient  schemed 
for  arranging  the  computations.  Throughout  this  report,  special  effort  has 
been  made  to  present  each  chapter  as  independently  of  the  others  es  possible 
and  to  dincuss  the  numerical  examples  adequately. 

2.  Sign  Convention  and  Rotation. 

The  following  sign  convention  is  used  throughout  this  report  with 
the  exception  of  Appendix  B.  Clockwise  rotations  are  positive.  Internal 
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bonding  moments  acting  at  the  ends  of  a member  (not  a joint)  and  external 
momenta,  except  for  the  restraining  moments  provided  by  rotational  springe, 
are  positive  clockwise.  The  restraining  moment  of  a rotational  spring  is 
positive  when  it  acts  in  s counter-clockwise  direction.  Deflections  are 
positive  downward.  Shears  acting  at  the  ends  of  a member  (not  a joint)  and 
external  forces,  except  for  the  forces  produced  by  deflectional  springs,  are 
positive  downward.  The  restraining  force  of  a deflectional  spring  is  posi- 
tive upward. 

The  letter  symbols  usee  are  defined  where  they  first  appear  in  the 
text  or  by  illustration,  and  they  are  assembled  in  this  section  for  cons 
venience  of  reference. 

Generali 

u)  = circular  frequency  of  vibration 

= natural  circular  frequency  of  vibration 
f = frequency  of  vibration,  in  cycles  per  second 
t = time 

E = modulus  of  elasticity 

For  structures  composed  of  b ars  t 
x = horizontal  coordinate 

I = moment  of  inertia  of  the  cross  section  of  a bar  about  its 
centroidal  axis 


span  length  of  a bar 

mass  per  unit  of  length  of  a bar 

magnitude  of  a concentrated  rigid  mass 


magnitude  of  an  equivalent  concentrated  rigid  mess 
J~E1 — ^ = a ^^®ension^9B8  parameter  for  a bar 

A value  corresponding  to  a netural  frequency 


'bt  fc  wfrifoilffi*''  ht'f  ftilir* ' 4 jfc wt 


e. 


M. 
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‘Ji 


M. 


'J 

?j 

U,  V 

e',5' 

en,5\K,?' 

K,  Q,  T 
k,  q,  t 

CK,CQ,CT 

I*,!*,!1 

K* 

% 

D 

(D) 

B 
C 
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eq 


= deflection  amplitude  at  a point  of  a bar  defined  by  the 
coordinate  i 

= rotation  amplitude  at  support  or  joint  Ji 
= deflection  amplitude  at  support  J^ 

= amplitude  of  internal  bending  moment  at  support  j_  of  a continuous 

beam 

= amplitude  of  internal  bending  moment  st  end  J^  of  a bar  ji_  in  a 
continuous  frame 

= external  moment  at  support  or  joint  J of  a continuous  beam  or 
frsae 

= amplitude  of  shear  at  support  J of  a continuous  beam 
= external  force  at  support  J 
= variable  parameters 

= values  of  9,  £ , M,  and  F due  to  u = 1.0C  and  v = 0 
= values  of  9,  5 , M,  and  F due  to  u = 0 and  v = 1.00 
= dynamic  stiffnesses  for  a bar,  defined  in  Section  5 
= dynamic  carry-over  factors  for  a bar,  defined  in  Section  5 
= dimensionless  coefficients  in  expressions  for  K,  £,  and  _T 

= modified  stiffnesses  for  a bar,  defined  by  Bqs.  (9),  (10)  and  (11) 
= effective  flexural  stiffness  for  a bar,  defined  in  Section  13  * 

= stiffnesses  K,  jP  of  all  bare  meeting  at  support  or  joint 

= effective  flexural  stiffness  of  all  bars  meeting  at  support  or 
joint  _J 

= stiffness  of  a deflectlonal  elastic  spring 

- stiffness  of  an  equivalent  deflectlonal  spring 

- stiffness  of  a rotational  restraint 

= dimensionless  coefficient  in  Eq.  (34)  for  the  deflection  of  a bar 
= axial  force  in  a bar 

= buckling  load  for  a bar  hinged  at  both  ends 
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modulus  of  a continuous  elastic  subgrade 

circular  vibration  frequency  for  a bsr  on  a continuous  elastic 
subgrade 

deflection  amplitude  of  a sprung  mass 

dimensionless  coefficients  defined,  respectively,  by  Eqs.  (24)? 
(73),  (74),  and  (50) 


For  continuous  plates t 
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horizontal  rectangular  coordinates.  The  y-axis  is  taken  parallel 
to  the  pair  of  simply  supported  edges 

span  lengtlis  in  the  x and  y directions,  respectively,  for  a penal 
of  a plate 

density  of  plate  material  in  a particular  panel 

thickness  of  plate  in  a particular  panel 
3 , 

h / 12  ~ moment  of  inertia,  per  unit  width,  for  a particular 
panel  of  the  plate 

Poisson^  ratio 

2l/l-  = flexural  rigidity  of  a particular  panel  of  a plate 

^ = a dimensionless  parameter  for  a panel  of  a pieta 

integer  representing  number  of  half  sine  waves  in  the  distribu- 
tion of  deflections,  slopes,  moments,  etc.,  across  a plate 
width  a 

maximum  rotation  amplitude  along  support  j 

maximum  amplitude  of  bending  moment  at  support 

flexural  stiffness  and  flexural  carry-over  factor  for  a panel 
of  a plate,  defined  in  Section  35 
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II.  METHOD  OF  ANALYSIS 
£•  Basis  of  Method  of  Analysis. 

The  method  used  in  this  report  Is  based  on  the  .fact  that  the  ex- 
citing eouple  or  the  exciting  force  which  is  necessary  to  maintain  a dynamical 
system  in  a steady-state  forced  vibration  with  finite  amplitudes  becomes  equal 
to  zero  at  any  one  of  the  natural  frequencies  of  the  system. 

Figure  1 shows  z members  of  a plane  framework  rigidly  connected  at 
their  common  joint  The  far  ends  of  the  members  may  be  considered  elasti- 
cally restrained  against  both  rotation  and  translation.  These  restraints  are 
furnished  by  the  portion  of  the  structure  not  shown  on  the  figure. 

Consider  that  the  frame  undergoes  a steady-state  forced  vibration 
under  the  action  of  a harmonically  varying  exciting  couple  applied  at  joint  a. 
Joint  z is  assumed  to  be  different  from  joint  o.  The  magnitude  of  the  exciting 
moment  is  assumed  to  be  such  that  the  amplitude  of  either  the  elope  or  of  the 
internal  bending  moment  at  a joint  of  the  structure  different  from  joint  z, 
say  at  joint  1*  has  a prescribed  finite  velue.  The  vibration  of  the  structure 
is  harmonic  and  its  frequency  is  equal  to  the  frequency  of  the  exciting  couple; 
since  the  effect  of  damping  is  neglected,  the  amplitudes  of  vibration  are 
constant  and  the  response  is  either  in  phase  with,  or  180  degrees  out  of 
phase  with,  the  exciting  couple. 

For  a given  system,  the  magnitude  of  the  exciting  moment  necessary 
to  maintain  the  prescribed  amplitude  of  vibration  at  joint  1 depends  on  the 
frequency  of  vibration.  For  the  limiting  case  of  no  vibration,  the  magni- 
tude of  tuw  moment  is  obviously  finite;  its  actual  value  may,  if  desired,  be 
calculated  by  any  of  the  available  methods  of  indeterminate  stress  analysis. 

As  the  frequency  of  vibration  increases  above  zero,  the  structure  is  acted  upon 
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by  Inertia  forces  of  increasingly  greater  magnitudes.  These  forces,  which 
are  distributed  along  the  length  of  the  individual  members,  bring  about  dis- 
tortions in  addition  to  those  produced  by  the  external  dement  acting  static- 
ally. Therefore,  the  amplitude  of  the  dynamic  moment  necessary  to  produce 
the  prescribed  distortion  at  joint  1 may  be  quite  different  from  the  magni- 
tude of  the  corresponding  atatic  moment.  As  the  vibration  frequency  approaches 
any  one  of  the  natural  frequencies  of  the  structure  considered,  the  inertia 
effects  predominate,  and  at  a natural  frequency  the  vibration  is  maintained 
without  any  exciting  moment  acting  permanently  on  the  structure. 

Briefly,  the  method  presented  herein  consists  of  (a)  choosing  a fre- 
quency of  vibration,  (b)  determining  the  magnitude  of  the  exciting  moment 
which,  when  applied  at  joint  z,  will  produce  a vibration  configuration  with  a 
fixed  amplitude  of  slope  or  bending  moment  at  joint  1,  (c)  repeating  the bo 
steps  for  a number  of  assumed  frequencies,  and  (d)  plotting  the  magnitude  of 
the  exciting  moment  es  a function  of  the  frequency  of  vibration.  The  fre- 
quencies for  which  the  exciting  moment  vanishes  represent  natural  frequen- 
cies of  the  system. 

In  the  application  of  this  procedure,  the  following  two  conditions 
must  be  satisfied*  (1)  The  Joint  to  which  a finite  amplitude  of  slope  is 
assigned  should  be  one  that  ia  known  to  rotate  for  all  the  natural  frequen- 
cies that  need  to  be  determined.  If  instead  of  fixing  the  amplitude  of 
slope  the  amplitude  of  moment  is  fixed,  it  must  be  known  that  this  moment 
amplitude  remains  finite.  (2)  The  exciting  couple  must  be  applied  at  a 
Joint  which  i8  known  to  rotate  for  all  the  natural  frequencies  that  need  to 
be  determined.  A couple  applied  at  a joint  which  does  not  rotate  acts 

through  zero  displacement  and  imparts  no  energy  to  the  structure;  therefore, 

•* 

it  does  not  influence  the  natural  frequencies  or  the  vibration  modes  of  the 
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system.  In  such  a case,  the  exciting  moment  may  not  vanish  at  a natural 
frequency.  t 

If  either  ' ? tnese  conditions  is  not  satisfied,  the  procedure 
will  fail  to  reveal  some  of  the  natural  frequencies  of  the  system.  It 

. J* 

will  be  shown  later  that,  for  certain  structures,  it  is  impossible  to 
satisfy  these  requirements.  In  such  cases  * in  order  to  obtain  the  natural 
frequencies  which  the  basic  procedure  fails  to  reveal,  it  becomes  necessary 
to  use  a supplementary  technique  as  described  in  Chapter  17. 

For  an  assumed  frequency  of  vibration,  the  magnitude  of  the  ex- 
citing moment  may  be  determined  by  a number  of  different  procedures.  The 
conditions  to  be  satisfied  are  simply  those  of  equilibrium  and  continuity 
for  each  joint  of  the  structure.  To  satisfy  the  condition  of  equilibrium, 
the  sum  of  the  moments  and  of  the  forces  at  the  ends  of  the  members  meeting 
at  a joint  must  be  respectively  equal  to  zero.  To  satisfy  the  condition 
of  continuity,  the  slopes  of  the  members  meeting  at  a joint  must  be  equal 
and  also  the  deflection  of  the  members  meeting  at  the  joint  must  have  the 
same  magnitude.  These  conditions  may  be  expressed  in  equation  form  in  a 
number  of  different  ways  and  the  equations  may  be  solved  by  a number  of  pro- 
cedures. In  the  method  adopted,  these  conditions  are  expressed  in  the  form 
of  a generalized  alope  deflection  equation,  and  the  distortions  of  the 
structure  at  the  supports  are  computed  by  the  repeated  application  of  this 
equation,  working  progressively  from  one  end  of  the  structure  to  the  other. 

5«  21a otic  Constants  for  a Vibrating  fta-ry 

The  various  quantities  necessary  to  express  the  resistance  to  de- 
formation of  a bar  undergoing  steady-state  forced  vibration  are  defined  in 
this  section. 
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Consider  a bar  fg  with  the  far  end  g fixed.  Let  the  near  end  be 
subjected  to  a harmonically  varying  bending  moment  of  a circular  frequency 
ui  producing  at  that  end  a steady-state  forced  rotation 

6Ct ) “ dc05LUt  . 

The  amplitude  of  the  impressed  moment  may  be  related  to  the  amplitude  of 
the  resulting  rotation  at  end  f by  the  equation 

Mf- Kd  (1) 

The  quantity  K represents  the  moment  required  to  produce  a rotation  of  unit 
amplitude  and  is  defined  herein  as  the  "dynamic  flexural  stiffneseB  of  the 
end  of  the  bar  being  rotated. 

The  moment  induced  at  the  fixed  end  g may  be  written  as 

t1g-«K6.  (2) 

The  quantity  k is  the  ratio  of  the  dynamic  moment  at  the  far  fixed  end  to 
the  moment  at  the  near  end  and  is  defined  as  the  "dynamic  flexural  carry- 
over factor". 

In  the  analysis  of  continuous  beams  on  rigid  supports  and  of  con- 
tinuous frames  for  which  the  joints  do  net  translate,  the  flexural  stiffness 
and  the  product  of  the  flexural  stiffness  and  the  flexural  carry-over  factor 
are  the  only  two  quantities  needed. 

The  foregoing  definitions  are  generalizations  of  .those .originally 
introduced  by  Hardy  Cross  (12)  for  the  analysis  of  frames  subjected  to  static 
loads,  and  they  were  first  used  by  Gaskell  (13),  who  extended  and  applied  the 
method  of  moment  distribution  to  the  problem  of  determining  the  steady-state 
forced  vibration  of  continuous  beams  and  frames  subjected  to  pulsating  loads. 

For  static  conditions,  the  end  reactions  of  the  bar  resulting  from 
the  rotation  of  one  end,  may  be  determined  from  the  end  moments  by  statics. 
For  dynamical  conditions,  this  is  not  possible,  since  tha  bar  is  acted  upon 
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by  inertia  forces  the  distribution  of  which  along  the  length  of  the  bar  lo 
generally  unknown.  The  reactions  must,  therefore,  be  defined  by  additional 
constants . 

The  amplitude  of  the  vertical  reaction  at  the  end  being  rotated  may 
be  written  as 

v,-ae . <3) 

and  the  reaction  at  the  fixod  eixi  may  be  written  as 

Vs  * -yGQ  - -yVf.  (4) 

The  quantity  (J  represents  the  reaction  at  end  f produced  by  a 
rotation  of  unit  amplitude  at  that  end  and  is  defined  as  the  "dynamic 
flexural  shear  stiffness".  The  quantity  £ represents  the  ratio  of  the  re- 
action induced  at  the  far  fixed  end  to  that  produced  at  the  near  end  and 
is  called  the  "dynamic  flexure-shear  carry-over  factor". 

Consider  now  that  end  f is  subjected  to  a harmonically  varying 
force  producing  at  that  end  a steady-state  forced  deflection  without 
rotation,  such  that  the  magnitude  of  the  deflection  is 

5(t ) ” Scoswt . 

The  amplitudes  of  the  force  and  of  the  deflection  at  tne  left  end  may  be 
related  by  the  expression 

Vf  m T5  (5) 

The  quantity  T denotes  the  force  necessary  to  cause  a deflection  of  unit 
amplitude  and  is  defined  as  the  "dynamic  shear  stiffness"  for  the  end  being 
deflected.  The  reaction  at  the  far  fixed  end  may  be  written  as 

Vj  ~iTo  - -i  Vf  (6) 

The  quantity  t shows  the  ratio  of  ihe  reaction  at  the  far  end  to  that  at  the 
near  end  and  is  called  the  "dynamic  shear  carry-over  factor".  The  amplitude 
of  the  ntoEGnt  induced  at  the  end  being  deflected  is 
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nf  = Q5.  (?) 

and  the  moment  at  the  far  fixed  end  1b 

n2  * yQ8  • ytlf  . (8) 

The  quantities  (j  and  a are  the  same  as  those  used  in  Eqa.  (3)  and  (4)*  That 
these  should  be  the  same  follows  from  a reciprocal  theorem  given  by  Lord 
Rayleigh  (14),  which  is  the  dynamic  equivalent  of  Maxwell's  Law  of  reciprocal 
relations. 

Throughout  this  presentation,  the  members  composing  the  structure 
are  considered  to  be  uniform.  The  stiffness  and  the  carry-over  factors  for 
both  ends  of  such  members  are  equal. 

The  notation  used  for  the  various  stiffnesses  and  carry-over  factors 
is  the  same  as  that  used  by  Newmsrk  (15)  in  his  static  analysis  of  slabs  con- 
tinuous over  flexible  supports.  The  notation  is  summarized  in  Fig.  2.  The 
derivation  of  the  algebraic  expressions  for  the  various  stiffnesses  and 
carry-over  factors  is  given  in  Appendix  B* 

For  certain  conditions  of  symmetry,  antisymmetry,  and  for  those 
cases  for  which  the  degree  of  restraint  at  the  far  end  of  a member  is  known, 
it  is  convenient  to  use  effective  stiffnesses.  The  pertinent  expressions  for 
these  stiffnesses  are  the  same  as  those  for  the  static  case.  Expressions  are 
given  here  only  for  effective  flexural  stiffnesses.  The  particular  cases  con- 
sidered end  the  symbols  used  to  identify  them  are:  g' , when  the  far  end  of 

the  bar  is  prevented  from  deflecting  and  is  elastically  restrained  against 

H C 

rotation;  Z , when  the  far  end  is  simply  supported;  when  the  bar  is  on 
rigid  supports  ard  its  deformation  is  symmetrical,  and  g\  when  the  bar  ia 
on  rigid  supports  end  its  deformation  is  antisymmetrical.  It  can  be  proved 
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l"  = X(l-k2)  , 


^ = K(l-k)  , 
X1  = K(l+k)  . 


The  stiffness  of  a given  bar  1b  also  equal  to  the  stiffness  K of  a similar 

i 

bar  one  half  as  long*  The  expression  for  X is  given  in  article  12,  where  it 
is  used  first* 

It  is  possible  to  derive  also  expressions  for  effective  shear  stiff- 
nesses* However,  these  are  not,  in  general,  as  simple  and  convenient  to  use 
as  the  effective  flexural  stiffnesses . 


(9) 


(10) 


(11) 


6.  Numerical  Values  of  Stiffness  and  Carry-Over  Factors* 

All  carry-over  factors  are  dimensionless  and  depend  on  a single  di- 
mensionless parameter 

imu)1  , , . 

* ~ J El  L’  ( 12 ) 

in  which  m = the  mass  per  unit  of  length  of  the  bar, 

u)  = the  circular  frequency  of  vibration,  as  previously  noted, 

B = the  modulus  of  elasticity  of  the  material  in  the  bar, 

I = the  moment  of  inertia  of  the  cross  section  of  the  bar  about  its 
centroidal  axis,  and 

L = the  Bpan  length  of  the  bar* 

The  various  stiffnesses  are  determined  from  the  expressions 

* Q - Ca- 71“'  T - Cr^rr  * 03) 

where  the  C‘s  are  dimensionless  coefficients  depending  on  the  parameter  A • 

A graphical  representation  of  the  variation  with  A of  the  various 
carry-over  factors,  stiffness  coefficients,  and  of  their  products  is  given  in 
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Figs.  3 through  12.  It  is  noted  that  the  curves  in  these  figures  range  be- 
tween minus  infinity  and  plus  infinity.  The  A values  corresponding  to  the 
zero  ordinates  and  to  the  discontinuities  of  the  curves,  represent  natural 
frequencies  of  bars  having  standard  boundary  c onditions.  For  example,  con- 
sider the  curves  in  Figs.  3 and  5 for  the  flexural  stiffness  and  the  flexural 
carry-over  factor.  Values  of  \ equal  to  3*927,  7.069  and  10.210  correspond, 
respectively,  to  the  first,  the  second,  and  the  third  natural  frequencies  of 
a hinged-fixed  bar.  At  these  frequencies,  no  exciting  moment  is  required  to 
maintain  the  vibration;  consequently,  the  value  of  dynamic  stiffness  is  equal 
to  zero-.  Furthermore,  3ince  the  moment  at  the  fixed  end  of  the  bar  has  a 
finite  magnitude,  the  carry-wer  factor  for  the  member  becomes  infinite  at 
these  frequencies.  Values  of  A equal  to  4*730  and  7.S53  correspond,  re- 
spectively, to  the  first  and  the  second  natural  frequencies  of  a bar  fixed 
at  both  ends.  At  these  frequencies,  the  end  moments  have  a finite  value 
while  the  rotations  of  the  ends  are  zero;  accordingly,  the  stiffness  of  the 
member  has  an  infinite  value.  For  the  case  of  no  vibration,  A = 0*  the 
various  quantities  in  Fig.  3 through  12  assume  the  well  known  static  values 
of 
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Numerical  values  of  the  carry-over  factors,  of  the  coefficients  of 
the  various  stiffnesses,  and  of  their  products  are  given  in  Table  I of  Ap- 
pendix A.  All  values  are  reported  to  seven  significant  figures  for  the  range 
of  A from  zero  to  10.20  at  increments  of  0.01.  In  some  cases,  the  accuracy 
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of  the  seventh  significant  figure  reported  is  uncertain.  These  values  were 
computed  by  use  of  the  Electronic  Digital  Computer  of  the  University  of 
Illinois . 

7.  Numerical  Values  of  Deflections  dua  to  Bid  Rotations t 

Consider  the  beam  shown  in  Fig.  2a  with  the  left  end  subjected  to 
a steady-state  forced  rotation  $•  The  deflection  amplitude  of  the  beam  at 
a distance  3 1 from  the  left  end  may  be  written  as 

Ix  = C6L  t (14) 

where  C is  a dimensionless  coefficient  dependent  on  the  value  of  2 and  the 
parameter  A «•  If  © represents  the  rotation  at  the  right  end  instead  of  at 
the  left  end  of  the  beam,  the  deflection  amplitude  at  a distance  x from  the 
right  end  will  be  equal  to  the  right  hand  side  of  Bq.  (14)  multiplied  by- 
minus  one.  The  minus  sign  is  a consequence  of  the  sign  convention  adopted* 
Numerical  values  of  C are  given  in  Table  II  of  Appendix  A for  suc- 
cessive twelveth  points  of  the  bean  for  values  of  A ranging  between  aero  and 
10.20,  These  Values  were  computed  from  Eq.  (B-32)  in  Appendix  B,  by  use  of 
the  Electronic  Digital  Computer  of  the  University  of  Illinois.  The  values 
are  reported  to  five  significant  figures,  but  to  no  more  than  six  decimal 
places. 

The  values  in  Table  II  may  also  be  interpreted  as  coefficients  of 
dynamic  fixed-end  moment  for  s beam  subjected  to  a puleating  concentrated 
force.  This  follows  from  Muller-Broslau } s principle,  which  it  can  be  shown 
to  hold  true  for  dynamical  systems  undergoing  steady-state  forced  vibrat-ion* 
This  principle,  as  applied  to  the  dynamical  case,  is  presented  in  Appendix  B* 


: A 


M i 

$ 

& 


i 


19 


III.  APPLICATION  OF  METHOD  TO  CONTINUOUS 

» 

BEAMS  ON  RIGID  SUPPORTS 


8.  General. 

The  beams  considered  are  assumed  to  be  straight  and  may  have  any 
number  of  spans  of  arbitrary  length.  At  their  extreme  ends  they  may  be 
hinged,  fixed,  or  only  partially  fixed  by  means  of  rotational  restraints 
which  are  assumed  to  be  proportional  to  the  end  rotations.  The  cross  section 
and  the  mass  per  unit  of  length  of  the  beam  may  vary  from  one  Bpan  to  the 
other,  but  in  any  one  span  these  quantities  are  considered  constant.  It  is 
assumed  that  vibration  is  restricted  to  one  of  the  principal  planes  of  flexure 
of  the  beam,  and  that  the  cross  sectional  dimensions  of  each  span  are  small 
in  comparison  to  its  length  so  that  the  effects  of  shearing  deformation  and 
rotatory  inertia  are  negligible. 

The  supports  of  the  beam  are  numbered  successively  from  left  to 
right  starting  with  1 at  the  extreme  left  end  and  terminating  with  z at  the 
extreme  right  end. 

The  portion  of  the  beam  between  two  consecutive  supports  and  J+l 
is  referred  to  as  th9  ^-th  sp“Ur  The  quantities  Lj,  E^,  1^,  K^,  and  kj 

refer  to  the  J-th  span. 

©.  denotes  the  amplitude  of  rotation  of  the  deflected  beam  over  the 
*1 

jh-th  support  and  Mj  denotes  the  amplitude  of  bending  moment  across  a section 

at  the  same  support.  The  subscripts  L and  R designate,  respectively,  sections 

just  to  the  left  and  just  to  the  right  of  the  support.  M,  denotes  the  ampli- 

•> 

tude  of  the  external  couple  at  support  J[. 


9«  Develoonent  of  the  Basic  Equations. 


Figure  33  shows  the  extreme  deflected  position  of  spane  i-1  end  £ of 
a continuous  beam  undergoing  & steady-state  forced  vibration.  The  vibration 
is  assumed  to  be  maintained  by  an  exciting  couple  applied  at  the  extreme  right 
end  of  the  beau.  There  is  no  other  exciting  force  or  moment  acting  on  the 
system. 

In  Fig.  13 t the  rotations  and  bending  moments  at  the  ends  of  each 
span  are  indicated  in  their  positive  directions.  The  slope  and  the  bending 
moment  at  a time  t for  eupport  £ are 

Oj(t)  - Oj  cos  u)t , ond  Mj(t)  ~ Hj  cos  wt  . (15) 


In  thG  equations  to  be  used  the  cosurt  appears  as  a common  factor;  for  con- 
venience, this  will  be  omitted,  and  in  the  remainder  of  this  discussion  the 
terms  "amplitude  of  slope"  and  "slope"  and  the  terms  "amplitude  of  moment*  and 


"moment"  will  be  used  interchangeably. 

To  insure  continuity  and  equilibrium  of  the  beam  over  the  interior 
support  2*  it  is  required  that 


(Oik 


m*  - e; 


(16) 


tlj  * (71;  k + (Mj)ft  - o - 


(17) 


The  moments  (M,)t  and  (M,)b  can  now  be  expressed  as  functions  of  the 
j u ] » 

end  rotations  of  the  two  spans  &s  follows*  Consider  span  First,  assume 

that  the  right  end  of  the  span  is  held  fixed  while  the  left  end  is  rotated 

through  an  angle  then,  the  moment  at  the  end  being  rotated  is  equal  to 

the  product  of  the  rotation  S,  and  the  flexure!  stiffness  cf  the  member  £,• 

J J 

Next,  imagine  that  the  left  end  of  the  span  is  kept  fixed  while  the  right  end 
is  rotated  through  the  moment  induced  at  the  fixed  left  end  is  equal  to 

the  product  of  the  rotation  and  the  product  of  the  flexural  stiffness 
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and  flexural  carry-over  factor  of  the  member  (kX)j.  Since  the  principle  of 

superposition  holds  true,  the  moment  (M,)n  corresponding  to  the  rotations  ©. 

j H J 

and  ®j+2  *s  the  sum  of  these  partial  momenta. 

(Hj)r  - Kjdj  + rttfOj  eJfl  . (X 


Considering  span  j-1,  one  obtains  in  a similar  manner: 

(rij)L  - Kj-,dj  + (*K)j-,6j-,  . 


(18a) 


(13b) 


Substituting  Bos.  (18a)  and  (18b)  in  Eq.  (17)  and  solving  for  oae  obtains 

the  following  equation  relating  the  slopes  over  three  consecutive  supports  of 


a continuous  beam: 


djti 


(Kj- , * Kj)dj  + (*K)j-,  dj-, 
(kK)j 


(19a) 


This  equation  is  a generalized  elope  deflection  equation  with  the  deflection 
term  missing*  It  will  be  referred  to  as  the  "three  slope  equation".  Eq.  (19a) 
Is  applicable  only  to  interior  supports;  the  appropriate  relations  for  the  end 
supports  are  given  in  the  following  paragraphs. 

It  is  assumed  that  the  extreme  ends  of  the  beam  are  elastically  re- 
strained agsinst  rotation.  The  relationship  between  the  end  moments  and  end 


rotations  are 


n,  « -R,e,  , 


nz  - -Rzdz , (2i 

where,  R^  and  are  the  known  stiffnesses  of  the  rotational  restraints  at 
the  left  and  the  right  ends,  respectively.  For  a hinged  end,  R = 0,  and  for 
a damned  end,  R = infinity.  The  negative  signs  in  these  expressions  follow 
fro*"  the  sign  convention  used  and  indicate  that  for  a positive  restraint. 


.J 
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the  moment  exerted  on  the  beam  by  the  restraint  ..ots  in  a direction  opposite 
to  the  direction  of  rotation  of  the  beam. 

The  moments  and  Ma  can  also  be  expressed  by  the  following  equa- 
tions, obtained  respectively  from  2qs.  (18a)  end  (18b). 

n,  - K,0,  + (*K),Qt  . (18a’) 

nt  - Kz.se:  +(*!<)'-, ez  . dsb’) 

Eliminating  between  Bqs.  (18a  ) and  (20)  and  between  Eqs.  (13b*)  and  (21), 


one  obtains 


(R,  + K,)d,  + (4tK),  e7  « 0 , 


(21a) 


(Kz.+RI)eI*<mt.,et.rO . (22.) 

At  a natural  frequency,  both  of  these  equations  must  be  satisfied  identically. 

Equations  (21a)  and  (22a)  apply  only  to  hinged  and  to  partially 
fixed  ends.  For  fixed  ends,  the  equations  are  specialized  as  follows*  For 
= 0,  the  relation  between  the  moment  at  the  fixed  end  and  the  rotation  of 
the  beam  over  thesBcond  support  is  obtained  from  Bq.  (18a  ) as 

n,  - (*JO,e,  . (21b) 


If  the  right  end  is  fixed, 


dz  - 0 , 


(22b) 


and  the  criterion  for  a natural  frequency  is  that  Eq.  (22b)  be  satisfied. 

The  magnitude  of  the  moment  at  the  fixed  end  is  of  no  interest  but,  should 
it  be  desired,  it  may  be  calculated  from  Eq.  (18b  ),  keeping  in  mind  that  9^=0. 

10.  Outline  of  the  Procedure. 

The  procedure  for  arriving  at  the  natural  frequencies  of  a continuous 


4 

"1 


f 

% 

% 


i 


beam  may  be  outlined  as  follows* 
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1.  A fixed  value  is  assigned  to  the  amplitude  of  elope  or  bending 
moment  at  the  first  support  of  the  beam.  Since  the  natural  fre- 
quencies of  a system  depend  only  on  the  relative  values  of  the  de- 
flection* any  arbitrary  amplitude  consistent  with  the  actual  bound- 
ary conditions  may  be  chosen.  i'or  a hinged  or  for  a partially  fix- 
ed end,  ©^  is  taken,  for  convenience,  equal  to  unityj  for  a clamped 
end,  ©^  is  equal  to  aero,  and  or  times  the  i/EI  of  some  re- 
ference span  is  taken  equal  to  unity  instead* 

2*  A trial  frequency  of  vibration,  in  , is  chosen  and  the  A values  for 
all  spans  are  evaluated.  These  calculations  are  carried  out  con- 
veniently in  a tabular  form,  as  illustrated  in  Example  2. 

3.  With  the  A values  available,  the  flexural  stiffness  and  the  product 
of  the  flexural  stiffness  and  flexural  carry-over  factor  for  each 
span  of  the  beam  are  found  from  Table  I in  Appendix  A* 

4*  The  rotation  of  the  beam  over  the  second  support  is  determined  from 
Bq.  (21a)  or  (21b). 

5*  By  successive  applications  of  Bq*  (19)*  the  rotations  ®aar® 

evaluated*  A convenient  tabular  scheme  for  arranging  the  computations 
is  described  in  Example  2* 

6*  If  support  z is  fixed,  the  determination  of  the  rotation  © completes 

z 

• > ene  cyel-a-  sf  the  procedure’  (bee’  jsq.  22b).  However,  if  this  support 
is  hinged  or  is  only  partially  fixed,  it  is  necessary  to  carry  out 
the  additional  step  of  evaluating  the  left  hand  aide  of  Eq*  (22a). 

7*  Steps  1 through  6 are  repeated  for  different  assumed  frequencies,  and 
the  values  calculated  for  the  left  hand  side  of  Bq*  (22a)  or  (22b) 
are  plotted  as  a function  of  the  assumed  frequencies,  or  what  is 
usually  more  convenient,  as  a function  of  the  corresponding  A 
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values  for  some  reference  span.  The  zero  intercepts  of  the  re- 
sulting, curve  which  is,  in  general,  similar  in  shape  to  that  shown 
in  Figs.  16  ana  IS,  correspond  to  the  natural  frequencies  of  th» 
system. 

11.  Determination  of  Modes  of  Vibration 

Since  the  rotations  of  the  beam  over  the  cupports  are  evaluated  in 
each  cycle  of  this  procedure,  the  deflection  configuration  of  the  beam  for 
any  desired  frequency  can  ordinarily  bo  sketched  from  these  rotations.  The 
natural  inodes  of  free  vibration  may  be  determined  .^rom  the  rotations  corre- 
sponding to  the  natural  frequencies  in  the  same  manner. 

If  it  is  desired  to  compute  these  deflections  accurately,  it  is 
necessary  to  use  the  numerical  coefficients  given  in  Table  II.  The  de- 
flection at  any  point  within  a span  may  be  obtained  by  adding  (a)  the  de- 
flection produced  by  the  rotation  of  the  left  end  of  the  spen,  assuming 
that  the  right  end  is  fixed  and  (b)  the  deflection  produced  by  the  rotation 
of  the  right  end  of  the  span,  assuming  that  the  left  end  is  fixed. 

12.  Illustrative  Examples. 

Example  1,  Consider  a uniform  beam  continuous  over  five  rigid 
supports  spaced  equidistantly.  The  beam  is  simply  supported  at  the  left 
end  sad  fixed  at  the  right  end,  as  shown  in  Fig.  14.  It  is  desired  to  cal- 
culate  its  first  eight  natural  frequencies  and  the  corresponding  natural  modes 
of  vibration.  It  is  assumed  that  the  beam  is  cut  at  •the  extreme  right  end  and 
then  an  exciting  moment  is  applied  there.  It  a natural  frequency,  the  magni- 
tude of  this  moment  must  be  ouch  that  the  condition  5^  = 0 is  satisfied. 

The  amplitude  of  slope  at  the  extreme  left  end  is  taken  equal  to  unity.  Since 
all  spans  are  identical,  rather  than  repeating  the  procedure  outlined  in  Sec- 
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tion  10  for  each  assumed  frequency  of  vibration,  it  is  more  convenient  to 
derive  a general  expression  for  6*.  and  determine  directly  from  this  expression 
the  natural  frequencies  of  the  beam. 

Let  K be  the  flexural  stiffness  and  k the  flexural  carry-over  factor 
for  each  span.  These  quantities  depend,  of  course,  on  the  parameter  A * 

From  Eq.  (21a)  one  obtains 


a. 


K 

kK 


I 

* 


Applying  Eq.  (19)  successively  to  joints  2,  3»  and  4»  one  obtains 

z-k* 


e , - t *j  + *k 

d4--K[z-*j£>*{-ip  WC- 

6S  ‘ -K  [z  ^r-  * * 2'*' 


-4 


*kK- 


£~8k*+8 


The  expression  for  was  evaluated  for  several  values  of  A and  the  results 
were  used  to  plot  the  curve  shown  in  Fig.  16.  The  values  of  k corresponding 
uc  the  assumed  values  of  A were  obtained  from  Table  X.  The  A values  cor- 
responding to  the  natural  frequencies  are 


Order  of 

Value  of 

Order  of 

* 

Value  of 

An 

An 

An 

An 

1 

3.21 

5 

6.36 

2 

3.65 

6 

6.79 

3 

4.21 

7 

7.34 

4 

4.655 

8 

7.78 

These  values  agree  with  those  reported  elsewhere  (16).  The  circular  natural 


frequencies  are  obtained  from  the  expression 
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and  the  natural  frequencies,  in  cycles  per  second,  are  computed  from 

Z w 

In  order  to  determine  the  natural  modes  of  vibration,  first,  the  k 
values  corresponding  to  the  natural  frequencies  were  determined,  and  then  the 
expressions  of  ©,,,  and  were  evaluated.  The  results  are  summarized  in 
the  following* 


i 

Order  of 
Mods 

ei 

e2 

«3 

°5 

1 

1.00 

— 

-.924 

.707 

—383 

0 

2 

1.00 

-.383 

-.707 

.924 

0 

3 

1.00 

-.707 

—924 

0 

4 

1.00 

. -1 

.707 

.383 

0 

5 

1.00 

.924 

.707 

.383 

0 

6 

1.00 

- .383 

-.707 

-.924 

0 

7 

1.00 

CT\ 

to 

o 

* 

1 

-.707 

.924 

0 



1.00 

-.924 

.707 

-.383 

0 

From  these  rotations,  the  shapes  of  the  natural  modes  of  vibration  can  be 
sketched.  For  this  particular  problem,  the  vibration  modes  were  computed 
by  use  of  the  numerical  values  given  in  Table  II,  following  the  procedure 
described  in  the  preceding  Article.  For  the  purpose  of  illustration,  the 
computations  involved  in  the  determination  of  the  first  natural  mode 
( A = 3*21)  are  presented  in  detail. 

Deflection  of  span  1 at  successive  1/6  - points: 


for  0/  = 1.00,  02  = 0; 

0 

.128 

.179 

.163 

.103 

.033 

0 

for  G,  = 0,  ©j  = -.924: 

0 

.031 

.095 

.151 

.166 

.118 

0 

total: 

0 

.159 

.274 

.314 

.269 

.151 

0 

Deflection  of  span  2 at  successive  1/6  - points: 
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for  0,  = 

—.924,  ©3  — 0: 

0 

-.118 

-»166 

-.151 

-.095 

-.031 

0 

for  82  = 

0,  %3  = .707: 

0 

-.024 

-.072 

-.115 

-.127 

-.090 

0 

total: 

0 

-.342 

-.238 

-.266 

-.121 

c 

Deflection  of 

span  3 at  successive  1/6  - points: 

for  ©j  = 

.707,  ©4  = 0: 

0 

.090 

.127 

.115 

.072 

.024 

0 

for  = 

0,  ©4  = -.383: 

0 

.013 

.039 

.062 

.069 

.049 

0 

total: 

0 

.103 

.166 

.177 

.141 

.073 

0 

Deflection  of 

span  4 at  successive  1/6  - points: 

for  ©4  = 

-.383,  = 0: 

0 

-.049 

-.069 

-.062 

-.039 

-.013 

< 

The  first  six  natural  modes  of  vibration  are  shown  in  Fig.  17 
Example  2.  In  order  to  illustrate  several  additional  details  of  the 
procedure  and  present  a convenient  tabular  scheme  for  recording  the  compu- 
tations for  the  general  case  in  which  the  dimensions  of  the  beam  may  vary 
from  span  to  span,  we  consider  the  four-span  continuous  beam  shown  in  Fig.  15* 
Only  the  first  five  natural  frequencies  will  be  evaluated.  The  beam  is  as- 
sumed to  be  elastically  restrained  at  the  left  end  and  hinged  at  the  right 
end.  The  stiffness  of  the  end  restraint  and  the  characteristics  of  the 
various  spans  are  shown  in  Fig.  15  . 

For  convenience  in  carrying  out  the  calculations,  the  natural  fre- 
quencied  of  the  system  are  expressed  in  terms  of  the  pertinent  properties  of 
some  reference  span,  say  span  r.  In  this  particular  example  we  take  r = 1. 

In  terms  of  the  A value  of  the  r-th  span,  the  A value  for  any  span  J is 

(23) 


El 


In  terms  of  the  — of  the  r-th  span,  the  stiffness  and  the  product  of  ths 
stiffness  and  of  the  carry-over  factor  for  any  span  j are  equal  to  the  values 
obtained  from  Table  I multiplied  by  the  dimensionless  factor 
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Equations  (23)  and  (24)  can  be  verified  readily* 

The  quantities  Aj/Ar  and  exj  are  evaluated  in  Table  1A.  It  should 
be  noted  that  the  calculations  in  this  table  are  independent  of  the  frequency 
of  vibration. 

The  ti*ial-and -error  produce  tor  determining  the  natural  frequencies 

of  the  system  is  carried  cut  in  Table  IB*  As  an  example  of  the  use  of  this 

table,  a complete  cycle  of  calculations  is  carried  out  for  a trial  value  of 

A r*  A,  - 2*40.  This  value,  shown  encircled  in  the  r-th  line  of  Column  (2), 

corresponds  to  a circular  frequency  of  vibration  w . The 

L,  J rn, 

arrangement  of  the  various  quantities  in  this  table  is  believed  to  facilitate 
the  computational  work  and  to  reduce  substantially  the  probability  for  errors. 
The  order  in  which  the  columns  in  this  table  are  filled  in  is  indicated  by 
the  following  sequence  of  column  numbers:  (l),  (3),  (2),  (4  and  8),  (5)» 

(7),  and  (6).  Columns  (l)  and  (3)  are  reproduced,  respectively,  from  Columns 
(5)  and  (6)  of  Table  A.  The  A values  for  the  various  spans  in  Column  (2) 
are  obtained  as  the  product  of  the  assumed  Ar  and  each  of  the  values  in 
Column  (l).  Columns  (4)  and  (8)  give,  respectively,  values  of  the  stiffness 
and  of  the  product  of  the  stiffness  and  the  carry-rover  factor  for  each  span, 


E.I. 


in  terms  of 


'i 


• these  quantities  are  obtained  directly  from  Table  I in 


Appendix  A,  using  the  A values  computed  in  Column  (2).  Column  (5)  gives  the 

V L 

total  stiffness  of  the  spans  adjoining  each  support  in  terms  of  . 

li 

r 

The  value  for  the  ,1-th  line  in  this  column  is  determined  by  taxing  the  srnm 
of  the  products  of  ihe  values  in  Columns  (3)  and  (4)  for  lines  j-1  and 
Column  (7)  gives  the  product  of  the  stiffness  and  of  the  carry-over  factor 
for  each  span  in  terms  of  — j the  entries  in  this  column  are  obtained 
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by  multiplying  the  entries  in  Column  (S)  by  those  in  Column  (3)»  Column  (6) 
gives  the  rotation  of  the  beam  over  the  supports.  The  first  value  in  this 
column  is  unity.  (Had  the  beam  been  fixed  at  the  left  end.  this  value  would 
have  been  zero).  Tho  second  value  in  the  column,  ©2»  is  evaluated  from  Eq. 
(21a) 

, . (o.?ooo  0.6649)  1,0000  = _1>8466 

2 2.2555 


This  operation  is  not  indicated  in  the  Table.  (Had  support  1 been  fixed, 
Eq.  (21b)  would  have  been  used  instead).  The  values  of  to  ©z  are  de- 
termined from  the  values  in  Columns  (5)  and  (7)  by  use  of  Eq.  (19a),  which, 
in  terms  of  column  members,  takes  the  form: 


= _ (5)j(6)j  + (Vj-,(7)j-, 

C?#Af  / 1 \ . 


</  + # 


(for  j £ 2). 


(19b) 


Thus, 


e s . 6Q061(-1.8466)  ♦ 1.0000(2.2^  s , 6l85 

2.0330  4.0X0? 


The  left  hand  side  of  Eq.  (22a),  evaluated  at  the  bottom  of  the  table,  is 

®1^1 

found  to  be  equal  to  17.55  . - 

n 

Since,  for  the  assumed  value  of  A,  = 2.40,  Eq.  (22a)  was  not  satis- 
fied, this  value  does  not  correspond  to  a natural  frequency  of  the  system. 


The  physical  significance  of  the  computed  value  of  17.55  ~ 


hh 

h. 

follows:  the  negative  of  this  value  divided  by  the  rotation  ©2* 


is  as 


_ n.fi  = 
21.463  L 


-0.8179 


represents  the  stiffness  of  a rotational  constraint  which,  if  it  were  im- 
posed at  the  right  end  of  the  beam,  would  have  made  the  assumed  frequency 
correspond  to  a natural  frequency  of  the  system. 

8y  repeating  several  such  cycles  of  computation  for  different  values 
of  X,  , the  curve  in  Fig.  18  was  obtained.  The  first  fivt.  critical  values 
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ere  recorder1  on  +he  figure.  The  corresponding  circular  natural  frequencies 


are 
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If  it  is  desired  to  evaluate  these  quantities  more  precisely,  the  compu- 
tations should  be  repeated  for  several  additional  values  of  A|  in  the 
neighborhood  of  the  critical  values,  and  the  results  Bhould  be  plotted  on 
a larger  scale. 

The  natural  inodes  of  vi bra tions determined  in  the  manner  described 
in  Section  11,  are  shown  in  Fig.  19.  It  should  be  stated  that,  in  general, 
for  the  fundamental  or  lowest  natural  frequency,  the  rotations  of  the  beam 
over  the  supports  are  not  very  sensitive  to  the  magnitude  of  the  frequency 
of  vibration.  For  some  of  the  higher  vibration  frequencies,  however,  a 
slight  variation  in  the  value  of  the  frequency  may  affect  the  rotations 
materially.  Accordingly,  the  accurate  evaluation  of  the  rotations  in  these 
latter  cases  may  become  somewhat  cumbersome. 


lit  Alternate  Methods  of  Analysis. 

As  applied  to  continuous  beams,  the  criterion  for  a natural- fre- 
quency is  that 

hz  - o . 

It  is  presumed  that  support  z ie  not  fixed. 

As  has  already  been  remarked,  the  method  used  to  evaluate  M is 

z 


(22a) 
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merely  one  of  a number  of  possible  methods*  It  is  the  purpose  of  the  follow- 
ing discussion  to  present  several  alternate  procedures  for  arriving  at  the 
same  result. 

The  Effective  Stiffness  Criterion,  The  moment  at  a Joint  of  a structure 

necessary  to  produce  at  that  joint  a rotation  of  unit  amplitude,  while  all 

other  joints  ere  in  their  actual  condition  of  restraint,  is  defined  as  the 

"effective  flexural  stiffness"  of  the  joint.  This  quantity  depends  on  the 

properties  of  all  the  members  of  the  structure,  and  it  will  be  denoted  by  K . 

Let  K*  represent  the  total  effective  stiffness  at  the  right  head 
z 

support  z of  a continuous  beam;  then,  Eq.  (22a)  may  be  written  as 

K i 4 - 0 . (25) 

Since  G^  is  assumed  to  be  different  from  zero,  this  equation  is  satisfied 
only  if 

K'  - O.  (26a) 

It  should  be  emphasized  that  Eqs.  (22a)  and  (26a)  express  identically  the 
sam~  condition,  only  in  slightly  different  forms. 

Equation  (26a)  represents  the  effective  stiffness  criterion  for 
determining  natural  frequencies.  This  criterion  will  now  be  applied  by  use 
of  the  moment  distribution  procedure  and  a procedure  which,  for  want  of 
any  better  term,  will  be  referred  to  as  the  "direct"  procedure. 

The  Moment  Distribution  Procedure.  Gaskell's  adaptation  of  the  method 
of  moment  distribution  (13)  may  be  applied  as  follows: 

f is  ananwArl  , .q ry^  fl  S"ti?f*2lSSS  j K ; 

and  the  flexural  carry-over  factor,  k,  for  each  span  of  the  beam 
are  computed  from  Table  I in  Appendix  A. 
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2.  With  all  joints  of  the  structure,  except  joint  z>  fixed  against 
rotation,  an  exciting  moment  is  applied  at  joint  z producing  a 
rotation  of  unit  amplitude  at  that  joint.  Obviously,  the  fre- 
quency of  this  moment  is  equal  to  the  assumed  frequency  of  vibra- 
tion and  the  magnitude  of  the  moment  is  equal  to  Kz.t  t Rz  . 

3.  This  moment  is  distributed  to  the  adjacent  members  in  proportion 
to  their  relative  stiffness,  and  the  proper  proportion  of  the 
balancing  moment  is  carried-over  to  joint  z-1 . 

4.  Joint  z is  then  locked,  and  the  unbalanced  moment  at  joint  z-1  is 
distributed,  carried  over,  and  balanced  through  the  rest  of  the 
structure.  During  this  process  of  moment  balancing,  joint  z is 
maintained  locked. 

5.  The  total  moment  carried  back  to  joint  z is  determined.  Finally, 
the  effective  stiffness  of  the  joint  is  computed  as  the  algebraic 
sum  of  the  moment  applied  initially  to  the  joint  and  the  moment- 
carried  back  after  all  the  other  joints  have  been  balanced. 

6.  Steps  1 through  5 are  repeated  for  several  frequencies  of  vibration, 

and  the  natural  frequencies  are  determined  as  those  frequencies  for 

« 

which  the  effective  stiffness  vanishes. 

The  "Direct11  Procedure.  The  second  procedure  for  applying  the  effective 
stiffness  criterion  is  presented  in  this  section.  Consider  a bar  on  unyielding 
supports  at  each  end  with  one  end  elastically  restrained  against  rotation. 

The  restraint  may  be  due  to  an  actual  coil  spring  or  it  may  symbolize  the  con- 
tinuity of  the  bar  with  adjoining  members.  The  stiffness  of  the  restraint  is 

*At  this  point,  attention  should  be  called  to  the  fact  that  the  method  of 
moment  distribution  does  net  converge  always  to  an  answer#  Therefore,  this 
method,  which  probably  would  appeal  to  many  engineers,  ia  restricted  in  its 
practical  application.  Thi8  fact  i£  considered  in  somewhat  greater  detail  in 
Section  14 • 
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denoted  by  R.  It  can  be  shown  (17)  that  the  effective  flexural  stiffness  of 
the  opposite  end  of  the  bar  is  given  by  the  expression 

(*Kf 


K'  - K~ 


(27) 


K+R 

This  equation  may  be  used  to  calculate  the  effective  Stiffness  of  a 
continuous  beam  as  follows: 

!•  A frequency  of  vibration  is  assumed,  and  the  corresponding  values 
of  K and  kK  for  each  span  of  the  beam  are  determined  from  Table  I 
in  Appendix  A. 

2*  With  the  stiffness  of  the  restraint  R^  at  the  extreme  left  end  of 
the  beam  known,  the  value  of  effective  stiffness  of  the  first  span, 

, is  computed  from  Bq«  (27).  This  value  represents  also  the 
stiffness  of  the  rotational  restraint  exerted  by  the  first  span 
on  the  left  end  of  the  second  span.  By  application  of  Eq.  (27)  to 
consecutive  spans,  the  effective  stiffness  of  spans  2 to  z-1  are 
evaluated. 

t 

3.  Having  determined  Kz-1’  the  effective  stiffness  at  joint  a is 

» 

comparted  as  K . + R • 
z-1  z 

4.  As  usual,  the  foregoing  steps  are  repeated  for  a number  of  frequencies 
and  the  natural  frequencies  are  determined  as  those  frequencies  for 
which  the  effective  stiffness  vanishes. 

In  the  foregoing  discussion  it  was  assumed  that  = 0.  Consider 
new  that  support  z is  fixed:  then,  R_  = infinite,  © = 0,  and  M is  finite. 

Zt  Z Z 

But,  since 


n. 


Kl-.G, 


K'z-i  " infinity 


(26b) 
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bnccmss  the  codified  criterion  for  a natural  frequency. 

The  curve  in  Fig.  20  shove  the  vr.iation  of  the  effective  stiffness 

— t 

of  a continuous  beam  as  a function  of  the  frequency  of  vibration.  The 
curve  was  determined  by  the  "direct"  method  end  is  applicable  to  the  partic- 
ular beam  considered  in  Example  2.  The  zeros  of  the  curve  correspond  to  the 
natural  frequencies  of  the  beam;  the  discontinuities  correspond  to  the  natural 
f frequencies  of  the  beam  assuming  that  its  right  end  is  fixed.  Tbs  abscissa 
of  any  other  point  of  the  curve  corresponds  to  the  natural  frequency  of  the 
beam,  provided  its  right  end  is  subjected  to  a restraint,  the  stiffness  of 
which  is  equal  to  the  negative  of  the  value  represented  by  the  ordinate  of 
the  curve. 


It  should  be  pointed  out  again  that  the  main  method,  which  was  pre- 
sented at  the  beginning  of  this  Chapter,  and  the  effective  stiffness  method 
presented  in  the  preceding  paragraphs,  are  fundamentally  alike.  In  the 
former  method,  the  moment  at  joint  z necessary  to  produce  a rotation  of  unit 
amplitude  at  joint  1 is  determined,  while  in  the  latter  method,  the  moment 
at  joint  z necessary  to  produce  a rotation  of  unit  amplitude  at  the  same 
Joint  z is  determined.  The  correspondence  of  the  two  methods  can  be  demon- 
strated further  by  noting  that,  if  each  ordinate  of  the  curve  in  Fig.  18  is 


divided  by  the  rotation  of  the  beam  corresponding  to  that  ordinate,  the 

curve  will  be  transformed  into  that  shown  in  Fig.  20.  For  example,  for 

= 2.40  the  ordinate  of  the  curve  in  Fig.  18  is  17.55  Si  , and  the 

17  55  ^*1  E-i  I. 

corresponding  rotation  ©2  = 6^  = 21.463 • The  ratio  21*463  jT  ‘ = 

E- 1^  x 

-ftT7Q  — m irtain+.l  f»»l  +x>  +.ba  on-rroonAnrflrm  3 £ thS  CUTV5  in 

h ----- 

Fig.  20. 

The  effective  stiffness  method  is  similar  to  Porter's  (18)  and 


Manley's  (19),  (20)  methods  of  determining  natural  frequencies  of  torsional 


vibration  of  shafts,  and  is  similar  elao  to  Lundquist's  stiffness  and  series 
criteria  for  determining  the  critical  buckling  loads  of  structures  (21).  The 
stiffness  criterion  has  been  applied  to  the  determination  of  the  natural  fre- 
quencies of  continuous  bears  previously  by  Mudrak  (1),  (2),  (3)*  However, 
Mudrak’s  method  differs  from  the  procedures  described  in  this  section  both  in 
its  development  and  in  the  form  of  its  application. 

The  Method  of  Three  Moments.  An  alternate  procedure  for  calculating  the 
magnitude  of  the  exciting  moment  M is  provided  by  the  use  of  the  equation  of 
three  moments,  first  applied  to  the  study  of  steady-state  forced  vibrations 
by  W.  Prager  (22).  Humerical  values  of  the  various  coefficients  appearing  in 
these  equations  have  been  published  (23),  (24)  but, unfortunately,  these 
references  are  not  readily  accessible.  In  general,  the  three-moment  equa- 
tion can  be  applied  in  the  same  manner  as  the  three- slope  equation. 

14.  Range  of  Applicability  and  Relative  Merita  of  VarlouB  Procedures. 

As  previously  remarked,  the  moment  distribution  procedure  is  of 
restricted  practical  value.  Convergence  of  this  procedure  can  be  insured 
only  for  vibration  frequencies  which  are  (smaller  that  the  (unknown)  funda- 
mental or  lowest  natural  frequency  of  the  system  considered  (13) • Con- 
sequently, the  method  can,  in  general,  be  used  to  determine  only  the  lowest 
natural  frequency  of  a structure.  AIbo,  it  might  be  important  to  note  that, 
even  for  vibration  frequencies  which  are  below  the  fundamental  natural  fre- 
quency of  a structure,  the  moment  dietri button  procedure  may  be  so  slow  to 
converge  that  it  may  be  necessary  to  carry  out  a large  number  of  distri- 
butions to  affect  a solution.  This  process  may  become  rather  time  consuming, 
especially  whon  applied  to  structures  involving  a large  number  of  members. 

The  "direct"  method  does  not  offer  any  difficulty  of  convergence. 
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It  can,  therefore,  bs  used  to  calculate  the  higher  natural  frequencies  of 
continuous  bean.  •.  In  general,  this  procedure  requires  a much  larger  number 

* 

; of  trials  than  the  main  procedure  of  this  report.  In  addition,  it  cannot  be 

% 

extended  to  continuous  frames  involving  closed  panels.  It  is,  therefore,  of 
restricted  applicability,  too. 

For  continuous  beams  only,  the  choice  between  the  main  method  of 
this  report  and  the  procedure  based  on  the  use  of  the  three  moment  equation 
depends,  to  a large  extent,  on  personal  preference  and  on  one's  familiarity 
with  the  particular  procedure.  One  major  advantage  of  the  use  of  the  three 
slope  equation  is  that  it  gives  a clear  picture  of  the  distortions  which  the 
structure  undergoes  during  vibration.  This  feature  is  particularly  important 

^ because?  in  practice,  it  is  frequently  desirable  to  have  a rapid  means  of 

sketching  the  vibration  configuration  corresponding  to  a given  frequency. 

For  the  analysis  of  continuous  frames,  equations  involving  the  rotation  of 
the  joints  as  unknowns  are  remarkably  better  suited  than  equations  involving 
moments  as  the  redundaAt -quantities.  The  extension  of  the  main  method  of 
this  report  to  the  determination  of  the  natural  frequencies  of  continuous 
frames  without  sidesway  is  presented  in  the  following  Chapter. 


I 
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* 17.  APPLICATION  OF  METHOD  TO  CONTINUOUS  FRAMES 

' WITHOUT  SILESVAI 


15.  General 

This  Chapter  Is  concerned  with  the  determination  of  the  natural 
frequencies  of  flexural  vibration  of  rigid  jointed  plane  frameworks  for 
which  the  joints  do  not  move.  The  extension  of  the  method  to  some  rela- 
tively simple  frames  with  sidesway  will  be  presented  in  Chapter  VI. 

The  frames  considered  may  have  any  number  of  members  of  arbitrary 
length*  the  mass  per  unit  of  length  and  the  flexural  rigidity  of  cross  section 
of  the  members  may  differ  from  one  member  to  the  other,  but  in  any  one  member, 
these  quantities  are  assumed  to  remain  constant.  The  simplifying  assump- 

i* 

tione  made  in  the  analysis  are  as  follows*  The  vibrations  are  assumed  to 
take  place  in  the  plane  of  the  framework.  The  change  in  length  of  the  mem- 
bers due  to  axial  deformation,  and  the  effect  of  the  axial  forces  on  the 
bending  moment  in  the  members  are  neglected.  In  addition,  no  account  is 
taken  of  the  influence  of  axial  vibrations.  As  before,  the  cross  sectional 
dimensions  of  the  members  are  considered  to  be  small  in  comparison  to  their 
length,  so  that  the  effects  of  shearing  deformation  and  rotatory  inertia  may 
be  neglected. 

16.  Basic  Relations 

' Figure  21  shows  s members  of  a structure  rigidly  connected  at  their 

common  intersection  o.  The  far  ends  of  the  members  are  assumed  to  be  fixed 
• against  translation,  but  free  to  rotate  faubject  to  the  restraint  imposed  by 

the  adjoining  members.  Assume  that  the  structure  is  in  a steady-state  forced 
» vibration  under  the  action  of  some  exciting  moment  applied  at  a joint  differ- 

ent from  joint  o. 


_ ,«k- 


Let  denote  the  amplitude  of  rotation  at  end  o of  member  oj,  and 
©jq  denote  the  amplitude  of  rotation  at  end  j of  the  aarae  member.  Similarly, 
let  Mqj  and  Mjq  be  the  corresponding  moment  amplitudes  at  the  same  ends* 

Since  all  members  are  rigidly  connected  st  their  joints, 


0OI  " 6ot  “ 0oj  = * 609  - 60  , (2 8a 

and  6jc  - Qj  . (281 

Furthermore,  since  no  external  moment  acts  at  joint  o, 

/%  -n„ ,+(%,+ *ti. j* ♦«„  - Xr%j  - o . (29) 

J-i 

The  moment  Mq  ^ may  be  expressed  in  termB  of  the  end  rotations  of 
member  oj  by  the  relation 

n.j  - KoA  f (*K)°j  eJ  ■ (30) 

Substituting  this  expression  into  Eq.  (29),  one  obtains  Eq.  (31a) 

Z Kojdo  + Z (*Khj  0jm  0 - (31a 

which  expresses  the  conditions  of  both  equilibrium  and  continuity  for  joint 
o of  the  structure*  If  only  two  members  meet  at  joint  o,  Eq.  (31a)  reduces 
to  Eq.  (19a)  for  continuous  beams. 

If  the  degree  of  restraint  at  the  far  ends  of  the  members  meeting 
at  a joint  are  known,  it  is  convenient  to  use  effective  stiffnesses.  Assume 
that  the  restraints  at  ends  1 and  2 of  the  portion  of  the  structure  shown  in 
Fig.  21  are  known.  Let  Kq^  and  represent  the  effective  stiffness  of 
members  ol  and  o2.  Then,  Eq.  (31a)  may  be  written  as 

(K‘„  + K„)  6 „ + t K< ,j  do  + t (*M.j  dj-0  . (3it 
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Equetions  (30)  and  (31)  are  the  only  two  relations  needed  in  the  analysis 
of  frames  without  sidesway. 


ill  Qpoa  Frames 


The  frames  considered  in  thla  section  do  not  involve  any  closed 
panels  and  have  known  conditions  of  restraint  at  all  external  terminals* 

It  is  assumed  that  the  joints  of  the  frame  do  not  translate* 

Simple  L-frames  and  portal  frames , such  as  those  shown  in  Fig*  22, 
act  as  continuous  beams  on  rigid  supports*  Their  natural  frequencies  can 
therefore  be  calculated  by  the  procedure  outlined  in  Section  10  of  the  pre- 
ceding Chapter. 

When  applied  to  the  analysis  of  continuous  frames , such  as  those 
shown  in  Fig*  23,  this  procedure  will,  in  general,  reveal  only  a portion  of 
the  natural  frequencies  of  the  frame  considered.  The  failure  of  the  procedure 
to  identify  the  complete  set  of  natural  frequencies  results  from  the  fact 
that,  for  certain  natural  frequencies,  only  a portion  of  the  frame  may  vibrate 
with  finite  amplitudes  while  the  rest  may  remain  stationary.  The  natural 
frequencies  corresponding  to  these  modes,  which  will  be  referred  to  as  "modes 
of  partial  vibration",  must  be  determined  by  a supplementary  procedure. 

Consider  any  of  the  frames  shown  in  Fig.  23*  Let  1 denote  the  joint 
of  the  frame  at  the  extreme  left  terminal  and  £ denote  the  joint  at  the  ex- 
treme right  terminal.  Without  loss  of  generality,  it  may  be  assumed  that 
joint  z is  either  hinged  or  elastically  restrained.  A fixed  end  maybe 
handled  in  the  manner  described  in  Illustrative  Example  1.  Assume  that  the 
structure  is  in  a steady-state  forced  vibration  under  the  action  of  an  ex- 
citing couple  Ma  applied  at  joint  s.  The  amplitude  of  the  slope  or  of  the 
bending  moment  at  joint  1 is  assumed  to  have  some  fixed  value.  If  the  joint 
is  hinged  or  elastically  restrained,  the  amplitude  of  slope  is  taken  equal 
to  unity.  If  the  joint  is  fixed,  the  amplitude  of  bending  moment  is  taken 
equal  to  unity  instead.  For  an  assumed  frequency  of  vibration,  it  is 
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generally  possible  to  calculate  the  magnitude  of  the  exciting  moment  Mg»  in 
a manner  entirely  analogous  to  that  used  for  continuous  beams,  by  working 
progressively  from  one  end  of  the  frame  to  the  other.  2y  repeating  this 
procedure  for  several  frequencies  of  vibration,  the  magnitude  of  M may  be 
plotted  as  a function  of  the  frequency.  All  frequencies  for  which  the 
magnitude  of  the  exciting  moment  vanishes  are  natural  frequencies  of  the 
frame. 

The  natural  frequencies  determined  by  the  previous  procedure  may 
not  represent  the  complete  set  of  natural  frequencies  of  the  frame*  The 
procedure  is  based  on  the  assumption  that  the  amplitude  of  slope  or  bending 
moment  at  joint  1 is  finite.  For  continuous  beams  this  condition  is  satis- 
fied for  all  non-trivial  natural  frequencies.  For  continuous  frames,  how- 
ever, bar  1-2  may  be  still  even  though  the  rest  of  the  structure,  or  some 
portion  of  it,  vibrates  with  finite  amplitudes.  Obviously  then,  the  pro- 
cedure fails  to  reveal  those  natural  frequencies  for  which  bar  1-2  remains 
still.  A second  assumption  implicit  in  the  procedure  described  is  that  the 
rotation  of  joint  z is  finite  for  the  natural  frequencies  to  be  determined. 

For  continuous  framss,  this  condition  is  not  satisfied  always.  Therefore, 
the  procedure  fails  also  to  reveal  the  natural  frequencies  for  which  the  bar 
meeting  at  joint  z is  stationary. 

Figure  24  presents  several  natural  modes  of  vibration  for  which 
either  bar  1-2  or  the  bar  meeting  at  Joint  z is  stationary.  The  modes  are 
applicable  to  the  particular  atructure8  8hown  in  Fig.  23  and  can,  of  course, 
exist  only  if  the  dimensions  of  the  various  members  composing  these  structures 
satisfy  certain  definite  relations.  It  should  be  emphasized  that  natural  modes 
of  partial  vibration  are  peculiar  to  frames  and  cannot  exist  in  the  case  of 
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The  technioue  for  determining  the  natural  frequencies  for  which 
either  bar  1-2  or  the  bar  meeting  at  joint  z is  stationary  consists  of 

(a)  calculating  the  frequencies  for  which  these  conditions  esn  occur,  and 

(b)  ascertaining  whether  or  not  these  frequencies  are  natural  frequencies  of 
the  system.  The  details  of  this  supplementary  technique  will  be  explained 
in  the  examples  to  be  presented. 


8 
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18^  Illustrative  Examples. 

Example  3.  The  simple  frame  shown  in  Fig.  25a  has  been  selected 
for  analysis.  To  illustrate  several  features  of  the  method,  the  bars 
identified  by  (1)  are  taken  identical  whi?*>  the  bar  designated  by  (2)  iB 
considered  to  have  such  dimensions  that 


and 


g2Z2  = , Lj  b 0.801^ 

h h 


K = f'lfn  = 0.6302  A, 


4.730 

The  aubscripts  1 and  2 refer  to  bars  (1)  and  (2),  respectively* 

For  the  sake  of  brevity,  only  one  cycle  of  the  procedure  is  pre- 
sented. The  computations  are  given  for  a value  of  A,  = 3*30;  this  corresponds 
to  a value  of  A2  = 2.74.  The  appropriate  values  of  K and  kK  are  obtained 
from  Table  I in  Appendix  A. 


LI  E.L. 

for  bars  (1)  Kx  = 2.5720  and  (kK)1  - 3.1375 


for  bar  (2)  ^ = 3.4051  and  (kX)2  = 2.4589 

The  data  necessary  for  the  analysis  are  compiled  on  the  disgram  in  Fig.  25b. 
The  number  in  parentheses  opposite  each  joint  gives  the  sum  of  the  stiffnesses 
of  the  members  meeting  at  that  Joint.  The  parenthesized  number  at  the  middle 
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of  each  member  gives  the  product  of  the  stiffness  and  the  carry-over  factor 

Vl 


for  the  member.  Both  quantities  are  expressed  in  terms  of  ~y~~  • The  numbers 
without  parentheses  denote  the  rotations  of  the  various  joints.  These  rota- 
tions are  evaluated  ir  the  manner  described  below,  and  they  are  recorded  on 
the  diagram  as  they  are  computed, 

Vi 

The  procedure  ia  started  by  taking  = 1.00  • Then,  ©,- 

computed  by  application  to  joint  1 of  Eq.  (30),  as 

*2  ' 33375  * °'31873 

The  rotations  ©^  and  ©^  are  determined  successively  by  application  of  Eq. 

(31a)  to  joints  2 and  4* 

04  * - MW  (0)..  -0.52256 

6 - 6.5491  (-.52256)  + 3.1375  f .31873)  _ 1_u01 

5 2.4589 

The  magnitude  of  the  exciting  moment  is  determined  from  Eq.  (30),  as 

B,Il  *,11  R.L 

Me  = (1.4101)3.4051  -«  + (-.52256)2.4589  -f-4  = 3*517 
•5  U n n 

It  should  be  noted  that  each  of  the  foregoing  operations  can  be  carried  out 
with  a single  set-up  on  a desk  calculstor. 

Repeating  this  procedure  for  several  values  of  \ , the  curve 
shown  in  Fig.  26  was  obtained.  The  A,  values  corresponding  to  the  natural 
frequencies  are  recorded  on  the  figure.  The  corresponding  natural  modes  of 
free  vibration  are  given  in  Figs.  27a  to  27f* 

The  foregoing  procedure  is  based  on  the  assumption  that  the  amplitudes 
of  both  the  bending  moment  at  joint  1 and  of  the  rotation  at  joint  z - 5 are 
finite.  To  obtain  the  natural  frequencies  for  which  = 0,  the  following 
reasoning  is  used.  In  order  for  to  be  equal  to  zero,  bar  1-2  must  be 


fljfft*'  Mri  mw**"**.* 


u 

stationary;  then*  the  rotation  of  joint  2 and  the  internal  bending  moment  at 
the  joint  are  also  equal  to  zero.  This  condition  requires  that  bar  2-4  be 
stationary  and  that  S4  = M42  = 0.  But  with  joint  4 remaining  fixed  against 
rotation,  bar  3-4  can  oscillate  freely  only  for  frequencies  represented  by 
values  of 

A,  = 4.730,  7.853*  (32) 

During  a natural  frequency,  every  member  of  the  structure  vibrates  with 
the  same  frequency.  Therefore,  in  order  for  these  frequencies  to  be  natural 
frequencies  of  the  entire  frame,  they  must  be  a!) so  natural  frequencies  of  the 
remaining  portion  of  the  frame  (bar  4-5,  in  thi  - particular  case).  The 
natural  frequencies  of  bar  4-5,  considering  that  its  left  end  is  fixed,  are 

Az  = 3.927  , 7.069 

these  correspond  to  values  of 

A(  ~ 4.730,  8.514,  ...... 

Comparing  the  latter  values  with  those  given  in  Eq.  (32),  one  concludes  that, 
within  the  range  of  frequencies  considered  in  Pig.  26,  A,  = 7.853  and 
A,  = 8*514  do  not  represent  a natural  frequency,  while  X>  = 4.730  does. 

The  natural  mode  of  vibration  for  A,  = 4.730  is  shown  by  the  solid  curve  in 
Fig.  27g. 

The  natural  frequencies,  if  any,  for  which  6^  = 0,  are  determined 
in  t similar  manner.  0^  can  be  equal  to  zero  only  if  bar  4-5  is  stationary. 

9 

Under  this  condition;  = 0;  then,  bar  3-4  say  vibrate  freely  only  at 

frequencies  r ©presented  by  the  A,  values  given  in  Eq.  (32).  Since  members 
1-2  and  2-4  of  the  remaining  portion  of  the  frame  are  identical  to  member 
3-4,  each  of  these  members  can  vibrate  with  its  ends  fixed  for  the  same  fre- 
quency; therefore,  the  A.  values  given  in  Eq.  (32)  correspond  to  natural 
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frequencies  of  the  frame . The  natural  node  corresponding  to  X,  ~ 4*730  is 
shewn  by  the  dotted  curve  in  Fig.  27g,  while  the  mode  corresponding  to 
A,  = 7.853  is  shown  in  Fig.  27h. 

It  should  be  observed  that  the  natural  nodes  shown  in  Figs*  27c  and 
27  g can  exist  for  the  same  frequency.  Of  these  three  modes,  however,  only 
the  two  are  independent;  the  third  is  a linear  combination  of  the  other  two. 

In  fact,  from  any  two  of  these  three  nodea,  one  oan  obtain  an  infinite  number 
of  combination  modes. 

Within  the  range  of  frequencies  considered  in  Fig.  26,  the  complete 
set  of  Aj  values  corresponding  to  natural  frequencies  is 

3*59,  4.22,  4.73  (double),  6.80,  7.44»  7.85,  8.35 

More  involved  frames  may  be  handled  by  the  same  procedure.  The 
technique  for  obtaining  the  natural  frequencies  corresponding  to  modes  of 
partial  vibration  ie  illustrated  further  by  the  examples  given  in  Section  21. 

Example  4.  A sketch  of  the  frame  considered  is  shown  in  Fig.  28. 

This  frame  ia  similar  to  that  analysed  in  the  preceding  example.  In  this 
case,  the  dimensions  of  the  structure  are  assumed  to  be  such  that 

_ Vs  _ iis  i V4 

' ^ ~ L3  ~ L4 
and  Xz-  0.75 A. , A,=  1.30A,,  \=  0.90 A,  . 

To  determine  the  natural  frequencies  of  this  frame,  we  proceed  in 
the  usual  manner  and  plot  the  magnitude  of  the  exciting  moment  at  joint  5 as 
a function  of  the  assumed  frequency  of  vibration.  The  curve  in  Fig.  28 
summarizes  the  results  obtained.  It  should  be  noted  that' this  curve,  un- 
like that  shown  in  Fig.  26,  1b  not  continuous.  The  values  of  A,  corresponding 
to  tha  first  few  natural  frequencies  are  recorded  on  Fig.  28.  It  can  be  shown 
that,  within  the  range  of  the  frequencies  considered,  there  are  no  natural 
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frequencies  corresponding  to  sodas  of  partial  vibration* 


19 » Closed  Erases* 

The  application  of  the  aethod  to  f rases  involving  closed  panels* 
such  as  those  sbovn  in  Fig.  29*  is  described  in  this  section  by  reference 
to  a simple  example.  The  hypothetical  two  celled  rectangular  frame  shown 
in  Fig.  29a  is  selected  for  this  ..purpose. 

The  first  step  in  the  analysis  is  to  assume  that  the  frame  is  cut 
at  some  convenient  joint.  In  the  example  considered*  the  cut  is  introduced 
at  joint  6.  Next,  it  is  assumed  that  the  structure  undergoes  a steady-state 
forced  vibration  with  finite  amplitudes  and  known  frequency.  The  vibration 
is  assumed  to  be  maintained  by  an  exciting  couple  applied  at  the  cut  joint 
(Joint  6).  If  the  frequency  of  vibration  is  equal  to  the  natural  frequency 
of  the  frame,  the  magnitude  of  the  exciting  moment  must  vanish  and  the 
amplitudes  of  slope  on  either  side  of  the  cut  must  be  equal. 

Fcr  an  assumed  frequency,  the  discontinuity  of  slope  and  the  magni- 
tude of  the  exciting  moment  may  be  determined  in  the  same  way  as  for  open 
frames,  by  working  progressively  from  joint  to  joint  across  the  structure. 

For  the  frame  considered,  Bq»  (31a)  is  first  applied  to  joint  1. 

It  is  noted  that  the  resulting  expression  involves  three  unknowns i ©2, 

and  9y  Therefore,  9^  an(*  ®3  cannot  be  solved  directly  in  terms  of  ©^  alone, 
as  it  was  possible  in  the  case  of  continuous  beams  and  continuous  open  frames. 
Instead,  it  is  necessary  to  express  ©^  in  terms  of  both  ©^  and  Gy  Next, 

Sq.  (31a)  is  applied  to  joint  2.  and  is  determined  i n terms  of  who  owxd 
two  permanent  unknowns  ©^  and  ©2*  By  successive  applications  of  the  same 
equation  to  joints  3»  4,  and  5»  the  rotations  6^,  6^,  and  645  may  also  be 
determined  in  terms  of  and  ©2*  having  computed  the  rotations  of  all 


Joints,  the  exciting  moment  at  Joint  6 may  also  be  expressed  in  terms  of 


47 


the  two  permanent  unknowns  sad  In  the  application  of  tale  technique* 
consecutive  joints  must  be  selected  in  such  an  order  that,  when  Bq.  (31a) 
is  applied  to  a joint,  the  resulting  expression  Involves  only  one  new  unknown* 
This  technique  is  an  adaption  of  Wilbur's  aoheroe  (2$)  of  solving  the  set  of 
simultaneous  equations  resulting  from  the  use  of  the  elope  deflection  equa^ 
tion. 


The  joint  for  which  the  rotation  ia  evaluated  last  in  this  procedure 
(joint  6 in  this  esse),  is  the  ons  st  which  tha  structure  ie  generally  as- 
sumed to  be  out*  The  number  of  members  meeting  at  this  joint  must  be  equal 
to  the  number  of  permanent  unknowns  used* 

At  a natural  frequency 


and 


£«♦  ~ 9gg  m 0 , 

• 0 . 

In  terms  of  the  two  permanent  unknowns  9^  and  these  conditions  may  be 
written  •• 


(33«) 


C„  S,  + Cit&g  m 0 , 


(?Sb) 


CttSi  f CnOt  “ 0 . 

where  the  o'e  are  constants,  the  magnitudes  of  whieh  depend  on  the  assumed 
frequency  of  vibration  and  on  the  characteristics  of  ell  the  members  in  the 
structure*  Equation  (33b)  represents  a aet  of  linear  and  homogeneous  equa- 
tions for  the  unknowns  Sv  and  these  equation*  have  a solution  different 
from  zero  when  the  determinant  of  their  ooeffioiant  is  sero, 


c„ 

C,? 

ca 

(34) 


This  criterion  of  vanishing  determinant  will  fail  to  reveal 
(a)  the  natural  frequencies  for  vhieh  6^  and  are  aimultsneoualy 

equal  to  zero*  and 


. vw*« 


4fi 

(b)  the  natural  frequencies  for  which  the  bers  Esetlsg  4t  the  cirt-  joint 
(joint  6)  are  stationary. 

The  natural  frequencies  corresponding  to  the  foregoing  two  conditions 
B87  be  calculated  by  the  supplementary  technique  described  in  connection  with 
continuous  open  frames c 

The  single-bay  multi-story  frame  shorn  in  fig*  29b  may  be  handled 
cy  the  ease  procedure.  The  rotations  of  its  joints  a ay  be  expressed  in  tanas 
of  9^  and  6.,)  the  cut  may  be  introduced  at  joint  9 or  10*  For  the  analysis 
of  the  two-bay  multi -story  frame  shown  in  Fig*  29c,  one  needs  to  take  three 
quantities,  say  9^,  and  9^»  98  permanent  unknowns*  The  cut  must  be 
introduced  at  joint  14*  The  conditions  of  continuity  and  equilibrium  for 
this  joint  may  then  be  expressed  as 

&*.n  ~ m C„&i  * Cn ft  * C, , “ 0 , 

~ &KJ9  “ C»Q  tcs^f  ^ 9%  m 0 > 
m CjtfQ  + C||0»  * CjS  0j  " 0 , 

where  the  £*a  are  numerloal  constants*  The  criterion  for  a natural  fre- 
quency is 
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It  is  seen  that,  for  this  problem,  it  becomes  necessary  to  evaluate  a third 
order  determinant.  For  the  generel  came,  the  order  of  the  determinant  is 
equal  to  the  number  of  the  permanent  unknowns  that  must  be  used  in  evalu- 
ating the  rotations  of  tha  joints* 

The  computational  work  required  to  oaloulate  the  exciting  moment 
and  the  discontinuity  of  slope  at  the  out,  may  get  fairly  involved,  particu- 
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larly  if  sore  than  two  quantities  must  be  used  as  permanent  unknowns*  One 


may  simplify  this  work  considerably  by  carrying  out  the  computations  in 


parts s Consider  egsin  the  frame  shown  in  Fig*  29«*  First,  assume  that 


©^  = 1*00  and  that  ©^  = 0*  Calculate  the  rotations  of  the  joints  in  the 


usual  manner,  and  designate  them  by  © . Calculate  also  the  magnitude  of 


the  exciting  moment  at  the  cut  .joint*  Designate  this  by  M . Next,  assume 


that  ©1  = C and  that  ©^  = 1.00,  and  calculate  the  corresponding  rotations 


and  moment.  Denote  these  by  © and  M , respectively.  Then,  the  actual 


rotation  at  a joint  j is 


®J  = °j°l  + ®j°2  » 


(35a) 


and  the  total  exciting  moment  at  the  cut  joint,  say  joint  j , is 


5J  = «A  * s>2  • 


(35b) 


The  direct  combination  of  these  partial  effects  is  justified  by  the  fact 


that  the  differential  equation  for  steady^state  forced  vibration  is  linear. 


The  procedure  for  determining  the  natural  frequencies  of  con- 


tinuous frames  Involving  closed  panels  may  be  outlined  as  follows i 


1.  For  some  member  of  the  frame,  say  member  r,  assume  e value  of 


4/  ffUUJ* 


A rmJ  £* ir  Lr  • This  is  equivalent  to  assuming  a frequency 


of  vibration  to  • 


..._2.  From  Eq.  (23)  compute  the  A values  for  the  remaining  members  of 


the  frame* 


3*  From  Table  I in  Appendix  A,  calculate  the  appropriate  values  of  K 


and  kg  for  each  member.  These  values, as  obtained  from  Table  I, 
are  expressed  in  terms  of  -S-*  , where  the  subscript  j_  refers  to 
the  particular  m«aber  considered. 
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* 


4*  repress  tbs  quantities  determined  in  step  (3)  in  terms  of  the 
rj?  of  the  reference  member  r,  by  multiplying  them  by  the  dimen- 
sionless coefficient  (Bq.  24). 

5 <■  Compute  the  sum  of  the  stiffnesses  of  the  members  meeting  at  the 
various  joints  and  record  these  values  on  a diagram  of  the  frame* 
On  the  same  diagram  record  also  the  product  of  the  stiffness  end 
the  carry-over  factor  for  each  member.  These  quantities  must  bo 
expressed  in  terms  of  the  of  the  same  reference  member  r.  A 
convenient  scheme  for  arranging  the  computations  is  shown  in  the 
illustrative  examples  presented  in  the  next  section. 

6.  Choose  the  unknowns  in  terms  of  which  the  distortions  of  the  frame 
will  be  expressed*  In  general*  the  number  of  unknowns  that  must 
be  selected  is  equal  to  the  number  of  the  main  longitudinal 
members  in  the  frame.  For  the  frame  shown  in  Fig*  29a,  one  may 
take  and  88  the  two  permanent  unknowns* 

7*  Consider  the  first  of  these  quantities  equal  to  unity  and  the 
other  equal  to  zero  (6-  =1,00  and  9^  = 0),  Working  across  the 
structure,  as  described  in  the  preceding  section*  compute  the 

i 

rotations  of  the  joints.  Denote  these  by  6 * Compute  also  the 
exciting  moment  and  denote  it  by  m'* 

8*  Repeat  step  (7)*  taking  the  second  quantity  equal  to  unity  and 
the  first  equal  to  zero  (0^  = 0 and  9^  = 1*00).  Denote  the 

K — « 

resulting  rotations  by  9 and  the  exciting  moment  by  M * In 

» 4 ann  /rj\  e *\A  ( &\  wu  ■+  V\a  waraaia^  o • TBBWtT  + A WA  a AM 

gWUVA  MA|  m0  VW \ I / w»t«*  \w/  ****•—  V Arw  a V|fWM  vw«*  MV  — «*I »iy  VAMW  MW 

there  are  permanent  unknowns* 

9*  Determine  the  total  discontinuity  of  slope  and  the  total  ex- 
citing moment  and  set  each  expression  equal  to  zero*  For  the 


■4 
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frame  considered,  these  expressions  will  be 


re«-sye,  f<e„-^s)e,  -o. 
r%e,  » n'et-o. 

10.  Evaluate  the  determinant  of  the  coefficients  of  and  in  the 
expressions  of  Eq.  (33c). 

11.  Repeat  steps  (l)  to  (10)  for  different  values  of  Ar  • 

12.  Plot  the  variation  of  the  determinant  evaluated  in  step  (10)  as  a 
function  of  the  Ar  values.  The  frequencies  for  which  the  deter- 
minant becomes  equal  to  zero  are  natural  frequencies  of  the  frame. 

The  foregoing  procedure  fails  to  reveal  the  natural  frequencies 
for  which  the  permanent  unknowns  (0^  and  ©^)  are  simultaneously  equal  to 
zero.  In  addition,  it  fails  to  reveal  the  natural  frequencies  for  which 
the  members  meeting  at  the  cut  joint  are  stationary.  A supplementary 
procedure  for  determining  these  natural  frequencies  has  been  described, 
and  its  details  are  illustrated  further  in  the  two  numerical  examples 
that  follow. 


(33c) 


21.  Illustrative  Examples . 

Example  5-  The  structure  considered  is  shown  in  Fig.  30a • All 
members  are  assumed  to  be  uniform  and  identical  to  each  other.  It  is  de- 
sired to  calculate  the  natural  frequencies  and  the  corresponding  natural 
modes  of  vibration  of  this  frame  for  a range  of  values  of  A less  than  6.50. 
Since  all  members  are  identical,  rather  than  repeating  the  procedure  out- 
lined in  Section  20  for  each  assumed  frequency,  it  is  more  convenient  to 
derive  a general  expression  for  the  criterion  for  a natural  frequency  and 
determine  directly  from  this  expression  the  desired  quantities.  This  is 
similar  to  what  was  done  in  Section  12  for  Example .1, 
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Let  X be  the  flexural  stiffness  and  k the  flexural  carry-over 
factor  for  each  member  of  the  frame.  These  quantities  are,  of  course, 
functions  of  the  parameter  A . The  sum  of  the  stiffnesses  for  the  vari- 
ous joints  and  the  product  of  the  stiffnees  and  the  carry-over  factor  for 
each  member  of  the  frame  are  shown  in  Fig.  30b.  The  rotations  cf  the 
joints  will  be  expressed,  not  in  terms  of  6^  and  ©2>as  ^ wefl  BU€SeB'te<* 
in  the  preceding  discussion,  but  rather,  in  terms  of  0^  and  the  internal 
bending  moment  M^.  The  reason  for  this  choice  will  become  apparent 
shortly.  The  frame  ia  assumed  to  be  cut  at  joint  6*  First,  it  is  assumed 
that  Mj^  = Mj.  = 1.00  and  ^ = 0. 

(1)  and  (2),  one  obtains 


Similarly,  applying  Eq.  (31a)  successively  to  joints  2,  3,  4»  and  5»  one 
obtains 

%-■ -«[$&)**(»]+  *K-  W- 

0-5“  ~k[3(*k)  + *W>  + *(t?k)] * *K  ' ~ *‘K  ' 

[3(lfc)+  *('**) f *(*%)]+  *K  * 

4»-  ~K  [ 2fw/t  ■ "w 


Applying  Bq.  (30)  to  ends  1 of  members 
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The  discontinuity  of  slope  at  joint  6 is 


004  - eL 


and  the  exciting  moment  at  joint  6 is 
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Ndztf  it  is  assumed  that  = 0 and  6^  - 1.00.  The  rotations  6 
are  obtained  in  a similar  manner.  The  results  are 

nl  I 

% ‘ TT' 

d’  * ' 

% ■ -*[*(-&*  +]*  *K‘^- 
C' -K[3jLj£-  * *(--£■)  * *(-&)]  - = - 

C-^-F  + *(-£)]+  *K 

The  discontinuity  of  slope  and  the  exciting  moment  at  joint  6 are 

0«-«5  -^r- 

/%•  - 1-  *-£rj  - %(l31t-IZ-#) e,  . 

The  total  discontinuity  of  slope  and  the  total  exciting  moment 
at  joint  6 are 


a 


~%r(4-+ *x)n,  +Jfeo3i?-iz-ir)e< 


(36) 


I • 


The  determinant  of  the  coefficients  of  and  is 

A*(n,,6,)  ■ 


144 

IT 


it-  41  It + 184 

It 


(37) 


The  curve  in  Fig.  30d  shows  the  variation  of  this  determinant  with 
X • The  values  of  k and  K corresponding  to  the  various  values  of  X were 
obtained  from  Table  I in  Appendix  A.  The  aero  intercepts  of  this  curve  rep- 
resent natural  frequencies  of  vibration.  These  values  are  indicated  on  the 
figure. 


The  vibration  modes  corresponding  to  a natural  frequency  are 
determined  as  follows*  First,  the  relationship  between  and  is  deter- 
mined by  setting  either  u-t  the  expressions  in  Eq.  (36)  equal  to  zero.  Then, 
the  rotations  of  the  joints  are  evaluated,  in  terms  of  either  6^  or  M^,  by 
use  of  the  relation 

ej  = * Vi 

In  this  particular  case,  the  rotations  were  expressed  in  terms  of  6^.  The 
results  are  summarized  in  the  following  table. 


Value  of 
An 

V*L 



8 A 

w 

8A 

V81 

7T 

-1.00 

1.00 

1.00 

-1.00 

3:556 

■5SS 

0 

0 

-1.00 

1.00 

3.805 

1.00 

-1.333 

-1.333 

1.00 

1.00 

4.048 

-1.00 

1.333 

1.00 

-1.00 

1.00 

0 

0 

-1.00 

-1.00 

2tt 

1.00 

1.00 

1.00 

1.00 

1.00 

It  should  be  noted  that  the  value  of  AN  = 4.730  is  not  included  in  this 
table.  For  this  value,  aa  it  will  follow  from  the  discussion  of  the  succeed- 
ing paragraphs,  there  is  an  infinite  number  of  possible  natural  modes. 

Such  modes  cannot  be  determined  by  the  procedure  described.  From  the 
rotations  in  the  above  table,  the  deflections  at  the  interior  points  of  the 
members  may  be  obtained  by  use  of  the  influence  coefficients  given  in  Table 
II  of  Appendix  A.  The  natural  modes  corresponding  to  the  natural  frequen- 
cies determined  in  Fig.  30d,  are  shown  in  Fig.  31a  through  31g. 

The  next  step  in  the  solution  is  the  determination  of  the  natural 
frequencies  for  which  ©^  and  vanish  simultaneously.  ©^  ®nd  IQay  6e 
equal  to  zero  only  if  both  bars  (1)  and  (2)  are  stationary.  Under  this 
condition,  bur  (3)  also  is  stationary  and  joints  3 and  4 remain  fixed  against 


rotation  ae  shown  in  Fig.  30c.  For  the  range  of  A values  considered,  bar 
(i  ) can  execute  free  vibrations  with  fixed,  eade  only  at  a frequency  represent- 
ed by 

A = 4.730. 

Since  all  members  of  panel  3-^-6- 5 are  identical,  each  member  can  vibrate 
at  this  frequency  with  fixed  endsj  therefore,  A = 4.730  corresponds  to  a 
natural  frequency  of  the  frame.  The  vibration  mode  corresponding  to  this 
naturel  frequency  is  shown  in  Fig.  31h. 

The  concluding  step  in  the  solution  of  the  problem  under  con- 
sideration ie  the  determination  of  the  natural  frequencies  for  which  the 
members  meeting  at  a joint  6 remain  etationejy.  Proceeding  in  the  manner 
described  in  the  preceding  case,  it  can  be  shown  that,  within  the  range  of 
A values  considered,  A = 4.730  represents  a natural  frequency  for  which 
barB  (6)  and  (7)  are  stationary.  The  vibration  mode  corresponding  to  thiB 
natural  frequency  ie  shown  in  Fig.  31i. 

It  Bhould  be  observed  that  the  natural  modes  shown  in  Figs.  31f » 

31h,  and  31i  can  exist  for  the  same  frequency  of  vibration.  However,  of 
theee  three  modes  only  two  are  independent.  The  third  ie  a linear  combin- 
ation of  the  other  two. 

Had  the  distortions  of  the  frame  and  the  exciting  moment  been 
expressed  in  terms  of  ©^  Bn<l  -ss^a-  terms  of  ©^  and  the  baBic 

procedure  would  have  failed  to  reveal  the  natural  frequencies  for  which 
©^  = = 0.  other  words,  the  curve  in  Fig.  30d  would  not  have  inter- 

sected the  A-axia  at  A = 4.730.  Of  course,  this  natural  frequency  would 
have  been  determined  by  the  supplementary  technique. 

Example  6.  As  a further  illustration  of  the  application  of  the 
method  to  continuous  closed  frames,  the  symmetrical  frame  shown  in  Fig.  32a 
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is  selected  for  analysis.  All  columns  are  considered  hinged  at  the  base* 
It  is  desired  to  calculate  natural  frequencies  corresponding  only  to 
symmetrical  modes  of  vibration.  Since  for  symmetrical  vibrations  the 
joints  of  the  frame,  even  though  free  from  external  restraining  forces, 
do  not-  translate;  the  method  of  this  Chapter  is  directly  applicable  to  the 
problem  considered.  The  effect  of  symmetry  is  taken  into  account  by  using 
for  the  girders  of  the  central  bay  modified  stiffnesses  K**  instead  of  the 
usual  stiffnesses  K, 

The  dimensions  of  the  frame  and  some  additional  data  pertinent 
in  the  analysis  are  assembled  in  the  table  below.  Member  (2)  of  the  frame 


is  taken  as  the  reference  member.  Columns  (5)  and  (6)  of  this  table  give. 


(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Member 

h 

% 

M 

l4 

-J. 

L 

r 

r1  = 00(4) 
V 

*4 

1 

0.50 

0.6’ 

0.9294 

1.00 

0.9294 

0.67 

2.04 

2=r 

1.00  - 

1.00 

j.  * 000 

1.00 

1.000 

1.00 

2.20 

3 

3.00 

1.75 

1.H4 

1.50 

1.716 

1.1667 

3.78 

MM 

3.00 

1.50 

1.125 

1.338 

1.3333 

2.94 

El 

0.40 

0.45 

0.97±: 

0.75 

0.7282 

0.600 

1.60 

6 

0.50 

0.625 

0.9457 

0.75 

0.7093 

0.83333 

1.56 

mm 

2.00 

1.25 

1.125 

1.  '0 

1.687 

0.83333 

3.71 

B 

2.00 

1.25 

1.125 

1.12. 

1.265 

1.13.11 

2.78 

respectively,  the  ratio  of  ^j/A z (Eq*  23;  end  the  ^ values  (Eq.  24)  for 
each  member  of  the  frame.  For  the  sake  of  brevity,  only  one  cycle  of  the 


procedure  is  presented.  The  computations  are  given  for  a value  of  ^ = 2.20$ 
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The  cK  values  of  the  various  members  are  recorded  also  on  the  dia- 
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gram  in  Fig.  32b.  The  numbers  without  parentheses  above  these  values  give 
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the  stiff ness  of  the  various  members,  while  the  parenthesized  numbers  below 
the  « values  give  the  product  of  the  stiffness  and  the  carry-over  factor  of 
the  members  s These  quantities  are  expressed  in  terms  of  the  vdf  the  mem- 
ber to  which  they  refer,  and  they  are  obtained  directly  from  Table  I in 
Appendix  A* 

On  the  diagram  in  Fig.  32c,  the  number  in  parenthesis  opposite 
each  joint  represents  the  total  stiffness  of  the  members  framing  into  that 
joint.  These  stiffnesses  are  expressed  in  terms  of  the  ^ of  the  reference 
member  (2).  Thus,  the  total  stiffness  of  the  members  connecting  at  joint  4 is 
K4  = 3.7675x1.00  + 0.91884x1.1667  + 3.9430x0.83333  + 0.54510x1.3333  = 8.8521  . 
She  parenthesised  number  at  the  middle  of  each  member  represents  the  product 
of  the  stiffness  and  the  carry-over  factor  for  that  members  these  values 
are  also  expressed  in  terms  of  the  IP  of  the  reference  member  (2). 

Th6  rotations  of  the  joints  are  expressed  in  terms  of  0^  and 
The  frame  is  assumed  to  be  cut  at  joint  5 (and  joint  5*).  First,  consider 
that 

©^  = 1.00  and  ©2  = 0. 

Applying  Bq.  (30)  to  joints  1 and  2,  one  obtains 


2.5661  x 1.00 
1.4262 


3.7675  x 0 
2.1764 


-1.7993  , 


0 . 


As  usual,  these  values  are  recorded  on  the  diagram  as  they  are  computed. 
Applying  Eq.  (31)  successively  to  joints  3,  4,  and  6,  one  obtains 


4 /•  i 

\J  + uuvj  yfi  » 


r*AAA  \ 

! 771.1. 


n t A / 1 AA\ 


1.2285 


7.6273  , 


P-.  + ±SL 

1.7024 


5.6359  ; 
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e«  = „ ♦ o = -8*4201  . 

56  3*5544  o.Viuj. 

The  unbalanced  exciting  moment  at  joint  5 is  computed  as 

M*  = 1.2285  (-1.7993)  + 2*3621(7.6273)  + 1.0522  (-8*4201)  + 3*5544(5*3103)  = 26*979 


Next,  consider  that 


01  = 0 and  e2  = 1.00  . 


The  rotations  © .are  . computed  in  the  same  manner  and  the  results  ars  recorded 
on  the  diagram  above  the  values  of  0*. 

The  total  slope  at  a joint,  say  at  joint  6,  is 

©6  = 5*635961  + 7.722962  . 

The  total  discontinuity  of  slope  and  the  total  exciting  moment  at 
joint  5 are 

e56  " e53  = “16.Q47©1  - 18.2 239 z , 

M5  = 26.979®1  + 33*931©2  • 

The  value  of  the  determinant  of  the  coefficients  of  6^  s-d  is 

1-16.047  -18.223 

26.979  33*931 


& s (&i  • &t)  ~ 


= -52.85 


Since  this  is  different  from  zero,  the  assumed  value  of  = 2.20  doeB  not 
correspond  to  a natural  frequency.  Repeating  this  procedure  for  several 
values  of  Xz  and  evaluating,  in  each  case,  the  resulting  determinant,  the 
curve  given  in  Fig*  33  was  obtained.  The  values  of  Xz  corresponding  to  the 
natural  frequencies  of  the  frame  are  recorded  on  the  figure. 

The  next  £.  -,ep  in  the  solution  is  the  determination  of  the  natural 
frequencies  for  which  and  ©2  are  simultaneously  equal  to  zero.  ©^  and  ©2 
can  be  equal  to  zero  only  when  bars  (1)  and  (2)  are  stationary  as  shown  in 
Fig.  34a. 
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Then,  M_-  = = 0 and  ^ c ^ = 0 

U&der  this  condition,  bars  (3)  and  (4)  can  vibrate  freely  only  at  frequencies 
equal  to  the  natural  frequencies  of  these  bars  for  the  condition  of  fixed 
ends.  For  bar  (3),  these  natural  frequencies  are  represented  by  values  of 

A,  = 4.730,  7.853,  10.996,  ....  (38a) 

The  equivalent  values  of  A;  are 

— 2.756  , 4.576,  6.406,  . . « '•  (38b) 

For  bar  (4),  natural  frequencies  corresponding  only  to  symmetrical  modes 
must  be  considered;  these  are  given  by  values  of 

A*  = 4.730,  10.996,  ....  (39a) 

which  ere  equivalent  to  - 

A2  = 3.535,  8.218,  ....  (39b) 

It  ia  now  necessary  to  ascertain  whether  or  not  these  frequencies 
are  natural  frequencies  of  the  portion  of  the  frame  composed  of  bars  (5), 

(6),  (7),  and  (8).  To  do  this,  it  is  necessary  to  carry  out  one  cycle  of 
the  basic  procedure  for  each  frequency  to  be  investigated.  For  the  pur- 
pose of  illustration,  three  different  frequencies  will  be  considered  in 
detail. 

(a)  A*=  2.756  (A , = 4*730).  Since  none  of  the  values  given  in 

Eq.  (39b)  is  equal  to  2.756,  bar  (4)  oust  be  stationary;  this  means  that 

M44!  c 0 • 

The  condition  of  equilibrium  at  Joints  3 and  4 requires  that 

**34  = ”^35  *®d  M46  = “^43  • 

Since,  for  the  A value  considered,  the  deflection  of  bar  (3)  is  symmetrical 
about-  its  mid-pas. 


-M43  = = -H« 
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The  joint  rotations  of  the  fraao  may  new  bo  expressed  in  tens 
of  Mjg,  which,  for  oonvenienoa,  may  bo  token  equal  to  1.0C  # starting 

from  joint  3 sad  oonoidoring  that  **  3,  one  may  deturni.no  the  rotationo  of 
joints  5»  6,  and  4 in  the  usual  manner*  If  A*  b 2*756  ia  a natural  frequency, 
the  computed  value  of  ©,  muat  bo  equal  to  aero  and  M.*  moot  be  equal  to 

V2 

-1.00  ~g-*  . 

The  data  neoeasary  in  oarrying  out  these  oomputatieae  are  asaembled 
in  the  following  table  and  in  Tigs*  34b  and  34o. 


Member 

2 

5 

7 

8 

6 

Xj/X, 

1.00 

0*7282 

1.687 

1.265 

0.7093 

h 

2*756 

2*01 

4.65 

3.49 

!e95 

In  fig*  34b  the  atiffheeaes  are  expressed  in  tens  of  the  *■  of  the  member 

* * T 

jm2*2 

to  which  they  refer,  while  in  Tig*  34o  they  are  expressed  in  terms  of  • 

**2 

The  rotation  ©j  is  computed  by  application  of  So*  (30)  to  joint  3,  while  ©^ 
and  ©^  are  oomputed  by  use  of  Sq«  (31)*  Sinoe  the  computed  value  of  ©^  ia 
different  from  aero,  A*  = 2*756  does  not  represent  s natural  frequency  of 

the  frame*  If  ©,  were  found  to  be  equal  to  aero,  it  would  have  been  nec- 

4 *2*2 
eaaary  to  investigate  also  if  vara  equal  to  -1*00  -£■*  * 

(b)  A*  = 4*576  ( Aj  * 7*853)*  la  this  oast  also,  bar  (4)  is  sta- 
tionary, but  ia  equal  to  instead  of  -Mjj#  In  all  other  reepeots, 
the  details  of  the  analysis  are  similar  to  those  of  the  previous  oaae*  The 
pertinent  computations  are  given  in  the  fblloylng  table  and  in  Hg«*  34d 
and  34a. 


Member 

2 

5 

7 

8 

6 

A • / X* 

1.00 

0*7282 

1.687 

1.265 

0.7093 

4.576 

3.33 

7.72 

5.79 

3.25 

I 
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Since  the  computed  rs lus  of  Is  different  from  aero,  A*  = 4*576  does  not 
correspond  to  * natural  frequency. 

(c)  A*  = 3.535  ( A*  x 4*730).  Since  none  of  the  value*  giver,  in 

Eq.  (38b)  ftseaaqual  to  3.535*  bar  (3)  in  Fig.  34a  stoat  be  stationary  in  this 
cess.  Tula  means  that  bars  (5)  and  (7)  are  stationary  also  and  that  joints 
4 and  6 remain  fixed  against  rotation.  * It  follows  that,  if  A*  = 3*535  is  a 
natural  frequency,  each  of  the  mashers  composing  penel  4-6-6  *-4'  must  be 
capable  of  vibrating  faeoly  with  the  ends  fixed.  For  bar  (6)  this  ic  poss- 
ible for  values  of 


\ - 4*730,  7.853,  .... 

(40a) 

while  for  bar  (8)  it 

is  possible  for  values  of 



\ s 4.730,  10.996,  .... 

(41a) 

The  corresponding  A* 

values  are,  respectively 

Aj  — 6.669,  11.071,  . . . . 

(40b) 

A*  — 3*739*  8.692,  .... 

(41b) 

Since  these  are  different  from  3*535*  A,  = 3*535  does  not  correspond  to  a 

A H 

natural  frequency.  It  is  worth  noting  also  that  tbs  A,  values  given  in  Eqs. 
(40b),  (41b),  and  (39b)  are  different.  It  is*  therefore,  concluded  that, 
within  the  range  of  the  frequencies  considered,  the  panel  formed  by  the 
members  (4),  (6),  and  (8)  cannot  vlbaN|te  freely  while  the  rest  of  the  frame 
remains  stationary. 

The  concluding  step  in  the  solution  of  the  problem  under  considera- 
tion consists  of  investigating  if  there  are  ary  natural  frequencies  for  which 
the  bars  seating  at  Joint  6 are  stationary- 

When  bars  (5)  and  (7)  are  stationary  ? 0 and  joints  3 snd 

6 remain  fixed  against  rotation  aa  shown  in  Fig.  35s.  Excluding  the  trivial 
case  of  no  vibration,  this  condition  eon  occur  only  at  frequencies  equal  to 
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the  natural  frequencies  of  bar  (8),  assuming  that  it  ia  flxad  at  the  ends, 
and  of  bar  (1),  assuming  that  joint  3 ia  fixed.  For  bar  (8)  the  firat  two 
of  theee  natural  frequencies  are  represented  by  the  A values  giren  in  Eq. 
(41)*  For  bar  (l)  the  first  two  natural  frequencies  are  giren  by  ralues  of 

X,  = 3.927,  7.069  (42a) 

which  are  equivalent  to  . ■ 

A*  = 4.225,  7,606  (42b) 

It  remains  now  to  investigate  whether  these  ralues  re  present  natural  fre- 
quencies of  the  frame  shown  in  fig.  35a.  Only  two  frequencies  will  be  con- 
sidered in  detail)  the  others  may  be  handled  in  a similar  manner. 

(a)  A t - 3.739  ( Aff  = 4.730).  Since  none  of  the  raluee  given  in 
Eq.  (42b)  are  equal  to  3.739,  bar  (1)  cannot  vibrate  freely  in  this  case) 
consequently,  both  bare  (1)  and  (3)  must  remain  still  and  joint  4 must 
remain  fixed  againat  rotation  aa  shown  In  Fig.  35b.  If  A*  = 3*739  represents 
a natural  f r»  luenoy,  each  member  of  the  frame  in  this  figure  must  be  capable 
of  vibrating  at  the  same  frequency.  It  has  already  been  shown  that  this 
condition  is  not  possible  for  panel  4-4-6 '-4 *•  It  remains,  therefore,  to 

i t 

ascertain  whether  panel  2-4-4  -2  can  vibrate  freely.  The  firat  two  natural 
frequencies  of  bars  (4)  and  (2)  are  represented  respectively  by  values  of 

A,  = 3.535,  8.218 
A,  = 3.927,  7.069 

Since  the  A{  values  for  the  two  bars  are  unequal,  it  ia  concluded  also  that 
panel  2-4-4' -2*  cannot  vibrate  freely  while  the  rest  of  the  frame  remains 
atill. 

(b)  A*  = 4.225  ( A,  = 3.927).  In  this  ease,  bar  (8)  cannot  vibrate 
freely)  therefore,  bars  (8)  and  (6)  in  Fig.  35a  are  stationery,  and  the  rota- 
tion of  joint  4 is  aero  as  shewn  In  Fig,  35c.  This  means  that  each  member  of 
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the  frame  shows  in  Fig*  35s  must  baht  re  as  if  it  vere  fixed  at  joints  3 and  4* 
The  A*  value a for  which  this  condition  can  be  Realised  are  different  for  the 
different  members}  consequently,  A * = 4*225  does  not  correspond  to  a natural 
frequency. 

In  summary,  it  should  be  stated  th-t  for  this  particular  problem 
and  for  the  range  of  frequencies  considered,  the  aero  Intercepts  of  the  curve 
in  Fig.  33  represent  the  complete  set  of  natural  frequencies  of  the  frame* 

22*  Comments  on  Hatnral  Modes  of  Partial  Vibration. 

The  determination  of  the  natural  frequencies  corresponding  to  modes 
of  partial  vibration  is  not  as  tedious  as  it  might  appear  from  the  quce  devoted 
to  its  discussion  in  the  preceding  sections*  Furthermore,  it  should  be  noted 
that,  for  most  practical  applications,  it  may  be  entirely  unnecessary  to  cal- 
culate these  natural  frequencies*  As  already  explained,  the  procedure  for 
determining  these  natural  frequencies  consists  of  (a)  computing  a number  of 
frequencies,  and  (b)  establishing  whether  or  not  these  frequencies  are  natural 
frequencies  of  the  frame*  In  most  actual  problems,  one  is  interested  in 
determining  natural  frequencies  comprised  within  specified  ranges  of  frequen- 
cies. Therefore,  if  the  values  calculated  in  the  first  step  of  this  procedure 
are  found  to  lie  outside  the  ranges  of  interest,  it  will  be  superfluous  to 
carry  out  the  second  step*  The  first  step  of  the  procedure  can  usually  be 
carried  out  almost  by  inspection. 

Natural  modes  of  partial  vibration  correspond  always  to  higher 
natural  frequencies.  Therefore,  no  consideration  need  be  given  to  these 
sodas,  if  only  the  Cutiuameii Lai  natural  frequency  of  a frame  is  to  oe  aeieiv 
mined.  The  fundamental  natural  frequency  of  a frame  may  be  determined  also 
by  the  moment  distribution  procedure  described  in  Section  13* 
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23 . Heed  Tor  Aroroadj**^  jfethgJg  of  Anelvsia. 

The  method  that  has  besn  described,  even  though  simple  bothoih 
principle  and  in  the  details  of  its  application,  may  beoome  time  consuming 
when  applied  to  structures  comprising  a very  large  number  of  members.  For 
very  complex  structures,  such  as  multi-story  multi-bay  building  frames,  it 
may  be  desirable  to  ha  be  a simpler,  though  less  accurate,  method  of  analysis. 

Strictly  speaking;  the  dynamic  response  of  a member  of  a framework 
and,  as  a consequence,  the  natural  frequencies  of  the  structure  ddpend  on  the 
properties  of  all  the  members  in  the  structure.  Intuition  leads  one  to  ex- 
pect, however,  that  the  importance  of  this  influence  diminishes  rapidly i ao 
the  distance  from  the  member  concerned  increases.  For  example,  the  natural 
frequencies  of  a two-span  beam  may  vary  greatly,  depending  on  whether  the 

extreme  ends  are  fixed  or  hinged;  cn  the  other  hand,  for  a multiple-span 

/ 

beam  the  natural  frequencies  may  be  almost  independent  of  the  condition  of 
restraint  at  the  extreme  ends.  For  example,  for  a uniform  beam  of  two 
equal  spans,  the  fundamental  natural  frequency  for  fixed  ends  is  56  percent 
higher  than  the  corresponding  natural  frequency  for  hinged  ends.  For  a 
uniform  beam,  of  seven  equal  spans,  the  fundamental  natural  frequency  for 
fixed  ends  is  only  6 percent  higher  than  for  hinged  ends.  (The  latter  value 
was  obtained  from  reference  (16)  )•  This  condition  indicates  that  It  may- 
be possible  to  determine  the  natural  frequencies  of  complicated  systems  by 
considering  only  a portion  of  the  structure  instead  of  the  entire  structure. 

It  is  believed  that  an  investigation  aimed  at  determining  the 
possibilities  of  such  an  approximate  procedure  will  be  moat  rewarding.  The 
procedure  described  herein  will  be  of  great  value  in  such  a study. 
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la  the  discussion  thus  far,  the  effect  of  permanent  axial  forces  on 
the  natural  frequencies  of  flexural  vibration  has  been  omitted.  This  effect, 
■which  may  be  quite  important  when  the  magnitude  of  the  external  force  on  the 
structure  is  a siseable  fraction  of  the  eritieel  buckling  load,  may  be  taken 
into  account  by  uae  of  modified  atiffneaa  end  cany-over  factors  which  in- 
clude the  influence  of  the  axial  thrust.  It  is  assumed  that  the  axial  loads 
are  independent  ef  time  and  that  they  are  applied  at  the  ends  of  the  members. 

The  modified  stiffnesses  and  carry-over  factors  may  be  expressed 
conveniently  in  tarsia  of  the  dimensionless  parameter  A , which  was  uaod 
previously,  and  the  ratio  P/P0,  where  £ is  the  axial  compressive  or  tensile 
foroe  and  is  the  fundamental  buckling  load  of  the  member  assuming  its 
ends  to  be  simply  supported.-  . 

Algebraic  expressions  for  the  dynamic  flexural  stiffness  and  tbs 
dynamic  flexural  carry-over  factor  are  given  in  Appendix  II.  Bumerieal 
values  of  these  two  quantities  and  of  their  product  have  been  computed  for  v 
values  of  A between  aero  and  4*75  at  increments  of  0.05  and  for  valuea  of 
P/P0  between  -4.0G  and  4.00  at  inerements  of  0.1.  It  is  expected  that  these 
results  will  be  made  available  soon. 

For  the  limiting  case  of  A *=  0,  the  values  of  stiffness  and  carry- 
over factors  are  those  for  an  axially  loaded  bar  which  does  not  vibrate.' 
Metalled  tabulations  ef  these  quantities  have  been  presented  previously  by 
James  (26),  by  Landquiet  and  Xroll  (27),  and  by  Eu  end  Libove  (28).  It 
becomes  apparent  that  the  problem  of  elastic  instability  is  • limiting  ease 
of  the  more  general  problem  of  the  vibration  of  structures  for  vbich  the 
members  are  subjected  to  axial  foreet. 

The  thrsa  slope  equatim  has  been  applied  to  the  investigation  of 
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the  Instability  of  ftraatsworks  by  Qbwslla  and  J old. ah  (2?)  aad  by  Vinter  and 
Ms  associates  (30).  The  procedure  described  in  this  report  baa  the  advantage 
that  it  can  be  need  to  obtain  solutions  with  considerably  leaa  effort* 
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V.  APPLICATION  OF  METHOD  TO  COHTHUOBS  BEAKS 
OH  FIXCIELB  SUPPORTS 


25.  General* 

Thus  far  application  of  the  method  haa  been  restricted  to  continuous 
beams  on  rigid  supports  and  to  continuous  frames  for  vhich  the  joint  do  not 
translate.  In  this  chapter,  the  nethed  is  extended  to  continuous  beans  on 
flexible  support*.  The  flexibility  of  the  supports  ia  represented  by  a set 
of  mutually  independent  deflectienal  and  rotational  springs.  To  start  with, 
it  ia. assumed  that  the  stiffnesses  of  the  restraints  are  independent  of  the 
vibration  frequency.  The  assumptions  made  in  the  analysis  are  similar  to 
those  made  in  the  case  of  continuous  beams  on  rigid  supports* 

The  continuous  beam  considered  is  shovn  in  Pig*  36*  As  in  the  case 
of  continuous  beams  on  rigid  supports,  consecutive  supports  are  numbered  from 
left  to  right,  starting  with  1 at  the  extreme  left  end  and  terminating  with 
z at  the  extreme  right  end* 


26.  Basic  Belations* 

— .ji. — - 

Figure  37  shovs  the  extreme  deflected  position  of  spans  J-l  and  j of 
a continuous  beam  undergoing  steady- state  forced  vibrations*  The  vibrations 
are  assumed  to  be  maintained  by  an  exciting  moment  applied  at  the  extreme 
right  end  of  the  beam. 

The  forces  acting  at  the  ends  of  each  span  and  the  deflections  and 
rotations  at  the  supports  ar«  indicated  in  their  positive  directions*  In 
addition  to  the  symbols  used  previously?  the  symbol  8 is  used  to  designate 

fa 

the  deflection  of  the  beam  at  a support;  the  stiffness  of  a defleetlonal 
spring  is  denoted  by  ]>,  while  that  of  a rotational  spring  is  denoted  by  £• 
Sinoe  the  beam  is  continuous  at  the  supports. 
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(Bj) l * ( Qj )*  * and  ( Sj)l  * (6j)p  - 5}  . (43) 

Also,  since  no  external  moment  or  force  acts  at  the  joint*  and  since  the  in- 

ternal moments  and  forces  must  be  in  equilibrium, 

Mj  - (Mjk  + (M.)*  + Rj6j  ' 0,  (44) 

Fj  - (Vj)L  + (Vj)K  * OjSj  - 0.  (45) 

xhs  ncaeata  and  shears  acting  at  the  ends  of  a span  may  be  ex- 
pressed in  terms  of  the  rotations  and  deflections  of  the  ends  of  the  span* 

For  example*, the  moment  or  shear  at  the  left  end  of  span  j m ay  be  obtained 
by  the  addition  of  the  following  four  component  effects! 

(1)  Moment  or  shear  produced  at  the  left  end*  when  that  end  is  rotated  by 
without  deflection  and  the  right  end  is  held  fixed* 

{2)  Moment  or  shear  produced  at  the  left  end,  when  that  end  is  held  fixed 
while  the  right  end  is  rotated  by  without  deflection* 

(3)  Moment  or  shear  produced  at  the  left  end*  when  that  end  is  displaced 
by  5 ^ without  rotation  and  right  end  is  held  fixed* 

(4)  Moment  or  shear  produced  at  the  left  end*  when  that  end  is  held  fixed 
and  the  other  end  la  deflected  by  5 without  rotation* 

The  direct  superposition  of  these  effects  is  justified  by  the  fact  that*  for 
a given  frequency  of  vibration*  the  moments  and  shears  are  linear  functions 
of  the  distortions*  Thus* 

aijh  - Kjdj  + OiK)j  6jn  +Qj5j  - (fQ)jSjH  , (46s) 

(Vj)R  - TjSj  - (tT)j  Sjh  * Qj0j  + (fQ)j  6j+,  . (47a) 

Similarly, 

(fljk  m fyidj  * (itK)j-i  6j.,  -Qj-tSj  + (yQ)j-,6j-t  , (46b) 

(Vjk  - Tj-iSj  - (tTJj-i  <5j-t  - Qj-,  Qj  - (^Qh-i  6j-i  . (47b) 


Subeti tuting  Sqs.  (46a)  and  (46b)  in  Eq*  (44)  and  Eqs« 
(471»)  In  Eq*  (45 )»  one  obtains 

rjj  m (qQ)j-i5j-i  + ( itK)j~i6j + CQj  - Qj-i)Sj 

f (Kj.,  tRj  * Kj)  6j  - (<fQ)j$j+,  + (‘kK)j  Ojti  m 0 , 

Fj  - ~(tT)j-,  Sy,  - (fQ)yi  dj-i  + (Tj-i  + + 7 })Sj 

+ (Qj  - 0j-;)6j  ~(tT)jfj+i  t (<jQ)jQj+,  * 0 . 
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(47a)  and 


(48a) 


(49) 


EHuu/aatlng  «j+1  from  Eq.  (49)  by  use  of  Eq*  (48a),  one  obtains 

[r/fQ h-t  * (tT)j-,]Sj.,  + [ifj(*K)j-i  *(fQ)j.i]Oj-i  * [ijSj -T}]6j 
- HiiFj-Sjlej  -hj(fQ)j-(tnj]Sj„  - o. 


in  which 


. (+Q)j 

(*K);  ' 


(50) 


Kj  - Kj-,  t Rj  * Kj  , 

- t im  Ti-i  * °j  + TJ  * 

Qj  m Qj  ~ Qf-i  • 

Equation  (50)  expresses  the  deflection  at  support  j+1  in  terms  of  the  de- 
flections and  rotations  at  the  two  preceding  supports*  The  rotation  at  sup- 
pert  J+l  may  be  obtained  in  terns  of  the  deflection  and  the  distortions 

at  supports  j and  j-1  from  Eq*  (48a) , which  may  be  rewritten  as 

-(*K)j.,ej.,  - QjSj  - KjSj  * (fQ)jSj„  - (*K)j  9jr,  - 0 . (MW 

Equation  (50)  assures  equilibrium  of  shears*  while  Sq*  (48)  assures  e qui- 
ll briun  of  moments  for  support  J.  Of  course*  both  equations  satisfy  the  con- 
ditions of  continuity*  If  support  £ is  rigid*  Eq.  (50)  is  satisfied  auto- 
matically* and  need  not  be  used.  If  * 5f  s ^l50’  that  is*  when 

three  consecutive  supports  ere  rigid*  then  Sq;  (48)  reduces  to  Eq*  (19a)* 
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7b*  foregoing  relation*  are  applieable  to  intermediate  -rapports  only* 
For  the  end  supports  the  following  epeeialised  relations  must  ha  need.  First, 
the  boundary  o onditions  for  the  extreme  left  support  are  considered*  Four 
different  eases  mat  be  distinguished! 

Case  I«  Both  &,  and  ls  are  eonsldored  finite  or  aero  in  this  ease*  Then,  the 
equilibrium  condition  fcr  moments  at  support  1 le 

ff,-  0,5,  + (R,  + K,)0,  - (<jQ),6t  + (*10,6:  - 0.  (51a) 

c ” 

The  corresponding  condition  for  shears  ia 

F,  - (Dl  + DS,  + Q,e,  - C t7},£  + (fQ),  & - 0.  (52a) 

Eliminating  #-  from  Eq*  (52),  one  ebtalna  an  equation  similar  to  Kq.  (50), 

[l,Q, - (£)  * T,)]Si  * faCRSKJ-Q,] 8,  - fafQjrUT),] & -0.  (53) 

From  this  equation,  i2  nay  be  determined  in  terms  of  3^  and  With  8 ^ 

deterBined,  ^ also  may  bo  eomputed  from  Eq*  (51a)  in  terms  of  ^ and 
Case  Dj  : , considered  infinite  (rigid  deflectienal  support)  and  finite 
or  aero*  Then,  the  equilihritai  condition  for  momenta  is  expressed  by  the 
equation 

M,  - (R,  + K,)e,  - (fQ),^  +(1<K)iez  - 0.  (51b) 

So  equilibrium  equation  for  shears  need  be  written,  since  at  a rigid  support 
this  ia  satisfied  automatically* 

Case  III.  B1  ia  considered  finite  or  sere,  but  Infinite  (clamped  end)*  In 
this  uaiv,  it  xe  ueoeseaxy  to  write  only  one  equilibrium  equation  for  shears! 
this  is 


F,  - (Q  + DS,  - (tT),Si  *(fQ)A  - 0. 


(52b) 


Casa  IT.  Bath  and  R,  sre  considered  infinite  la  thia  eaee  (fixed  end)* 

The  relation  between  the  moment  or  ahear  at  the  fixed  end  and  the  distortions 
at  the  second  support  are 


(*K),  % - (fQ),S,  , 

(54) 

(qQ),  dt  ~ (t  T),  6t  • 

(55) 

It  should  be  noted  that  Bio.  (53),  (51b),  (52b),' and-  (54)  or  (5$) 
involve  only  three  tudchoim  quantities* 

The  conditions  that  oust  be  satisfied  at  tha  right  end  of  the  beam 
are  as  follows t 


— •-  / • 

fl,  - 

(Kz.,  + R,)e.  + (•*/<),., - Q^J S,  * (fQl.6',  -0, 

(56) 

F,  - 

a,.,  *o,)5t  - (t  t),.,  - g,„  e,  - e,.,  - o. 

(57) 

For  Caro  H. 

~ ■ ■ - 

— — - 

6*  • 0 , 

(58) 

A, 

; - (KX.^RZ)QZ  + t (^Q)Z.,5Z.,  - 0 . 

(59) 

r°r  yap,  a. 

0,-0. 

(60)- 

■ “ (Tz.,  + DZ)SZ  - (tT)t.,Sz.,  - (fQ)z.,dz.,  - 0. 

(61) 

For  Case  17. 

* 

o 

■ 

(62) 

x - n 

(bj) 

27.  Outline  of  Procedure. 


The  rotation  and  deflection  of  the  bean  at  support  j may  be  written 


.(*•> 


Oj  - dju  + &jv  , 

ift  “ <5j  u * $ v , 

where  3 and  v are  dimensionless  par aae tore  which  represent  two  of  the  three 

unknowns  in  the  equations  expressing  tbs  boundary  conditions  for  the  left  end 

of  the  bean.  6*  and  <?j  are,  respectively,  the  rotation  and  deflection  at 

* .1 

support  j when  u = 1*00  and  v = 0r  and  9^  and  0 ^ are  the  corresponding  rota- 
tion and  deflection  when  u = 0 and  v = 1.00.  Since  the  natural  frequencies 
of  a system  depend  on  the  relative  values  of  the  deflection,  any  arbitrary 
amplitude  consistent  with  the  actual  boundary  condition  may  bo  chosen  for 
either  u or  v.  For  convenience,  the  following  values  ere  selected 1 


For  Case  lit 


II 

u a = 1.00 

and 

TC>* 

III 

U as  9^  sc  1*00 

sad  ~ 

•* « «2 

— r * 

nil 

u = = 1.00 

u = H~  * 1*00 

and 

t = #2 

17: 

and 

The  details  of  the  procedure  are: 

— 

1.  For  some  reference  span  of  the  beam,  say  span  £,  assume  a value 
of  A this  is  equivalent  to  assuming  a frequency  of  vibration 


u> 


■4/ 


frit 

mr 


2*  From  Eq.  (23)  compute  the  A values  for  the  remaining  spans  of 
the  beam* 


3*  From  Table  I in  Appendix  A and  -tits  A values  computed  in  the 
previous  step,  determine  the  stiffnesses  and  the  product  of 
the  stiffnesses  and  the  carry-over  factors  for  each  Bpan* 

4.  Identify  the  parameters  jg  and  v. 

5*  Consider  that  u = 1*00  and  v = 0*  Progressing  from  support 


A 


* 


to  support  across  the  beam,  determine  the  deflections  and  rota- 
tions at  all  supports.  Denote  these  by  3*  and  5*.  If  necessary, 
evaluate  also  the  unbalanced  moment  or  Bhear  at  the  extreme  right 
hand  support.  Denote  this  by  M*  or  The  distortions  of  the 
second  support  are  determined  from  the  appropriate  expressions 
given  in  Eqs.  (51)  through  (55 )•  The  distortions  of  the  remain- 
ing supports  are  evaluated  by  the  repeated  application  to  each 
support  of  Eqs.  (50)  and  (48b). 

6.  Repeat  the  preceding  step  by  considering  that  u - 0 and  v = 1.00. 

« fn 

Denote  the  resulting  distortions  by  6 and  o • If  necessary, 

determine  also  the  magnitude  of  the  unbalanced  moment  or  shear 

_n  _n 

at  the  extreme  right  hand  end.  Denote  thiB  by  M or  F • - 

0 8 

7.  The  actual  distortions  at  a support,  say  at  Bupport  j , are 


6j  = 6j  + Ojv  , 

Sj  = 5j  + 5j  y . 


(64b) 


Similarly  the  total  unbalanced  moment  or  shear  at  support  z is 


n.-  m + Kv , 
F.  - K*  Kv  ■ 


(65) 


8.  From  one  of  the  two  equations  expressing  the  boundary  conditions 
for  the  right  end  of  the  beam,  determine  the  unknown  parameter  v. 

9.  Evaluate  the  second  boundary  equation. 

10.  Repeat  Btepe  1 through  9 for  different  assumed  values  of  Ar  and 
plot  the  value  calculated  in  step  (9)  against  the  assumed  values 
of  Xr  . The  values  of  Ar  for  which  the  ordinates  of  the  result- 
ing curve  are  equal  to  zero  represent  the  natural  frequencies 
of  the  beam. 


74 


root  oi 


The  foregoing  procedure  can  be  modified  readily  to  include  concen- 
trated rigid  masses,  concentrated  sprung  masses,  intermediate  rigid  supports, 
flexible  supports  for  which  the  stiffness  depends  on  the  frequency  of  vib- 
ration, and  a continuous  elastic  aubgrade  of  tbs  Winkler  type.  For  con- 
venience, these  effects  are  discussed  separately. 

Rigid  Concentrated  Masses.  Assume  that  at  station  j the  beam  has  no  de- 
flectional  support,  instead  it  carries  a concentrated  mass  as  shown  in  Fig. 
38a.  Assume  that  the  beam  undergoes  a steady-state  forced  vibration  with  a 
frequency  • Then  for  a downward  deflection  8 y the  force  exerted  cn  the 
beam  by  the  mass  is  mjuf&j  , downward.  Had  the  beam  been  elastically  supported 
against  deflection  at  that  point,  the  corresponding  force  would  have  been 
DjSj  , upward.  Thus,  the  effect  of  the  concentrated  maea  is  equivalent  to 


that  of  a deflectional  spring  of  stiffness 

/ \ — * i4  Erlr  i 

(Qj)*,.-  t ^ £ 

If  the  mass  is  applied  at  a point  on  the  beam  that  is  supported  by  a de- 


(66a) 


flectional  spring  of  stiffness  as  shown  in  Fig.  38b,  the  stiffness  of 
the  equivalent  spring  is 


(Dj)c„  * Dj  - mjuf  . 


(66b) 


After  the  concentrated  mass  has  been  replaced  by  the  equivalent 
deflectional  spring,  the  natural  frequencies  of  the  system  may  be  determined 
in  the  usual  manner.  It  must  be  noted,  however,  that  the  stiffness  of  the 
equivalent  spring  is  a function  of  the  assumed  frequency  of  vibration  and 
must  be  evaluated  for  each  cycle  of  the  procedure. 


Concentrated  Sprung  Masses.  Assume  that  the  mass  m^  is  spring  borne  as  shown 
in  Fig.  38c.  Let  bo  the  stiffness  of  the  spring.  The  influence  of  the 


•.m; ‘ ^ 
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flaxible  support  may  be  taken  into  account  by  replacing  the  actual  mass  by 
an  equivalent  rigid  mass  of  magnitude 

^ 


/ - 


m jj  to* 

Sj 


(67) 


The  amplification  factor — -L, — j—is  determined  as  follows*  Let  Sj  be  the 
deflection  of  the  mass  while  Sj  is,  as  before*  the  deflection  of  the  elastic 
axis  of  the  beam.  The  forces  acting  on  the  mass  are  the  inertia  force  and 
the  spring  reaction j these  must  be  in  equilibrium;  therefore, 

/rijufSj  “ Sj  (6j  - Sj)  , 

whence,  th.9  amplification  factor  is 

JL  , _ 


I 


I - 


Sj 


It  should  be  noted  that,  for  vibration  frequencies  greater  than  the 
natural  frequency  of  the  spring  borne  mass,  the  amplification  factor  becomes 
negative,  and  Bq.  (67)  results  in  a negative  equivalent  mass. 

Rigid  Intermediate  Supports.  Consider  the  beam  shown  in  Fig * 3od  for  which 
support  4 is  rigid  instead  of  flexib.1*.  Let  the  outer  supports  be  flexible. 
The  natural  frequencies  of  this  beam  may  be  determined  by  a combination  of 
the  principles  that  have  been  described.  As  usual,  the  procedure  may  be 
initiated  at  support  1 by  taking  6-^  = 1.00.  For  the  portion  of  the  beam  be- 
tween supports  1 and  4 the  distortions  at  the  supports  may  be  expressed  as 

r 

the  sum  of  a constant  term  and  a term  involving  as  unknown.  The  magnl- 

•~r 

tude  of  this  unknown  may  be  determined  from  the  condition  that  4*  = 0.  For 
the  portion  of  the  beam  between  supports  4 and  7,  0^  may  be  selected  as  the 
now  unknown,  uuu  w*  .ivuo  a v £*uppOjri/s  mcy  De  c ompune d cy  x>he  re- 

peated application  of  Bqs.  (50)  and  (48b).  As  usual,  the  value  of  0^  may  be 
determined  from  one  of  the  two  equations  describing  the  boundary  condition 
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for  the  right  end  of  the  beam.  The  natural  frequencies  ere  those  fre- 
quencies for  which  the  second  boundary  condition  is  satisfied  identically* 

In  the  application  cf  this  procedure,  one  must,  in  general,  change 
temporary  unknowns  as  many  times  as  there  are  rigid  intermediate  supports* 
Abrupt  changes  in  the  magnitude  of  the  distributed  mass  or  in  the 
flexural  rigidity  within  a span  may  be  handled  by  assuming  that  the  beam  is 
supported  by  a flexible  support  of  zero  stiffness  at  the  point  of  the  dis- 
continuity. 

Supports  Having  Mass.  Consider  a construction  consisting  of  a cross  beam 
supported  cs  a series  of  longitudinal  girders  as  shown  in  Fig*  39a*  For 
simplicity's  sake,  the  connections  between  the  beam  and  the  girders  are  con- 
sidered hinged.  It  is  assumed  that  the  mass  of  the  supports  cannot  be 
neglected j consequently,  the  stiffnesses  of  the  supports  are  functions  of 
the  frequency  of  vibration. 

The  natural  frequencies  of  this  system  are  the  same  as  those  of  the 
crocs  beam  assuming  that  it  is  supported  on  deflectional  springe  as  shown  in 
Fig.  39b.  The  stiffness  of  aach  spring  is  equal  to  the  stiffness  of  the 
corresponding  girder.  The  stiffness  of  each  girder  is,  in  turn,  equal  to  the 
magnitude  of  a harmonically  varying  concentrated  force  which,  when  applied 
at  the  point  of  intersection  of  the  girder  and  the  beam,  will  deflect  the 
girder-  ty  a unit  amount.  The  magnitude  of  this  force  may  be  calculated  in 
terms  of  the  numerical  values  tabulated  in  Table  I,  Appendix  A*  As  an  ex- 
ample, consider  a girder  fixed  at  both  ends.  For  the  point  of  application 
of  the  force,  the  conditions  of  equilibrium  and  continuity  are  expressed  by 
the  equations 

(Qt -Q,)w„  + (K,*Kt)6.  - 0 , 

(T,-Ti)w„  + (GU-Q,)&,  - F,  . 
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where  the  subscripts  1 and  2 refer,  respectively,  to  the  portions  of  the 
girdsr  to  the  left  and  the  right  of  the  concentrated  force,  end  the  sub- 
script x refers  to  the  point  of  application  of  the  force.  From  the  simul- 


taneous solution  of  these  equations  one  obtains 

A - n t 

*v,  D " Tl  T*  K,  t K, 

For  a force  at  midspan,  this  equation  reduces  to 

D = 27T, 


(68  a) 


(68b) 


The  stiffness  of  a girder  with  different  boundary  conditions  may  be  de- 
termined in  a similar  manner. 

With  the  expression  for  the  effective  stiffness  of  the  deflectional 
spring  available,  the  natural  frequencies  of  the  system  may  be  determined  in 
the  usual  manner.  It  is  only  necessary  to  evaluate  the  stiffness  of  the 
springs  for  each  assumed  frequency  of  vibration  and  to  use  this  value  in 


the  calculations. 


•ade.  Consider  that  the  beam  is  supported  along  a 


portion  of  its  length  on  a continuous  elastic  subgrade  of  the  Winkler  type, 
as  shown  in  Fig.  39c « Let  d be  the  foundation  modulus,  defined  as  the  force 
per  unit  length  necessary  to  compress  the  foundation  by  unit  amount.  The  mass 
of  the  foundation  is  assumed  to  be  negligible. 

The  natural  frequencies  of  such  a system  may  be  determined  by  the 
same  procedure,  except  that  the  influence  of  the  elastic  foundation  must  be 
taken  into  account  in  evaluating  the  stiffnesses  and  the  carry-over  factors 
of  the  elastically  supported  spans. 

A vibrating  beam  which  In  not-  supported  by  sny  foundation,  is  acted 
upon  only  by  distributed  inertia  forces,  the  intensity  of  which  is 
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The  intensity  of  the  reactive  forces  produced  by  t':e  elastic  foundation  is 
« These  forces  act  opposite  to  the  inertia  forces  and,  in  effect,  reduce 
the  Intensity  of  the  latter  to 

- dw  . 

Nov-  the  ccoffloienus  of  dynamic  stiffnesses  and  the  carry-over  factors  depend 
solely  on  the  dimensionless  parameter  A*  f_  . therefore,  the  effect  of 

the  foundation  can  be  taken  into  account  by  replacing  in  this  expression  the 
quantity  mu?  by  the  reduced  value  of  mu?-d  a 

If  muf-d  is  positive,  A is  real,  and  the  various  stiffness  and 
carry-over  factors  may  be  obtained  directly  from  Table  I in  Appendix  A.  If, 
however,  mu?-d  is  negative,  A is  imaginary  and  the  values  of  Table  I can 
no  longer  be  used.  In  this  case,  the  expressions  for  dynamic  stiffness  and 
dymanic  carry-over  factors  may  be  obtained  from  the  expressions  for  static 
stiffness  and  static  carry-over  factors  of  a bar  on  elastic  foundation.  It 
is  only  necessary  to  replace  in  the  latter  expressions  the  quantity  d by  the 
reduced  value  d -mu? , These  expressions  may  be  obtained  readily  from  sol- 
utions given  by  Hetenyi  (31) • 

Consider  now" the  case  fbr  tfhich  the  foundation  extends  along- the 
efitire  length  of  the  beam  as  shown  in  Fig,  39d.  The  modulus  of  the  foundation 
may  differ  from  one  span  to  the  other.  Let  d^  be  the  foundation  modulus  for 
span  j.  In  what  follows,  it  is  shown  that,  if  the  natural  frequencies  of 
the  system  without  the  subgrade  are  known,  under  certain  conditions,  the 
natural  frequencies  of  the  system  with  the  subgrade  may  be  computed  directly. 

Assume  that  the  beam  without  subgrade  and  the  beam  with  the  sub- 
grade are  in  a steady-state  of  vibration.  Let  co  be  the  circular  vibrat5.on 
frequency  of  the  beam  without  the  subgrade  and  u>5  be  the  frequency  of  the  beam 
with  the  subgrade.  Then,  for  span  r the  A value  of  the  beam  without  tl 
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eubgrade  is 


5 - Lr 

V£r-  Lr  J 

and  tihO  value  for  "the  bean  with  the  subgrede  is 

(*')*  '=\/ fft (*U ~ %)  L" 

For  any  other  span  the  values  for  the  two  cases  are 


lj  ■ 

^ -JW^W  g 

“d  “ «t  = \ = cen^nt , (70) 

Ws  - .rmj  E777  Lr  Jj 

(K  h V mr  e}  Lj  Lr  2r  • 

Now,  if  and  Qr)3  are  equal,  the  numerical  calculations  involved  in  carry- 
ing out  a cycle  of  the  procedure  will  be  Identical  for  the  two  systems s 
Therefore,  if  3^  corresponds  to  a natural  frequency  of  the  beam  without  the 
subgrade,  then  G)^  will  correspond  to  a natural  frequency  of  the  beam  with 
the  subgrade o Equating  Eqs.  (12)  and  (69),  one  obtains  the  following  rela- 
tionship between  the  natural  frequencies  of  the  two  systems. 

(***  )n  = + (Jj/ftij)  - (71) 

The  subscript  N designates  natural  frequencies.  Equation  (71)  can  be  applied 
to  beams  having  any  number  of  spans  and  any  boundary  conditions,  provided  that 
the  relationship  given  in  Ea.  (70)  is  satisfied. 


After  the  distortions  of  the  beam  at  the  aids  of  a span  have  been 


ib  lk'4  j < 1L»  - 
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evaluated,  the  deflection  configuration  of  the  span  may  be  determined  by 
superimposing  the  relieving  four  component  configurations! 

(1)  Deflection  configuration  produced  by  the  rotation  of  the  left-  end 
of  the  span  when  that  end  is  restrained  against  deflection  and  the 
right  end  is  fixed. 

(2)  Deflection  configuration  produced  by  the  rotation  of  the  right  end 
of  the  span  when  that  end  is  restrained  against  deflection  and  the 
left  end  is  fixed. 

(3)  Deflection  configuration  produced  by  the  deflection  of  the  left  end 
of  the  span  when  that  end  is  restrained  against  rotation  and  the 
right  end  is  fixed,  end 

(4)  Deflection  configuration  produced  by  the  deflection  of  the  right 
end  of  the  span  when  that  end  i?  restrained  against  rotation  and 
the  left  end  is  fixed. 

Effects  (1)  and  (2)  can  be  obtained  readily  by  multiplying  the  end 

b 

rotations  by  the  appropriate  values  given  in  Table  II,  paying  proper  atten- 
tion to  signs.  Effects  (3)  and  (4)  may  be  calculated  from  Eq.  (B— 34)  in 
Appendix  B.  These  effects  may  also  be  expressed  in  terms  of  the  quantities 
given  in  Table  I as  follows*  Consider  a clamped  ended  uniform  beam,  the  left 
end  of  which  is  displaced  by  unit  amount*  The  deflection  amplitude  at  a point, 
a distance  x from  the  end  being  displaced,  may  be  obtained  by  writing  for  the 
point  the  two  equilibrium  equations  (see  Eqs.  48a  and  49); 

a * (<2z-G,)>v,  i-(K,+Kt)e.  + (<jQ),  - 0, 

F,  - a,  +t2k  + (Qt-a,)e.  - an  ■■  o . 

Th9  subscripts  1 and  2 denote,  respectively,  the  left  and  the  right  portions 
of  the  beam.  Eliminating  from  the  first  of  the  foregoing  two  equations  the 
unknown  8 , and  solving  for  the  deflection  *v  , one  obtains 

X X 


flfii-waS*********  r 


'try.  + 

/T  r > (QrQ,? 

(T,  + Tz)  tc.  +*- 


(72a) 


The  deflection  at  the  same  point  of  the  beam  due  to  a unit  deflection  ct  the 
right  end  may  be  obtained  from  Eq.  (72a)  by  interchanging  the  subscripts  1 


and  2.  When  x = jj,  this  equation  reduces  to 

CtT), 

" IT, 


(72b) 


Note  that  for  1=0  this  equation  yields  ■£,  which  is*  of  course,  the  correct 
value  for  the  static  deflection  at  midspans 

30«  Illustrative  Example* 

Ebrwrnple  7.  In  order  to  illustrate  the  details  of  the  procedure  and 
present  a tabular  scheme  for  arranging  the  computations,  we  shall  calculate 
the  first  three  natural  frequencies  of  the  hypothetical  continuous  beam  shown 
in  Fig.  40.  The  dimensions  of  this  beam  and  the  stiffnesses  of  the  rotation- 
al and  aeflectional  restraints  ere  shown  in  Fig.  40.  The  various  quantities 
are  expressed  in  terms  of  the  physical  properties  of  a reference  span  r.  In 
this  particular  case,  we  choose  r = 4» 

It  is  convenient  to  express  the  stiffnesses  and  the  carry-over 
effects  for  each  span  in  terms  of  the  El  and  L of  the  reference  span  r.  To 
do  this,  the  values  obtained  from  Table  I in  Appendix  A must  be  mu  ltiplied  by 
certain  dimeusionless  factors  as  follows:  The  coefficients  of  K and  kK  for 

a span  ) raue:t  be  multiplied  by  the  factor  (Eq.  24).  The  coefficients 
of  Q and  aQ  must  b«  multiplied  by 

E.Xj  L 2 

h’tlfz  ’ (73> 

trlJ 

and  the  coefficients  of  T and  tT  hy 
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These  relatione  aay  be  verified  readily* 

For  the  beam  considered,  the  ratio  A;/Arand  the  factors  0(^5 
and  ifj  are  evaluated  in  Table  21.  This  table  is  independent  of  the  fre- 
quency of  vibration  and  is  computed  but  once* 

Table  2B  presents  one  complete  cycle  of  the  procedure  for  a trial 
value  of  A*=  \=  1*30*  This  value,  shown  encircled  in  the  r-th  Cl  Column 
(1),  corresponds  to  a circular  frequency  of  vibration  U)  _ is|2  /fl  # Columns 

(1)  through  (18)  in  this  table  are  conveniently  filled  in  the  following  orders 

(2) ,  (9),  (1),  (3  and  8),  (4  and  7),  (10  and  12),  (5),  (6),  (11),  (13),  04), 

(15),  (16),  (17),  and  (18).  It  should  be  noted  that  all  quantities  in  Col- 
umns (2)  through  (12)  are  expressed  in  terms  of  the  gj  and  L of  the  reference 
span  r.  The  value  of  -5*7122  in  Column  (9),  which  represents  the  stiffness 
of  an  equivalent  deflectional  spring  having  the  same  effect  as  the  concen- 
trated mass  m^,  is  computed  from  Eq*  (66a),  The  partial  distortions  6«, 

5n,  and  9!  are  evaluated  in 'Columns  (19) .through  (22) it  The  parameter  r 
is  determined  from  tha -condition  that 


-j*-  = -87.5781  + 3436*36  v = 0 ; 


whence , 


v = 0.0254357  . 


The  total  distortions  at  supports  4 and  5 are  given  in  Columns  (23)  and  (24). 
The  magnitude  of  the  exciting  moment 

K5  = [4.473(2.327)  + 2.021  (-1*177)  + 6.089(0*3438)]  ~ = 10.12  ~ . 


Under  tho  action  of  this  moment  the  distortions  at  the  extreme  left 


support  are 


= 1.00  and  - 0.02548571. 
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Had  8,/Lf  instead  of  61t  been  taken  equal  to  unity,  the  magnitude  of  "Ghe  ex- 

•*» 

citing  moment  would  have  been 

M - H - o«r-  l S 

5 “ 0.0254857  L “ l 

Since  is  different  frcm  zero,  the  assumed  value  of  \ - 1*30 
does  not  correspond  to  a natural  frequency.  By  repeating  several  such  cycles 
of  computation  for  different  values  of  \ , the  curves  shown  in  Fig.  41  were 
obtained.  The  sclid  curve  in  this  figure  shows  the  moment  necessary  to 
produce  a rotation  8-,  = 1.00,  while  the  dashed  curve  shows  the  moment 
necessary  to  produce  a deflection  § ^ = 1,00 L The  A intercepts  of  these  . 
curves  correspond  to  natural  frequencies  of  vibration.  The  first  three  crit- 
ical values  are 

( A4)x  = 0.85  , ( A4)2  = 2.01  , ( A4)3  *=  2.55  . 

It  should  be  noted  that  the  distortions  5 and  8 are  obtained  as 
small  differences  between  large  quantities.  It  becomes  necessary,  therefore, 
to  use  In  the  computations  a relatively  large  number  of  significal  figures, 
particularly  if  the  higher  natural  frequencies  of  the  beam  are  to  be  evalu- 
ated. It  1b  recommended  that  at  least  6 or  7 significant  figures  be  retained 
throughout  the  computations. 
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71.  APPLICATION  OP  METHOD  TO  FRAMES 
WITH  SIDESWAI 

Hi  SggetgJ.oa.la  Single-Bay « Multl-Stonr  Frames, 

The  method  described  in  the  preceding  Chapter  can  also  be  used  to 
determine  the  natural  frequencies  of  lateral  vibration  of  frames  for  which 
the  joints  are  free  to  translate.  Though  this  method  can  be  applied  to  frames 
of  any  complexity,  it  can  efficiently  be  used  only  for  symmetrical,  b ingle- 
bay,  multi-story  frames.  For  such  frames,  the  details  of  application  of  the 
procedure  are  similar  to  those  for  e continuous  beam  on  flexible  aupports. 

Consider  the  symmetrical  frame  shown  in  Fig.  42a.  Because  of  symmetry, 
the  end  rotations  of  each  girder  are  algebraically  equal,  and  the  midpoints  of 
the  girders  are  points  of  inflection.  Consequently,  it  is  possible  to  consider 
in  the  analysis  only  one  half  of  the  structure  as  shown  in  Fig.  42b.  In  this 
figure,  the  right  ends  of  the  girders  can  rotate  and  slide  freely  in  the  hori- 
zontal direction  but  can  not  deflect  in  the  vertical  direction.  The  influence 
of  these  horizontal  members  may  be  represented  by  a set  of  concentrated  masses 
and  concentrated  rotational  springe,  attached  at  the  points  of  intersection  of 
the  horizontal  members  and  the  main  column  as  shown  in  Fig.  42o.  The  magnitude 
of  each  concentrated  mass  in  Fig.  42c  is  equal  to  the  total  distributed  mass  of 
the  corresponding  girder  in  Fig.  42b.  Similarly,  the  stiffness  of  each  rotav 
tional  spring  in  Fig.  42c  is  equal  to  the  flexural  stiffness  of  the  corres- 
ponding girder  in  Fig.  42b. 

The  natural  frequencies  of  thia  system  may  be  determined  by  the 
procedure  described  in  Seotiona  27  and  28.  It  must  be  borne  in  mind,  how- 
ever, that  the  stiffnesses  of  the  rotational  restraints  are  sot  constant  in 
this  case,  but  rather  depend  on  the  frequency  of  vibration,  and  nnust  be  eval- 
uated for  each  cycle  of  the  procedure. 
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VII.  EXTENSION  OF  METHOD  TO  CONTINUOIE  PLATES 

32.  General. 

This  chapter  is  concerned  with  the  determination  of  the  natural 
frequencies  of  flexural  vibration  of  rectangular  plates  which  are  simply 
supported  along  two  opposite  edges  and  which,  in  one  direction,  are  con- 
tinuous over  rigid  supports  transverse  to  the  simply  supported  edges.  The 
plate  may  have  any  number  of  panels  of  arbitrary  length.  The  mass  per  unit 
of  area  and  the  flexural  rigidity  of  the  plate  may  vary  from  one  panel  to 
the  other,  but,  in  any  one  panel,  these  quantities  are  assumed  to  remain 
constant. 

The  method  of  analysis  is  similar  to  that  described  in  Chapter  III 
for  the  case  of  continuous  beams  on  rigid  supports.  The  assumptions  made  in 
the  analysis  are  those  embodied  in  the  ordinary  flexure  theory  of  medium  thick 
plates  composed  of  an  elastic,  homogeneous,  and  isotropic  material.  In  add- 
ition, it  la  assumed  that  the  supports  can  offer  no  torsional  restraint  and 
that  no  frictional  forces  or  horizontal  sheering  forces  may  develop  between 
the  plate  and  the  supports.  Throughout  this  discussion.  the  effects  of  damp- 
ing, rotatory  Inertia,  and  shearing  deformation  are  disregarded. 

/ . 

33.  Basie  of  the  Method. 

Figure  43  shows  the  type  of  the  continvous  plate  considered.  The 
coordinate  axes  are  taken  parallel  to  the  edges  as  shown  in  Fig.  43*  The  sides 
parallel  to  the  y-axis  are  assumed  to  be  simply  supported.  411  transverse 
supports  4ire  considered  to  be  rigid.  Along  the  extreme  edges  the  plate  may  be 
hinged,  fixed  or  subjected  to  a rotational  elastic  restraint  of  constant  stiff- 
ness. I 'or  simplicity,  it  will  first  be  assumed  i.hat  the  extreme  right  edge 
is  either  hinged  or  elastically  restrained  against  rotation.  A fixed  edge 
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will  be  treated  separately  at  the  end* 

For  the  plate  shove  in  Fig.  43*  the  solution  of  the  differential 
equation  for  transverse  vibration  of  plates  reveals  that;  during  free  vib- 
ration, the  deflections,  rotations,  shears,  and  bending  momenta  along  sections 
perpendicular  to  the  simply  supported  edges  are  proportional  to 

Sin-2~-  C05  u)Ht  . (75) 

In  this  expression, represents  a circular  natural  frequency,  t time,  and 
a the  width  of  the  plate;  n is  an  integer  which  designates  the  number  of  half- 
sine waves,  alternately  upward  and  downward,  in  the  distribution  of  the  de- 
flection, slope,  shear,  or  bending  moment  along  the  plate  width*  Corresponding 
to  each  value  of  n there  is  an  infinite  number  of  natural  frequencies* 

Consider  now  that  the  plate  undergoes  a Bteady-state  f orced  vib- 
ration, such  that  the  rotation  of  the  plate  at  the  extreme  left  support  is 

6,  (x.t)  “ 0,  sin  -----  cos  u>t  , 

a 

where,  cv  ia  the  circular  frequency  of  vibration,  3 ia  the  same  integer  used 
in  Eq.  (75),  and  is  the  maximum  amplitude  of  elope  at  the  left  support* 

The  magnitude  of  is  assumed  to  have  a prescribed  finite  value.  If  support 

1 ia  fixed  instead  of  hinged  or  elastically  restrained,  6^  ia  equal  to  zero, 
and  it  becomes  necessary  to  assume  that  the  bending  moment  at  that  support  is 

M,  (x,t)  - n sin -n~~-  cos  wt  . 

ia  assumed  to  have  a fixed  value.  In  either  case,  the  vibration  of  the 
plate  13  presumed  to  be  maintained  by  an  exciting  couple  applied  »t  the  ex- 
treme right  hand  support  z.  Then,  the  solution  for  steady-state  forced  vib- 
ration of  the  governing  differential  equation  reveals  (a)  that  the  deflections, 
elopes,  shears*  and  bending  momenta  along  sections  perpendicular  to  the  simply 
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supported  edges  are  proportional  to,  and  in  phase  with,  the  rotation  (or 
bending  moment)  at  the  extreme  left  edge,  and  (b)  that  the  exciting  moment 
may  be  expressed  as 

nr(x,t ) - Mt  sin  ~~~  cosuit  . 

Obviously,  the  magnitude  of  K_  depends  upon  the  frequency  of  vib~ 

ration.  For  a frequency  equal  to  a natural  frequency  of  the  plate  ( U)  =U/^)» 

M becomes  equal  to  zero,  but  the  deflection  of  the  plate  remains  finite.  The 
z 

deflected  surface  of  the  plate,  which  represents  a natural  mode  of  vibration, 
is  then  expressed  as 

w(x,y.t)-  Y„  sin  -22*.  cos  tvt  . 

a 

where,  Yn  is  a function  of  the  y-coordinate  onlyj  its  absolute  value  depends 
on  the  value  assigned  to  the  amplitude  of  slope  (or  bending  moment)  at 
support  1. 

In  the  discussion  thus  far,  it  was  assumed  that  the  extreme  right 
support  of  the  plate  was  either  hinged  or  elastically  restrained.  If  support 
z is  fixed,  the  criterion  for  a natural  frequency  ia  that  = 0, 

Details  of  the  Method. 

The  foregoing  considerations  suggest  the  following  procedure  for 
calculating  the  natural  vibration  frequencies  of  continuous  plates  having  two 
opposite  edges  simply  supported, 

1.  Assume  that  the  amplitude  of  slope  or  bending  moment  at  the  extreme 
left  edge  of  the  plate  is  distributed  Binusoidally,  and  that  it  has 
a fixed  value.  Since  the  naturel  frequencies  of  a system  depend  only 
on  the  relative  values  of  the  deflection,  any  arbitrary  amplitude. 
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consistent  with  the  actual  boundary  condition,  may  be  chosen.  For  • 


-TW-}  m*T7 1 T jnVT*! 


* 

»- 


! - 

I 

* 


♦ 


90 

hinged  edge  or  for  a partially  fixed  edge,  take 

9^(x)  = 1.00  sin  22S 

For  a fixed  end,  8^  = 0;  therefore,  take 

M^x)  = 1.00  sin  3-2 

In  the  application  of  thia  procedure,  it  ia  necessary  to  consider 
a specific  value  of  £ in  each  case.  Let  n = nQ« 

2,  Choose  a frequency  of  vibration  and  determine  the  rotations  of  the 
plate  over  the  supports,  and  from  these  determine  the  magnitude  of 
une  exciting  moment  at  the  extreme  right  edge.  (For  a fixed  edge 
the  moment  at  that  edge  need  not  be  evaluated).  These  quantities 
are  determined  by  identically  the  same  procedure  that  was  used  for 
continuous  beams.  The  pertinent  relations  are  reviewed  in  sub- 
sequent paragraphs. 

3*  Repeat  the  previous  step  for  different  frequencies  and,  from  the 
results  obtained,  plot  a curve  of  the  Variation  of  the  exciting 
moment  or  of  the  slope  at  the  right  edge  as  a function  of  the  fre- 
quency of  vibration.  If  the  right  edge  is  hinged  or  elastically 
restrained,  plot  the  variation  of  the  exciting  moment;  then,  the 
natural  frequencies  are  thoae  frequencies  for  which  the  exciting 
moment  vanishes.  If  the  right  end  is  fixed,  establish  the  variation 

of  the  rotation  € ; the  frequencies  for  which  8.  vanishes  are 

» s 

natural  frequencies* 

4.  For  the  natural  frequencies  determined  in  step  (3) - the  deflection 
of  the  plate  along  sections  perpendicular  to  the  hinged  edges  has 
n^  half-sine  waves.  In  general,  steps  (1)  to  (3)  must  be  repeated 
for  as  many  values  of  n as  may  be  necessary  for  a given  problem. 
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For  a given  frequency  of  vibration,  the  rotations  of  the  plate  over 
the  supports  and  the  magnitude  of  the  exciting  moment  may  be  determined  from 
the  same  equations  that  were  used  to  anelyze  continuous  beams  on  rigid  sup- 
ports, namely  from  Eqs.  (19),  (21),  and  (22).  Equation  (19)  expresses  the 
condition  of  equilibrium  and  continuity  at  an  interior  support  of  a continuous 

beam,  by  relating  the  rotations  of  the  beam  over  three  consecutive  supports. 

♦ 

Equations  (21)  and  (22)  express,  respectively,  the  boundary  conditions  for 
the  extreme  left  and  the  extreme  right  ends  of  the  beam. 

In  the  application  of  theae  equations  to  continuous  plates,  it  is 
only  necessary  to  interpret  as  the  maximum  amplitude  of  elope  along  the  3-th 
support,  and  as  the  maximum  amplitude  of  moment  at  that  support.  In  add-*- 
ition,  it  is  necessary  to  replace  the  stiffness  and  carry-over  factors  for 
beams  by  equivalent  quantities  for  plates.  Such  quantities  are  defined  in 
the  aection  that  follova. 


35.  Hastio  Constants  for  a Panel  of  a Plate. 

Consider  a rectangular  plate  simply  supported  on  three  edges  and 
fixed  on  the  fourth  as  shown  in  Fig.  44.  Let  edge  f be  subjected  to  a steady- 
state  rotation  without  dsf lection,  such  that  the  magnitude  of  the  rotation  is 
given  by  the  relation 

■ • Of  (x,ty-  sin  cos  ujt 

The  amplitude  of  the  moment  at  edge  f necessary  to  produce  this 
rotation  is  proportional  to  sizr-jp  and  may  be  written  as 

Mf  = K0f  . (76) 

The  quantity  S is  a measure  of  the  resistance  to  steady-state  forced  rota- 
tion of  edge  f of  the  plate,  when  edge  g is  fixed,  and  is  referred  to  as 
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the  n dynamic  flexural  stiffness*  of  edge  f. 

Tbs  amplitude  of  the  moment  induced  at  the  fixed  edge  Is  also 
proportional  to  sinSK  , and  may  be  exp1'0  seed  as 

fl 


M - WCG*  = kJ* 
g it 


(77) 


The  quantity  k,  which  represents  the  ratio  of  the  moment  at  the  far  fixed 
edge  to  that  at  the  edge  being  rotated,  is  defined  as  the  "dynamic  flexural 
carry-over  factor"* 

These  quantities  are  strictly  analogous  to  those  defined  for  beams, 
and  they  are  used  in  the  analysis  of  plates  in  identically  the  same  way  as 
the  quantities  for  beams  are  used  in  the  analysis  of  continuous  beams* 

The  quantity  K may  be  expressed  as 

(78) 


where  H is  the  flexural  rigidity  of  the  plate  per  'unit  width  and  b is  its 
span  length  in  a direction  parallel  to  the  simply  supported  edges*  The 
coefficient  Cg  is  dimensionless  and  depends  on  the  ratio  and  the  di- 
mensionless parameter 


/ Pku)1 

A “ V / N 


(79) 


where  p denotes  the  density  and  h the  thickness  of  the  plate*  The  carry- 
over factor  k is  dimensionless  and  depends  pa  the.  ne  tva  quantities  as 
the  coefficient  Cg* 

The  pertinent  analytical  expressions  for  K and  k are  given  in 
Appendix  B.  For  = 0,  that  is,  when  the  plate  does  not  vibrate,  the 
values  of  K and  k reduce  to  those  given  by  Newraark  (15)* 

Assume  now  that  the  plate  is  not  vibrating  but  is  instead  sub- 
jected along  its  simply  supported  edges  to  uniformly  distributed  compress- 
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ire  force*  px.  Let 


*'■ 


b1  Pa 
7T*  N 


(80) 


Thea¥  it  can  be  shown  that  the  values  of  & and  k of  a vibrating  plate  for 
given  v aluee  of  and  A are  numerically  equal  to  the  stiffness  and  carry- 
over factor  of  an  equivalent  compressed  plate  having  the  same  value  but 
a value  of 

„ vO-  ») 

Numerical  values  of  K and  k for  compressed  plates  have  been  tab- 
ulated by  U.  D.  Kroll  (32)*  as  a function  of  and  k>  • It  should  be  noted 
that  Kroll  defines  stiffness  as  the  moment  necessary  to  produce  a maximum 
rotation  amplitude  of  1/4  radian  instead  of  one  radian.  Consequently,  the 
stiffness  values  obtained  from  Kroll  *8  tables  mist  be  multiplied  by  4 to 
make  them  conform  to  the  definition  given  in  this  report. 

As  a simple  illustration,  we  shall  determine  the  dynamic  flexural 
stiffness  and  dynamic  flexural  carry-over  factor  for  a rectangular  plate 
having  a ratio  of  sides  a/b  = 0.5  when  n = 1 and  A*  = 4*00.  The  equivalent 
value  of  k*  , 

k«  = (165(0.25)  « 4.00  . 

From  page  14  of  reference  (32)  we  find 


I = 


N 


~ = 2.54385  J and  k = 0.126499  , 


whence 


H 


K = 10.1754  £ and  k = 0.126499  . 


See  Section  2c  of  Appendix  B. 

+ Kroll  uses  the  symbols  A for  , k for  k*  , £ for  K,  and  C for  k. 
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36.  Hluitratlvo  Example* 

Example  8;  This  example  Involves  a rectangular  plate  Bimply 
supported  along  two  opposite  edges  and  continuous  over  five  equally  spaced 
rigid  supports  as  shown  by  the  sketch  in  Fig*  45*  The  plate  is  simply 
supported  at  the  extreme  left  edge  and  fixed  at  the  extreme  right  edge. 

All  panels  are  square.  It  is  desired  to  calculate  the  first  few  natural 
frequencies  of  this  structure. 

The  plate  is  analyzed  in  identically  the  same  manner  as  the  con- 
tinuous beam  considered  in  Example  1.  In  fact  the  expra salon  for  9*.  is  the 


same  as  that  for  the  beam ; 
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r4  - 8k2  + 8 


The  carry-over  factor  k,  which  depends  on  the  parameters  end  A , was 
determined  from  the  numerical  values  reported  in  reference  (32)  as  described 
in  the  preceding  section.  In  Fig.  45  9^  has  been  plotted  as  a function  of 
(A*)z  (’nb'f  for  values  of  n = 1 and  n = 2.  The  values  of  (A*)* (7^) 
corresponding  to  natural  frequencies  are  recorded  on  the  figure. 

The  lowest  circular  natural  frequency  for  n - 1 is 


D2  Ph 


The  lowest  frequency UJH for  n = 2 is 


For  purposes  of  comparison,  the  corresponding  circular  natural  frequencies 
of  o square  plate  haying  (a)  all  four  edges  simply  supported  and  (b)  three 
edges  simply  supported  and  one  fixed  are  given  in  the  following.  For  edge 
condition  (a): 


1 


when  n s i 


vhsn  n = 2 

Pftr  wtm  J4  _ /«  % 

r -—gw  Cvuui  UUQ  (0,»S 

vhen  n = 1 
when  n = 2 


Theae  valuea  hare  been  obtained  from  reference  (33) 
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Till  SUMMABI 


37.  General. 

A mathod  hra  been  or«**nted  in  this  report  for  the  determination 
of  the  undamped  natural  frequencies  and  of  the  corresponding  natural  nodes 
of  flexural  vibration  of  elastic  structures.  The  method  has  been  applied 
to  continuous  beams  on  rigid  or  flexible  supports » to  continuous  frames 
without  side sway,  to  symmetrical  single- bay  multi-story  frames  for  which 
the  joints  are  free  to  translate , and  to  continuous  plates  having  two  op- 
posite edges  simply  supported. 

The  method  is  a generalization  of  Holser'a  method  for  calculating 
the  natural  frequencies  of  toraional  vibration  of  shafts  end,  like  Holzer * s 
method,  it  has  been  reduced  to  a routine  scheme  of  computation  which,  when 
repeated  a sufficient  number  of  times,  will  give  the  natural  frequencies  of 
the  system  to  any  desired  degree  of  accuracy.  The  method  la  baaed  on  the 
fact  that  the  exciting  couple  neceaaary  to  maintain  a dynamical  system  in  a 
ateady-atate  of  forced  vibration  with  finite  amplitudes  becomes  equal  to 
zero  at  any  one  of  the  natural  frequencies  of  the  system- 

Extensive  tables  of  numerical  values  for  the  various  quantities 
entering  in  the  analysis  are  presented  in  Appendix  A of  this  report.  With 
these  tables  the  calculations  required  in  the  application  of  the  method  to 
particular  problems  are  simplified  greatly.  The  tabulated  values  may  be 
used  alao  with  other  analytical  techniques  as  well  as  for  the  analysis  of 
the  ateady-state  forced  vibration  of  structures. 
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APPENDIX  B.  DERIVATION  OF  FORMULAS 


Formulae  for  Barg  Without  Axial  Forces 
a.  General  Solution  of  Fundamental  Equation  for  Vibration  of  Bare 


Using  Bernoulli-Euler's  beam  theory,  the  differential  equation 
for  transverse  vibration  of  s uniform  bar  which  le  free  from  external 
loads  is 


El 


dV 

dx* 


+ m 


y_ w 


- o , 


(Bl) 


where  w(x,t)  is  the  deflection  of  any  point  of  the  bar  at  a time  t, 
m is  the  mass  per  unit  of  length  of  the  bar,  and  El  is  the  flexural 
rigidity  of  its  cross  section.  The  effects  of  damping,  shearing  de- 
formation, and  rotatory  inertia  are  disregarded. 

We  can  let 


w(x.t)  * Yfx)cosuut  , (B2) 

where  the  amplitude  of  the  motion  l(x)  is  a function  of  x only,  and 
ui  is  the  circular  frequency  of  the  motion. 

Substituting  (B2)  in  (Bl).  , one  obtains  the  following  expression 
for  determining  the  function  X(x): 

Y"  - (-£-)*  Y * 0 , (B3) 

where  L is  the  span  length  of  the  bar  and 

>-yWL-  «> 

Primes  in  Eq.  (B3)  indicate  differentiations  with  respect  to  x. 

The  solution  of  Eq.  (B3)  1b 

x x k v 

Y (x)  - c,  cosh  A — + c,  sinh  A -j-  -t  cscosAj^  + C^sin  A -j-  . (B5} 
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The  integration  ronatsnt*  e^.  r_^;  22Ci  auot  be  determined  from 
the  boundary  conditions  of  the  particular  problem  considered. 

b.  Formulas  for  Elastic  Constanta 

Consider  the  bar  shown  in  Fig.  B-l.  At  x = 0 the  bar  ia  fixed , 
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•a 
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I * 


n-*Kd  ( 


,t= 


V-^QQ' 

rM  r»  ^ 

rxg.  d-a 

while  at  x = L it  is  subjected  to  a moment 


^ )n-Kd 

-V-Q6 


f! 

- >. 

K(x,t)  = Mcoswt 

(B6) 

-• 

producing  at  that  end  a steady-state  forced  rotation 

_ 

■Jt 

- yc 

©(x,t)  = OcosuA  . 

(B7) 

y. 
• “ 

With  the  notation  and  sign  convention  given  in  Section  2, 

the 

ITS 

1 

boundary  conditions  may  be  stated  as  follows: 
At  x = 0 J = Y = 0, 

1 
- jj 

•t 

e = i»=  o , 

and  at  x = L 5 = 1 = 0 > 

(B8) 

4 

•■is 

& 

m = Eirn  . 

(B9) 

% 

From  these  boundary  conditions,  one  finds  the  following  values  for 
the  integration  constants  in  Eq.  (B5): 

Hi1  sinhX-sinX 

l cosh 

(BIO) 


C,  - -Cj 


Ctm  ~C*  - 


2 A*  El  coshX  sin>.  - sinhX  cosX 

Ml 2 


cczhX-ccs  A 


ZX El  coshX  sinX  - sinh  A cosX 
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The 

moments  M and 

the  reactions  V at  the  ends  of  the  bar  are 

found 

-4 

to  be  as 

follows  ? 

■ 

At  x = L 

M = K©. 

(Bll) 

where 

K « 

cIY"  l . cosh X sinX  - sinh  Xc-osX  El 

(B12) 

- 

Y'  1 - cosh  X cosX  L 

> 

and 

-v  = oe , 

(B13) 

where 

Q 

EIY'"  l2  sinh  A sinX  El 

(B14) 

4 

Y"  'v  / - cosh  A cosA  L* 

■4- 

At  x = 0 

m » kie , 

(B15) 

where 

k Y*]*.o  _ Sinh  A~  sin  A 

(B16) 

- 

Y'Jx.l  cosh X sinX -sinh XcosX 

5 

and 

V = qQ6, 

(217) 

where 

Y"']x.o  cosh  X ~ cos X 

Y"']k.l  sinhX  sinX 

(BIS) 

•> 

iS 

For  a simply  supported  beam,  the  moment  at  x = h may  be  expressed  as 

M = K"©  . (B19) 

It  can  be  shown  that 

2 sink  A sinX 


K"=  A 


(B20) 


cosh  A sin  A - sinh  A cos  A 
Equations  (BA2),  (B16),  and  (B20)  have  bsen  presented  previously  by 
Gaskell  (13) - 

Now  consider  the  case  In  which  the  right  end  of  the  beam  is  dis- 
placed without  rotation  by 

S (r.t)  = $ ccsut  , (£21 ) 

as  shown  in  Fig*  B2. 


M 
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rc-*as  (f 
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V-  tTS 


4p 

j ' -M  = Q5 

]/’  78 


Fig.  32 

In  this  case,  the  boundary  ccnditiona  may  be  stated  as  follows: 

At  x = 0 6=1=0. 

6 = Y'  =0. 

At  x = L 6 = Y»  = 0 , 

5 = Y . 

For  these  conditions,  one  finds  the  following  values  for  the  integration 
constants  in  Eq.  (B5 ) - 


(B22) 


(B23) 


c, --C,  - o 


******  I*  1 — me  a 


3 ~ 2 / - cosh  A cos  A 


cz  * -c4 


jL.  sinhX  + 5in\ 
2 1-coshXcosX 


(B24) 


The  moments  M and  the  reactions  V at  the  ends  of  the  bar  are  found 
to  be  as  follows: 

At  x = L M = ‘-Q<5.  (B25) 

where  (J  is  given  by  Eq.  (B14),  and 


where 


T s- 


V = T<5, 

E1Y'"  cash X sin X + sinhX  cos X 


8 I ~ coshX  casA 

At  x = 0 M = -qQ5, 

where  q is  given  by  Eq.  (B18),  and 


(B26) 

(B27) 

(B28) 
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where 


V = -t?5, 

Y'"]x.o  sinh  X + sin  X 


(B29) 

(B30) 


Y "j,mL  cosh  A smX  + sinh  X cos X 

As  already  remarked,  the  quantities  £ in  Eqa.  (B13)  and  (B25)  are  equal; 
likewise  the  quantities  q in  Eqs.  (B17)  and  (B28)  are  equal.  These 
equalities  follow  from  Lord  Rayleigh’s  reciprocal  relations  which,  for 
convenience,  are  reviewed  in  Section  Id  of  this  Appendix* 

It  should  be  stated  that  numerical  values  of  the  trigonometric  and 
hyperbolic  functions  appearing  in  the  numerator  and  denominator  of  most 
of  the  formulas  given  in  this  section  have  been  tabulated  previously  in 
Reference  (24).  These  functions  were  tabulated  for  values  of  A be- 
tween zero  and  10.00  at  increments  of  0.02.  The  numerical  values  pre- 
sented in  this  report  have  been  tabulated  at  increments  of  A of  0.01# 

Since  all  quantities  were  computed  on  the  Electronic  Digital  Computer , 
there  was  no  need  to  use  the  previously  tabulated  functions. 

c . Formulas  for  Deflections  of  £ Bar  Due  to  Distortions  at  the  Ends 

It  is  desired  to  find  the  deflections  of  the  bars  shown  in  Figs.  B1 
and  B2.  The  distortions  9 and  5 may  be  introduced  at  either  end  of  the 
tar. 

First  consider  the  bar  shown  in  Fig.  Bl.  Let  x denote  the  distance 
of  any  point  of  the  bar  from  the  end  being  rotated.  Then,  if  9 represents 
the  rotation  at  the  left  end,  the  deflection  amplitude 


I*  = C9L  , 


(B31a) 


and  If  9 represents  the  rotation  at  the  right  end 


I,  = -C9L  , 


(B31b) 


■'JJ'rW'f  tr*  mr  * •-  ♦ ■ i 
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where  C is  a dimensionless  coefficient  dependent  on  the  coordinate  z and 
the  parameter  A • 

Substituting  the  integration  constants  given  in  Eq.  (BIO)  into  Eq» 
(B5),  the  following  expression  is  found  for  C; 

C'  mF^hT^T)  {&>hX-sinX)[ash(l-$)X  -*»(!-&*] 


- [cosh  A - cosX] [sinh  Cl- 8) A - sin  Cl ~3)X]J  , 


where 


S =? 


(B32) 


(B33) 


Consider  now  the  bar  shown  in  Fig.  £2.  Let  x denote  the  distance  of  a 
point  of  the  bar  from  the  deflected  end.  Then,  for  a deflection  ampli- 


tude 8 at  either  end 


1-  = C*<5, 


(B34) 


where  C'  is  a dimensionless  coefficient. 


Substituting  the  integration  constants  given  in  Eq.  (B24)  into  Eq. 
(B5),  the  following  expression  is  found  for  C's 


d‘  TUI-'mh AcosA)  /MA-cosAjfcoshO-JH -obs<I-$)X] 

1 j (B35) 

-[sinhX*  sjnX][sinh(IS)A  ~ sin(l-§)\]j  . 


Lis  Recinroca] 


Rayleigh's  reciprocal  relations  for  dynamics  (14)  are  strictly  ana- 
logous to  Maxwell* e law  of  reciprocal  relations  for  statics.  Rayleigh's 
relations  may  be  stated  as  follows:  If  A and  B denote  two  points  in  a 

given  structure,  the  steady-state  forced  displacement  at  A produced  by 
a harmonically  varying  load  applied  at  B is  equal  to  the  steady-state 
displacement  at  B due  to  the  same  load  applied  at  The  term  dis- 
placement may  be  interpreted  in  a general  sense  as  either  linear  or 
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A 


* 


angular  displacement;  similarly,  the  term  load  may  be  interpreted  in  the 
same  general  sense  as  either  force  or  couple.  Couples  correspond  to  ro- 
tations while  forces  correspond  to  linear  displacements. 

Using  Rayleigh’s  reciprocal  relations,  one  can  prove  readily  that 
MBller-Breslau’s  principle  is  also  applicable  in  the  case  of  steady- 
state  forced  vibrations.  The  proof  is  identical  to  that  for  the  static 
case  and  is,  therefore,  omitted.  According  to  this  principle*  The 
influence  line  for  SDy  dynamic  function,  such  as  reaction,  shear,  bend- 
ing moment,  torque,  at  some  point  A of  * structure  due  to  harmonically 
varying  load  can  be  obtained  as  the  deflected  shape  of  the  structure  due 
to  a very  small  unit  displacement,  linear  or  angular,  introduced  at 
point  £. 

It  follows  from  this  principle  that  the  deflection  of  a fixed  ended 
beam  subjected  to  a unit  rotation  at  one  end  represents  an  influence 
line  for  fixed  end  moment  due  to  a unit  concentrated  force  on  the  beam. 
Similarly,  the  deflection  of  a beam  resulting  from  a unit  deflection 
without  rotation  of  one  end  represents  an  influence  line  for  fixed  end 
shear  due  to  a concentrated  unit  force. 


2,  Formulas  for  Bars  with  Axial  Forces 

a.  General  Solution  of  Governing  Differential  Equation 

The  differential  equation  for  transverse  deflection  w(x,t)  of  a uni- 
form bar  which  is  free  from  lateral  loads  but  is  acted  upon  by  constant 
axial  forcea  P is 


- , d fy 

CLJ?--r 


_ 7)\ 


dx: 


m 


A*  u/ 

dtl 


= 0. 


A positive  £ indicates  a tensile  force. 


(B36) 
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The  deflection  w(x,t)  may  be  expressed  as 

w(x,t)  = Y(x) cos  wt  . 


<Bj57  ) 


where  Y(x)  is  the  amplitude  of  the  harmonic  motion  and  u>  is  its  cir- 
cular frequency. 

Substituting  Sq.  (B37)  into  Eq.  (B36)  and  using  the  symbols 


p.=  - 

t r2EI 
L*  ’ 

(B38) 

A-4/ 

v' 

f muj*  . 

El  L ■ 

(B4) 

one  obtains  the  following  expression  for  determining  the  function  Y(x)s 

r'TT  r' {~bJ r'°  (B39) 

The  roots  of  the  corresponding  characteristic  equation  are 

f <t> 

n.z  * L - 

x 


(E40) 


r c i 


where 


’•-itm***'* 

x- 


A* 


(B41) 

(B42) 


Then,  the  solution  of  Bq.  (B39)  becomes 


Y(x)  = C,C05ht‘j~  +CzSinh  f-j-  +C)Cos\-j-+CtSmX-j-  ■ (B43) 


X 


The  integration  constants  c^»  c2»  and  c^  must  be  determined  from  the 
boundary  conditions  of  the  particular  problem  considered. 


b.  Formulas  for  Flexural  Stiffness  and  Flexural  Carry-Over  Factor 

Consider  a bar,  such  as  that  shown  in  Fig.  Bl,  fixed  at  x = 0 and 


-i 
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'i 


irl 

■ jj 

• 4. 
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subjected  to  a periodic  moment 

M(x,t)  = Mcosut 

at  x = L.  Assume  that  the  bar  is  acted  upon  by  tensile  end  forces* 
The  boundary  conditions  in  this  case  are  the  same  as  those  given 
in  Eqs.  (B8)  and  (B9).  For  these  boundary  conditions,  one  finds  the 
following  values  for  the  integration  constants  in  Eq.  (B43)* 


_ ML*  / Xsinh9~  9$m  X 

ci“  c*  El  9*+Xl  fcoshf  smX  - XsinhfcosX 

_ X Hi 1 X cosh  9 - cos  X 

c*~  El  9l+  X*  fcosh9sm\-  \smhfcosX 


(B44) 


The  moments  at  the  ends  of  the  bar  are  found  to  be  as  follows  t 


At  x = L 

M = KG  = KY«, 

(B45) 

where  K “ 

fcosh<PsinX~  XsinhfcosX  El 

(B46) 

Q*f\x[l~cosh9cosX]+  swh4>$inX 

At  x = 0 

M = kKe  , 

y»  ] sinh  9 — v~ 5,n  X 

4,  D - ’ Itro  - 2: 

y "1  a 

1 J*~L  cosh  9 sin  X ~ sinh  <P  cos  \ 

(B47) 

where 

(B48) 

Now,  assume  that  the  bar  is  subjected  to  compressive  forces.  In 
this  case,  it  is  necessary  to  replace  +P  by  -P.  Then,  Eq*  (BA1)  is 
changed  to  Eq*  (B42)  and  vice  versa.  Therefore,  the  expressions  of  K 
and  k for  a compressive  force  can  be  obtained  from  the  corresponding 
expressions  for  a tensile  force  simply  by  interchanging  in  the  latter 


oipieBeione  the  quantities  T and  « 


For  the  special  case  in  which  no  axial  force  is  present,  9“  \m  A , 
and  Eqs.  (B46)  and  (B48)  reduce  respectively  to  Eqs.  (B12)  and  (B16). 
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3f  formulas 


£t  Gsao^al  Solution  of  Fundamental  Equation  for  Vibration  of  Plateg 

Using  the  ordinary  flexure  theory  of  medium  thick  plates,  the  diff- 
erential equation  for  lateral  deflection  v(x,y,t)  of  a uniform  plate 


which  is  free  from  lateral  loads  is 


y»v  , - 34>v  d*>v  , rh  _ n 

ax4  a/dy-  ay*  /v  at** 


(B49) 


where  = the  density  of  the  plate  material,  h = the  thickness  of  the 
plate,  assumed  to  be  constant,  and  N = the  flexural  rigidity  of  the 


plate. 

The  coordinate  axes  are  taken  in  the  middle  plane  of  the  plate 
parallel  to  the  sides  of  the  plate.  The  origin  of  the  axes  is  taken  at 
the  upper  left  hand  corner  of  the  plate. 

If  the  plate  is  simply  supported  along  two  opposite  edges  (say, 
along  x = 0 and  x = a)  the  solution  of  Eq.  (B49)  may  be  taken  in  the 
form 


w(x,y,t)  = Yn  co5u)t , 


(B50) 


where  is  a function  of  the  y variable  only.  Substituting  this  equa- 
tion into  Eq.  (B49),  one  obtains  the  following  equation  for  determining 


the  function  I : 
n 


•XT'"  >?  (-^  - --jf1)  Y„-0  (Bji) 


The  roots  of  the  corresponding  characteristic  equation  ere 


* 


(B52) 


where 


[Thu? 
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(B53) 
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Two  different  cases  must  now  be  considered i 
Case  (l)  ; vh'jr  )*=*  (^fj2  t and 

Case  (II).  when  % < I^J2  . 

The  particular  case  in  which  A - 0 is  omitted. 

In  the  first  case  two  of  the  roots  are  real  while  the  remaining  two  are 
imaginary.  In  the  second  case  all  four  roots  are  real. 

First  consider  Case  (l).  Using  the  notation 


~9h 

1! 

(B54) 

and 

irl 

(B55) 

nJb_Y_  -x* 
a J A 


the  solution  of  Eq.  (B51)  becomes 

Yn(y)  * c,  cosh  $ Jjj-  + ct  sinh  f -y  + c3  cos  V^-  + c+sm \ (B56) 

In  Case  (II)  the  \ values  In  Eq.  (B55)  are  imaginary  and,  using 
the  notation 

**  rfW~h  * ■ (W7) 

the  solution  of  Eq*  (B51)  becomes 

X, (y)  * c\  cosh  <P  + c*  sinh + d,cosh  \ + c4sinh  X^j~  ■ (B58) 

The  integration  constants  In  Eqs.  (B56)  and  (B58)  must  be  determined 
from  the  boundary.  Qonditlqus  of.  thp  particular  problem  considered. 

b.  Formulas  for  Flexural  Stiffness  and  Flexural  Carry-Over  Factor 

Consider  a rectangular  plate  simply  supported  at  x = 0 and  x = a 
and  fixed  ct  y = C.  Let  an  ejtuitiiig  moment 

M(x,t)  = Main-'—  coswt  (B59) 

be  applied  along  the  edge  y = b 


The  boundary  conditions  in  this  case  a ret 


at  y = 0 


5 = \ • 0 ■ 


0*o) 


O-I'  = 0 . 
n 


it.  v ?=  b 


<5  = 1=0, 

XX 


(B61) 


M = HI” 
n 


The  moment  amplitudes  at  the  ends  nay  he  expressed  in  ter: 
rotation  amplitude  9 as  follows t 

At  y = 0 M = K©  , (B62) 

At  y = b M = kKO  . (363) 

where  K and  k are  respectively  the  flexural  stiffness  and  the  flexural 
carry-over  factor. 

It  should  be  noted  that  when  A*  Eq.  (B51)  for  plates  is  ‘ 

similar  to  Eq.  (B39)  for  bars  subjected  to  fixed  axial  forces.  And 
since,  for  the  particular  case  considered,  the  boundary  conditions  given 
in  Eqs.  (360)  and  (B61)  are  similar  to  those  given  in  Bqs.  (B8)  and  (B9)» 
the  expressions  of  K and  k for  a plate  panel  can  be  obtained  directly 
from  the  corresponding  expressions  given  in  Eqs.  (346)  and  (348).  It 
is  only  necessary  to  replace  in  Eqs.  (B46)  and  (B48)  the  quantities 
9 . X , El,  and  L by  9 , X , N,  and  b,  respectively.  The  results  are 
as  follows  j 


^7 \\  [ I"  cosh ? cos \]  + sinh  ¥sm  \ 


9 cosh  $ sin  X - X sinh  9 cos.X 

V 


N 


(B64) 


9 


sinh  9 — sin  \ 


'k  m —g — 

cosh  9 sin  X - sinh  $ cos  X 


(B65) 
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These  relatione  apply  to  rallies  of  Af  >'(^)^» 

For  A*  < the  expressions  for  K and  k are  obtained  from  those 
given  in  Eqs.  (B64)  and  (B65)  by  replacing  trigonometric  functions  by 
hyperbolic  functions  and  the  quantity  (A)1  by  “CR')1  . The  results  aret 


y _ ^ cosh  ^ sinh  X - \ sinh?  coshX N_  (&&) 

[l-co,h?cosh?]t3p^,sinhfsinhX  b 


$ - _ 
sinh  f - -*pr  X' 

* - r 7 . (067) 

Y cosh  ? sinh  X'  - sinh  f cosh  V 

c.  Correlation  Between  KLastlc  Constants  for  Vibrating  Plates  and 

GoftHreaagfl  nates 

Consider  a flat  plate  loaded  on  tvo  edges  parallel  to  the  y-sxls  by 
uniformly  distributed  compressive  forces  p^.  Assume  that  no  lateral 
load  acts  on  the  plate.  _ 

The  differential  equation  for  the  deflection  w(x,y)  of  the  plate  is 


(B68) 


For  the  case  in  which  the  loaded  edges  are  simply  supported;  the 


solution  of  Bq.  (B68)  may  be  taken  in  the  form 

v(xfy)  = t 


(B69) 


vhsrsTn  is  a function  of  y only.  The  unloaded  edges  may  have  any  con- 
dition of  restraint. 

Substituting  Bq.  (B69)  into  Bq.  (368);  one  obtains  the  following 

*»***•■>  4-4  AVt  *4  Ie4  e»  V • 
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(H7C) 
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where  4?  * — tt (t~  ■ (ET71) 

7T  N 

Comparing  Eq.  (B70)  with  Eq*  (B51)  it  can  readily  be  shown  that  if 

• (B72) 

the  two  equations  will  be  identical.  Accordingly*  the  quantities  K and 
k of  a vibrating  plate  for  a given  value  of  and  a*  are  numerically 
equal  to  the  corresponding  quantities  of  a compressed  plate  for  the  same 
•jg"  value  but  for  a value  of  k'm  nb ) * Since  tabulated  values  of 
stiffness  and  carry-over  factor  for  compressed  plates  do  exist  (see 
Reference  32),  these  values  can  be  used  to  determine  the  corresponding 
dynamic  quantities,  ft  should  be  noted,  however,  that  in  Reference  32 
stiffness  has  been  defined  as  the  moment  necessary  to  produce  a rota- 
tion of  a maximum  amplitude  of  a quarter  radian  rather  than  one  radian. 
Therefore,  the  stiffness  values  obtained  from  Reference  32  must  be  multi- 
plied by  4 to  make  them  conform  to  the  definition  given  in  this  report. 
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APPENDS  C.  ANALI3IS  OF  STEADI-STATE  FORCED  VIBRATIONS 

ix  General  Description  of  Method  of.  Analygle 

The  information  presented  in  this  report  may  be  used  also  to  ana- 
lyse the  steady-state  forced  vibration  of  continuous  beams  and  frames  sub- 
jected to  harmonically  varying  forces.  The  forces  are  presumed  to  have  the 
same  frequency  and  bhe  same  phase  angle.  The  effect  of  damping  is  neglected. 

The  determination  of  the  steady-state  forced  vibration  of  a system 
is  a much  simpler  problem  than  the  calculation  of  its  natural  frequencies, 
since  in  the  former  case  it  is  only  necessary  to  go  through  a single  cycle 
of  the  trial-and-error  procedure  used  in  the  determination  of  natural  fre- 
quencies. 

An  analysis  for  steady-state  forced  vibrations  can  be  carried  out  in 
substantially  the  same  way  as  an  analysis  for  static  conditions.  In  either 
case,  the  analysis  Involves  two  basic  steps: 

a.  Determination  of  the  redundant  quantities  at  the  joints  or 
supports  of  the  continuous  system.  In  general,  the  .redundant  quantities 
may  be  momenta,  rotations,  deflections  or  any  combination  of  these. 

b.  Determination  of  the  moments,  shears,  deflections  at  points 
between  joints  or  supports*  The  principal  difference  between  a static  and 
dynamic  analysis  is  that  while  in  a static  analydis  momenta  and  shears  with- 
in a member  may  be  determined  by  statics  from  the  moments  and  shears  at  the 
ends  of  the  member,  in  a dynamic  analysis  these  effects  must  be  computed 
independently. 

We  shall  now  outline  one  of  a number  of  possible  procedures  for 
executing  these  steps.  For  illustration,  we  shall  consider  the  relatively 
simple  frame  shown  in  Fig.  B3  and,  for  simplicity’s  sake,  Bhall  assume  that 


Zrs.'t*  ■ 
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the  joints  of  the  frame  do  not  move:.  Joints  i and  9 are  assumed  to  be 
simply  supported. 


i F,cosu)t 

4 I Q 


\F?cosu)t 


Fig.  B3 

Although  this  procedure,  in  the  form  presented  herein,  is  appli- 
cable to  continuous  open  frames  only,  it  can  readily  be  extended  to  frames 
involving  closed  panels*  The  details  of  the  procedure  are  as  follows: 

Step  a*  (1)  Replace  the  concentrated  forces  by  equivalent  fixed 
end  moments.  The  magnitude  of  these  moments  may  be  obtained  directly  from 
Tabic  II  in  Appendix  A.  (If  the  ends  of  the  loaded  members  were  free  to 
deflect,  it  would  have  been  necessary  to  calculate  also  fixed  end  shears. 
These  could  have  been  computed  either  from  Eq.-,  (72)  or  from  Eq.  (B3l). 

(2)  For  each  span  compute  the  quantities  K and  k£.  (If  the  joints 
of  the  frame  were  free  to  translate,  it  would  have  been  necessary  to  compute 
also  the  quantities  £>  £2’  I.  and  tT. ) These  quantities  are  obtained  direct- 
ly from  Table  I in  Appendix  A. 

(3)  Assume  a slope  at  support  1.  Progressing  from  joint  to 

joint  across  the  frame  in  the  manner  described  in  the  body  of  this  report, 

evaluate  the  rotations  of  the  joints.  It  should  be  remembered  that  in  this 

case  the  external  moment  M,  is,  in  general,  different  from  zero. 

«! 

(4)  From  the  rotations  determined  in  the  preceding  suep,  calculate 
the  magnitude  of  the  moment  at  the  extreme  right  end  of  the  frame.  This 
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mossnt  will,  in  general,  be  different  from  the  aoifeal  noaenx  at  tbnt  and,  a 
condition  indicating  that  the  assessed  slope  d-,  was  in  error.  It  ie  therefore 
necessary  to  add  a correction  configuration. 

(5)  Disregarding  ths  sstsrs&l  moments , find  the  influence  of  a 
unit  rotation  at  joint  1 on  the  distortions  of  the  regaining  joints.  Cal- 
culate also  the  magnitude  of  the  corresponding  moment  at  the  right  sad. 

(6)  Combine  the  configurations  determined  in  (4)  and  (5)  so  as  to 
eliminate  the  unbalanced  moment  at  the  right  end. 

(7)  From  the  rotations  determined  in  (6),  compute  the  aagnitude 
of  the  moments  at  the  joints  or  supports  of  the  frame  using  the  relationship 

M0j=W0  + k »j 

where  MoJ  ie  the  moment  at  end  o of  the  member  oj . This  concludes  Step  a* 

Step  b.  Any  desired  effect  (moment,  shear,  deflection,  ete.)  with- 
in a member  oj  may  be  determined  by  superimposing  tbs  following  effeotai 

(1)  The  effect  of  the  load  on  a simply  supported  member  * 

(2)  The  effect  of  the  moment  at  end  o 

(3)  The  effect  of  the  moment  at  end  j 

m 

If  the  ends  of  the  member  deflect,  the  following  effects  must  also  be  added s 

(4)  The  effect  of  the  deflection  at  end  o 

(5)  The  effect  of  the  deflection  at  end  £ 

Saoh  effect  is  computed  on  the  assumption  that  the  member  la  simply  supported 
at  the  ends. 

The  quantities  added  may,  if  desired,  be  different  from  those  out* 
lined.  For  example,  one  may  add  the  effect  of  the  load  on  s member  fixed  at 
both  ends,  the  effects  of  the  end  deflections,  and  the  effeote  of  the  tod 
rotations.  Irrespective  of  the  manner  in  which  the  computations  are  carried 


out 9 the  facility  with  which  results  are  computed  depends  on  the  information 
available  for  the  various  quantities  added* 

The  appropriate  expressions  for  effects  (2)  through  (5)  are  given 
in  Reference  (24).  Included  in  this  reference  are  also  tabulated  numerical 
values  of  functions  in  terms  of  which  these  effects  can  easily  be  computed. 
The  appropriate  expression  for  effect  (1)  is  given  in  the  next  section* 


2.  Formulas  for  Steady-State  Deflection  of  a Simply  Supported  Uniform  Beam 
Carrying  a Concentrated  Force  Fcoswt 

The  force  is  assumed  to  be  applied  at  a distance  'i  L from  the  left 

end  of  the  beam  as  shown  in  Fig.  B4»  The  subscripts  1 and  2 will  be  vised  to 

i F cos  tot 


•/S 


( I ) 


(2) 

U~t)L 


X,' 


Fig.  B4 


designate  quantities  for  the  left  and  the  right  portions  of  the  beam;  The 
distanced  and  are  measured  from  the  ends  of  the  beam  as  shown  on 
the  figure. 

The  deflections  w^  and  may  be  written  as 
w,  * [c,  cosh  X -jr-  + czsinh  X -jp  + c3  cos  X + c4  sin  X -jjJ  cos  tot  , 

iv2  * jj:,  cosh  X-j*-  + cz  sinh  X ^ + C3cosX  -j-*  + c+smX  ^ j cos  wt  . 

The  eight  integration  constants  have  been  evaluated  from  the  four  boundary 
conditions  end  from  the  four  continuity  and  equilibrium  conditions  for  the 
point  of  application  of  the  force.  The  results  ares 
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Pros?  these  equations,  expressions  for  bending  Moment  end  shear  can 
be  obtained  readily  by  differentiation. 


'.f 

t. 

j 

: .13 


Vi»**  '»<'">• 


'»  A.  ■ V.V-. 


!'  - : M'  < h<\V  ‘jft  Atf,  tlNk/-u  'fc^iir  i-iriji; 


BIBLIQGFAPHI 


(1)  "Sui  Emiilluxig  uer  Grvmdschwingtmgszaklen  von  durchlaufenaen  Tr&gern", 

by  W.  Mudrak,  Ingenieur-Archiv,  vol.  7,  1936-,  pp.  i>l— 55 

(2)  "Ermittlung  der  Eigonschvingzahlen  von  dnrchlaufenden  TrsgSiiA  sit 

feldweise  veranderlicher  Langskraf t" , by  W.  MudTak,xUgenieur-Archiv, 
vol. 7,  1936,  pp.  293-297 

(3)  "Bestimmung  der  Eigenschwingungszahlen  von  durchlauienden  Tragera  und 

Rahmen",  by  W.  Mudrak,  Zeitschrift  fur  Angewandte  Mathematik  and  Mechanik, 
vol.  28,  1948,  pp.  258-263 

(4)  "Vibration  Frequencies  oi  Cy«Linuouo  Beans”,  by  Edward  Saibel,  Journal  of 

the  Aeronautical  Sciences,  vol.  11,  1944*  pp.  88-9C 

(5)  "Free  Vibrations  of  Constrained  Beams",  by  W.  F.  Z.  Lbs  and  E.  Saibel, 

Journal  of  Applied  Mechanics,  vol.  19,  1952,  pp.  471-477 

(6)  "Natural  Modes  of  Vibration  of  Simple  Frames",  by  S.  Bennon,  Journal  of 

the  Franklin  institute,  vol.  243,  1947,  pp.  13-39 

(7)  "Die  Berechnung  der  Drechschvingungen",  by  H.  Holzer,  J.  Springer, 

Berlin,  1921. 

Holzer' s method  is  also  described  in  "Mechanical  Vibrations",  by 
J.  F.  Den  Her tog,  McGraw-Hill  Book  Company,  Inc.,  N.I.,  2nd  Edition, 

1940,  pp.  228-233,  or  in  "Vibration  Analysis",  by  N.  0.  Myklestad, 
McGraw-Hill  Book  Compary.  w.v» . J9AA,  po.  16.5-168 

(8)  nt\  New  Method  of  Calculating  Natural  Modes  of  Uncoupled  Bending  Vibration 

of  Airplane  Wings  and  Other  Types  of  Beams",  by  N.  0.  l^yklestad.  Journal 
of  the  Aeronautical  Sciences,  vol.  11,  1944,  pp.  153-162 

(9)  "A  general  Method  for  Calculating  Critical  Speeds  of  Flexible  Rotors",  by 

h.  H.  Pro hi.  Journal  of  Applied  Mechanics,  vol.  12,  1945,  pp.  142-148 

(10)  "Calculation  of  the  Multiple-Span  Critical  Speeds  of  Flexible  Shafts  by 

Means  of  Pusched-Cnrd  Machines",  by  A.  W.  Rankin,  Journal  of  Applied 
Mechanics,  vol.  13,  1946,  pp.  117-126 

(11)  "Determination  of  the  Natural  Frequencies  of  the.  Bending  Vibrations  of 

Beams",  by  A.  I.  Beilin,  Journal  of  Applied  Mechanics,  vol.  14,  1947, 

pp.  1-6 

(12)  "Analysis  of  Continuous  Frames  by  Distributing  Fixed-Eni  Moments",  by 

Hardy  Cross,  Transactions  A.S.C.E.,  vol.  96,  1932,  pp.  1-10 

(13)  "On  Moment  Balancing  in  Structural  Dynamics",  by  R.  E.  Gaskell,  Quarterly 

of  Applied  Mathematics,  vol.  1,  1943,  pp.  237-249 


172 


K 


I- 


« 


(u) 

(15) 

(16) 

(17) 

(18) 

(3-9) 

(20) 

(21) 

(22) 

(23) 

(24) 

(25) 


"Some  General  Theorems  Relating  "to  Vibrations",  by  Lord  Rayleigh, 
London  Proceedings  Mathematical  Society,  1873,  pp.  357-368}  also 
"The  Theory  of  Sound",  by  Lord  Rayleigh,  Rover  Publications,  1945 
vol s 1.  pp.  15 0-157 


» 

* 


"A  Distribution  Procedure  for  the  Analysis  of  Slabs  Continuous  over 
Flexible  Beams",  by  N.  M.  Nevmark,  University  of  Illinois  Engineering 
Experiment  Station  Bulletin  304,  1938.  p.  25 


"Natural  Frequencies  of  Continuous  Beams  of  Uniform  Span  Length",  by 
R.  Sr  Ayre  and  L.  S,  Jacobsen,  Journal  of  Applied  Mechanics,  vol.  17, 
1950,  pp.  391-395 

"A  Direct  Method  of  Moment  Distribution",  by  T.  Y.  ELn,  Transactions 
A.S.C.E.,  vol.  102,  19.37,  pp.  561-571 

"A  Simple  Method  for  the  Calculation  of  Natural  Frequencies  of  Torsional 
Vibration",  by  F.  P.  Porter,  Transactions  A.S.M.E. , vol.  53,  Paper 
OGP-53-2,  1931,  pp.  17-32. 

A detailed  description  of  Porter’s  method  is  given  also  by  N.  0? 
fyklestad  "Vibration  Analysis",  McGraw-Hill  Book  Company,  Inc.  N.Y., 
1944,  pp.  168-177 

"Fundamentals  of  Vibration  Study",  by  R.  G.  Manley,  John  Wiley  and  Sons, 
Inc.,  N.Y.,  1943 

"Free  and  Forced  Vibrations  of  Composite  Systems",  by  K,  Saibel,  Pro- 
ceedings of  the  Symposium  on  Spectral  Theory  and  Differential  Problems, 
Oklahoma  Agr.  ana  hech.  College,  1951,  pp.  33.3-34 3 

"Stability  of  Structural  Members  Under  Axial  Load",  by  E.  E.  Lundqulsx, 
Technical  Note  No.  617,  National  Advisory  Committee  for  Aeronautics, 
Washington,  D.C.,  October,  1937 


"Die  Beanspruchung  von  Tragverken  durch  schuingende  Lasten",  by  W.  Prager, 
Ingenieur- Archiv , vol;  1,  193 9,  pp-  527—532 

" Be  an  a nruchun  p und  Formanderuncr  von  Siabwerken  bet  onawunffi&nan  Sahwine— 

■ + ' ~ w ~ * ' “ 'W  _ w ■— 

■ungen".  by  S.  Gradstein  and  W.  Prager , Xxgeni eur~Arohiv,  1932,  vol.  2, 

rr  - - ^ - 


"Dynamic  der  Stabwerke",  by  E.  EoheitsOtser  and  V. 

Berlin,  1933 


Pr&gAr,  Julius  Springer, 


"Successive  Elimination  of  Unknowns  in 
J.  B.  Wilbur,  Transactions  A.S.C^E*,  ’ 


Slope  Deflection  Mathod",  by 
102,  1937,  pp,  346-351 


(26)  "Principal  Effects  of  Axial  Load  on  Moment  Distribution  Analysis  of  Rigid 
Structures" , b y 3,  W.  Jsme®,  Technical  Note  No  - 334 , National  Advisory 
Committee  for  Aeronautics,  Washington,  D.C.,  July  1935 


'■a 


\ 


173 

(27)  "Extended  Tables  of  Stiffness  and  Carry-Over  Faclor  for  Structural 

x4  Members  Ubder  Axial  Load",  by  E.  E.  Lundquiet  and  W.  D.  Kroll, 

Wartime  Report  L-255  (ARR  4B  24),  National  Advisory  Committee  for 
Aeronautics,  Washington.  D.C . .February-  1944 

(28)  !'A  Relaxation  Procedure  for  the  Strwas  Analysis  of  a Continuous  Beam- 

Column  Elastically  Restrained  Against  Deflection  and  Rotation  at  the 
Supports '-V  by  Pai  C.  Hu  and  Charles  Libers,  Technical  Note  No.  1150, 
National  Advisory  Committee  for  Aeronautics,  Washington,  D.Co$ 

October  1946 

(29)  "Uber  das  ebene  Knickproblem  dea  Stockwerkrahaens",  by  E.  Chwslla  and 

F.  Jokisch,  Der  Stahlbau,  vol.  14,  p»  33,  1941 

(30)  "Buckling  of  Trusses  and  Rigid  Frames'’,  by  Q.  Winter,  P.  T.  Hsu,  B.  Koo, 

and  Me  H*  Loh,  Cornell  University  Engineering  Experiment  Station  Bull- 
etin 36,  Ithaca,  N.T, 

(31)  "Beams  on  Elastic  Foundation",  by  M.  Hetenyi,  University  of  Michigan 

Press,  Ann  Arbor,  1946 

(32)  "Tables  of  Stiffness  and  Carry-Over  Factor  for  Flat  Rectangular  Plates 

Under  Compression",  by  W.  Do  Kroll,  Wartime  Report  L-398  (ARR  N0.3E27), 
National  Advisory  Committee  for  Aeronautics,  Washington,  D.C.,  Nov.1943 

(33)  "The  Frequency  of  Vibration  of  Rectangular  Isotropic  Plates",  by  R.  F.  S. 

Eearmon,  Journal  cf  Applied  Mechanics,  vol.  19,  1952,  pp.  402-403 

(34)  "On  the  Fundamental  Frequencies  of  Vibration  of  Rigid  Frame* " , by  E,  F» 

Kasur j paper  presented  at  the  First  Midwestern  Conference  on  Solid 
Mechanics,  University  of  Illinois,  April  24,  1953 


! : • ♦Wu-.Afij  a&feb t*-  ,!  &3tt(9£?.k 


M ivi/ tfifU ISHB®I*8!I ? >. tvwv ^ ^ ; sr-<; ‘ ^ ^ **:# wt&R&W  ^ -4 w«tM 


*+L**ki*>i&iZki  $fkit\  'JnVSti  i>  s*£rJE«*wVvi*v-;  3 .' . Hz *;i . <*«i l«e¥« . i.>  ^.-*-j.k  --ki&Juhi  : -i.»  * <1*..  . - -J i,:.' 


COEFFICIENTS  OF  DYNAMIC  FLEXURAL  STIFF 
NESS  Ji  FOR  A UNIFORM  BAR  FIXED  AT 
THE  FAR  END 


CARRY-OVER  FACTOR 


t 

? 


l-Ifl 


Vi 


6.000 


w 


PiG.  6 


COEFFICIENTS  OF  THE  PRODUCT  OF  DYNAMIC 
FLEXURE  - SHEAR  STIFFNESS  AND  DYNAMIC 
FLEXURE- SHEAR  CARRY-OVER  FACTOR 
FOR  A UNIFORM  BAR  FIXED  AT  THE  FAR  END 


1600 


1200 


800 


400 


tE  PRODUCT  OF  DYNAMIC 
AND  DYNAMIC  SHEAR 
•OR  II  FOR  A UNIFORM 
I FAR  END 


DYNAMIC  FLEXURE  - SHEAR  CARRY  - OVER 
FACTOR  q FOR  A UNIFORM  BAR  FIXED 
AT  THE  FAR  END 


SHEAR  CARRY-OVER  FACTOR 


g&asm 


185 


, ♦ 


FIG.  13  DEFLECTED  SHAPE  OF  TWO  ADJACENT  SPANS 

OF  A CONTINUOUS  BEAM  ON  RIGID  SUPPORTS 
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FIG.  14  BEAM  CONSIDERED  IN  EXAMPLE  I 


I 

I! 


F, I,  1.35  E,l, 

0.80  m.  S.20  m. 


1.25  L,- 


L, 


1.35  E,l, 

ni. 

1.50  L, *1 


F»G.  15  BEAM  CONSIDERED  IN  EXAMPLE  2 
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FIG.  IT  FIRST  SIX  NATURAL  VIBRATION  MODES  OF  A UNIFORM,  4- SPAN  CONTI- 
NUOUS BEAM  HINGED  AT  LEFT  END  AND  FIXED  AT  RIGHT  END 
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fWO  NATURAL  VIBRATION  MODES 
CONSIDERED  IN  EXAMPLE  2 
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FIG.  27  NATURAL  VIBRATION  MODES  FOR 

CONSIDERED  IN  EXAMPLE  3 
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FIG.  31  NATURAL  VIBRATION  MODES  FOR  THE  FRAME  CONSIDERED  IN 

EXAMPLE  S 


201 


1 

g ~2-!>20ll) 

A _* 

” « « 


| 

^ ojr 

• I (VI 

k — 


l 91IZ)* 


5 CM 
(O 

« » 

■ m _ m 
12  12 


(S8ZZ  I) 


•o 

*T- 
/ <0 

CVl 

<S> 

k 

<D 

SO 

m 

Cf» 

so 

00 

rO 

1 

<D~ 

+ 

<d“ 

r» 

art 

■«r 

t- 

u 

o> 

i£ 

U> 

1 

(M 

! (2$2fr‘l)' 


j® 
— Cf  lO 


rO  ip 

JP  i* 

CD 

d> 

1 

<xf* 

® 

in 

© 

>n 

<3 

(62*0*2) 
rrrrcj  *A 


- 09  0 - 
6956  5 


(*91!  2) 
- 00  i — 
5191  5 


I OC)7l‘?t 

\ «#  %»W  • VI 

- 19  0 — 
00585 


II 

J 


FIS.  32  ( CONT.'O ) ILLUSTRATIVE  EXAMPLE 
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FIG.  34 


ILLUSTRATIVE  EXAMPLE  6 (CONTINUED) 
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FIG.  42  SYMMETRICAL  SINGLE- BAY  FRAMES  WITH  SIDESWAY 
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ELASTIC  CONSTANTS  FOR  A PANEL  OF 
A SLAB  ON  RIGID  SUPPORTS 


FIG.  45  VARIATION  OF  ROTATION  AMPLITUDE  WITH  O),  FOR  EXAMPLE 


